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FicurE 1. Reconstructions of tree lycopsids, listed 
from left to right: Diaphorodendron scleroticum (modi- 
fied after Wnuk, 1985, fig. 19), Lepidophloios hallii 
(modified after DiMichele & Phillips, 1985, fig. 2), An- 
abathra pulcherrima (redrawn from DiMichele & Phil- 
lips, 1985, fig. 2), Synchysidendron dicentricum (mod- 
ified after Wnuk, 1985, fig. 19), Sigillaria approximata 
(modified after Hirmer, 1927, fig. 284, and Stewart, 
1983, fig. 11.19). Diaphorodendron phillipsii, Lepi- 
dodendron hickii (after Wnuk, 1985, fig. 12, and Tho- 
mas & Watson, 1976, textual description). All x 0.003. 


A. r z 


1989a, b; Donoghue, 1989; see also Hill & Camus, 
1986, on marattialean “‘ferns’’). Conceptually, our 
study owes much to Doyle & Donoghue (1986b) 
in particular, but differs from all the above studies 
in focusing on relationships of taxa within a widely 
accepted order (Class Lycopsida, Order Lepidoden- 
drales) that may lack extant descendants; certainly, 
all of the genera analyzed are extinct. The lycopsids 
are of particular phylogenetic interest as a potential 
sister group of the remainder of the tracheophyte 
clade (Doyle & Donoghue, 1986b, fig. 1; Di- 
Michele & Skog, 1992), 


Our purpose was not merely to unravel the his- 


torical relationships of various arborescent lycop- 
sids. We chose to analyze our data within the now 
well-established framework of evolutionary cladis- 
tics (e.g., Wiley, 1981; Farris, 1983; Funk & 
Brooks, 1990; Wiley et al., 1991) in order to test 
scenarios concerning patterns and underlying 
mechanisms of evolution within the group. In par- 
ticular, we wished to assess preconceived hypoth- 
eses concerning the phylogenetic distributions, 
functional roles, and ecological significance of as- 
pects of growth architecture and reproductive bi- 
ology. 

The structure and content of this paper reflect 
the philosophical framework outlined by Neff (1986) 
and elaborated by Bryant (1989). Bryant (1989, 
fig. 1b) emphasized the creative, deductive nature 
of a priori character analysis and a posteriori phy- 
logenetic interpretation relative to the purely syn- 
thetic, empirical, inductive procedure of tree con- 
struction. Our phylogenetic analysis investigates all 
three of these phases in detail, attempting to exploit 
the main benefit of cladistics: conceptual and meth- 
odological explicitness. 


SELECTION AND PARTITIONING OF WHOLE-PLANT 
SPECIES 


In order to qualify for inclusion in this study, 
plants had to be (1) either members or ee 
outgroups of the Order Lepidodendrales (lycopst 
possessing rhizomorphs, secondary ae 
periderm, ligules, and heterospory: Stewart, 198 4 
(2) anatomically preserved, and (3) known m sub 
ficient detail that all disarticulated component e 
gans could be reconstructed to form a connie 
whole plant (Chaloner, 1986; Bateman & Rothwe , 
1990; Bateman, 1991a); only whole-plant species 
can be thoroughly characterized. In practice, = 
three prerequisites confined our study to the 5 
boniferous of Euramerica (Fig. 3), specifically 4. 
two species of Oxroadia from EEA 
canigenic terrains (Bateman, 1988, 1992) jè a 
species of nine genera (Appendix 1C) spline i Fi 
sylvanian coal swamps (e.g., Hirmer, 1927; ree 
lips, 1979; DiMichele & Phillips, 1985). (We ha” 
deliberately avoided formal reclassification oe 
paper, though several recommendations for z 
nomic revision are outlined in Appendix 1C. P a 
derived from this study segregated Hizemod? 
dron from Lepidodendron sens. str. (Hamna 
DiMichele, 1991) and Synchysidendron a 
Diaphorodendron, as well as erecting pet 
family Diaphorodendraceae (DiMichele & 2 
man, 1992). Our use of the generic name oa 
bathra rather than Paralycopodites follows 
son (1986) and is justified in Appendix 1B). 
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Ficure 2. 
(redrawn from Pigg & Rothwell, 1983a, fig. 1 


1969, fig. 13).—d. Oxroadia gracilis (modified after Bateman, 
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Reconstructions of small-bodied lycopsids.—a. Hizemodendron serratum. —b. Chaloneria cormosa 
).—c. Paurodendron fraipontii (redrawn from Schlanker & Leisman, 


1988, fig. 7.11). Aspects of the Hizemodendron 


peconstruction in particular are speculative, and it is more appropriately described as a restoration (Bateman & 


DiMichele, 1991). a, b = x0.03, c, d = x0.12. 


These conceptual whole-plant species of arbo- 
rescent lycopsid are listed in Table 1, together with 
the bibliogr aphic sources of much of our data. Each 
conceptual whole plant encompasses at least nine 
readily distinguished organs (rootlet, rhizomorph, 
‘tem, branch, leaf, megasporophyll, megaspore, 
microsporophyll, microspore) that are formally 
named, either individually or in aggregates, as or- 
San-species, Table 2 correlates the more important 
of the organ-species that have been awarded Lin- 
nean binomials. Some organ-species binomials en- 
“ompass more than one organ; most of the larger- 

ed whole-plant species encompass five named 
Ln epecies (rootlet/rhizomorph, stem/ branch/ 
eaf, strobilus/ megasporophyll/microsporophyll, 
megaspore, microspore) and the smaller-bodied four 
(rootlet / rhizomorph is not nomenclaturally distin- 
Suished from stem/branch/leaf). Other binomials 
are applied to homologous organs of more than 


one whole-plant species (i.e., form-species sensu 
Bateman & Rothwell, 1990: an organ whose mor- 
phological expression is indistinguishable in two or 
more whole-plant species). Confining our study to 
anatomically preserved material avoided the fur- 
ther complication of correlating the same organs 
of the same whole-plant species in different pres- 
ervation states (e.g., Galtier, 1986; Bateman, 
199 1a). 

Evidence for the reconstruction of these organ- 
species into whole-plant species can be ascribed to 
three main categories (in order of increasing prob- 
ability of correct correlation): association/dissoci- 
ation (co-occurrence in space and time), anatomical 
similarity, and organic connection (this has tradi- 
tionally been regarded as proof of successful re- 
construction, though Bateman & Rothwell (1990) 
argued that at best it constitutes only a strong 
hypothesis). In practice, anatomical similarity is 
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FicurRE 3. Reported time ranges of OTUs (see caption to Table 3 for key to abbreviations). In general, s 
appearances are more reliable than first appearances. Inset shows first recorded occurrences of the genera (x 2 
vegetative remains only, + = reproductive remains only, * = vegetative and reproductive remains). Superr 
denote bibliographic sources; Montagne Noire, south-central France: Rowe & Galtier, 1989, Bateman, 1992 wA 
Meyer-Berthaud, 1984 (2); Burnmouth, northeastern England: Long, 1964 (3), Long, 1968 (4); Glenarbuck, TO 
western Scotland: Smith, 1962 (5), Smith, 1962, 1964 (6); Laggan, Arran, southwestern Scotland: Fry, 195 920, 
Walton, 1935, Scott, 1990 (8); Pettycur, southeastern Scotland: Williamson, 1872, 1893, Scott, 1900, ! ‘ie 
Gordon, 1910, Jongmans, 1930, DiMichele, 1980, Pearson, 1986 (9), Benson, 1908 (10), Scott, 1901, wes ee 
1908 (11); Kingswood, southeastern Scotland: Scott et al., 1986 (12) (see also Scott et al., 1984). Time scale follo 
Leeder (1988). * Flemingites schopfii cone. * F. diversus cone. 


usually employed in conjunction with association/ giate cone (cf. Baxter, 1965; Leisman & ee 
dissociation, and both are used as an adjunct to 1974; Bateman & DiMichele, 1991), and Sigt 
organic connection, especially when attempts to laria “sp. nov.,” which is strictly a Westphalian 
reconstruct whole plants by organic connection are A composite from more than one locality. > 
only partially successful. In the case of our tree- It should be emphasized that the inevitable dis- 
lycopsids, organic connection is particularly diffi- articulation of their constituent individual orga" 
cult to demonstrate between reproductive organs isms into organ-aggregates or single organs renders 
and vegetative axes. It is easier to achieve cor- the reconstructed whole-plant species both con 
relations between different organs within these two ceptual and typological. Each species is conceptu 
main categories; for example, by extracting in situ in the sense that it is a summation of aie 
megaspores and microspores from cones (e.g., Tho- probabilities of correlation of pairs of organs rath 
mas, 1987; Willard, 1989a, b). Additional infor- than an integral, demonstrable fact. It is typologie 
mation acquired since our analyses were performed in the sense that intraspecific variation can s 
suggests a closer relationship between Lepido- be documented among different specimens of | 
phloios johnsonii and L. hallii. Otherwise, the same organ; complete data sets cannot be —_— 
least conclusive reconstructions in our analysis are for specific individuals, in contrast with morp & 
Hizemodendron serratum, where doubts surround metric studies of extant plants (e-g-, Bateman 
its habit and tentatively correlated microsporan- Denholm, 1989a, b). This prevents objective 
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TABLE 1. List of OTUs (1-16) with selected synonyms and bibliographic sources. Oxroadia sp. nov. (2a) was 
not used as an OTU as it did not differ from O. gracilis (2) in the qualitative characters scored. 


Num- 
ber Taxon 


l  Paurodendron fraipontii (Leclercq) Fry 
(syn. Botryopteris fraipontii, Selaginella fraipontii) 
Fry (1954); Phillips & Leisman (1966); Schlanker & Leisman (1969); Rothwell & Erwin (1985) 
2 Oxroadia gracilis Alvin 
Alvin (1965, 1966); Long (1964, 1971, 1986); Bateman (1988, 1992) 
2a  Oxroadia sp. nov. 
Bateman (1988, 1992) 
3 Anabathra pulcherrima Witham 
(syn. Paralycopodites brevifolius, Lepidodendron brevifolium pro parte) 
Felix (1954); Brack (1970); Morey & Morey (1977); DiMichele (1980) 
4 Chaloneria cormosa (Newberry) Pigg & Rothwell 
(syn. Polysporia mirabilis) 
DiMichele et al. (1979); Pigg & Rothwell (1979, 1983a, b, 1985) 
5  Sigillaria approximata Fontaine & White 
Schopf (1941); Delevoryas (1957); Eggert (1972) 
6 Sigillaria sp. nov. 
Brongniart (1836); Benson (1918); Lemoigne (1961) 
7T Synchysidendron sp. nov. 
DiMichele (1979b, 1981); DiMichele & Bateman (1992) 
8 Synchysidendron dicentricum (Felix) DiMichele & Bateman 
(syn. Lepidodendron dicentricum, Diaphorodendron dicentricum) 
Arnold (1960); DiMichele (1979b, 1981, 1985); DiMichele & Bateman (1992) 
9 Diaphorodendron phillipsii DiMichele 
(syn. Lepidodendron phillipsii) 
DiMichele (1981, 1985) 
l0 Diaphorodendron vasculare (Binney) DiMichele 
(syn. Lepidodendron vasculare) 
Carruthers (1869); Hovelacque (1892); Seward (1910); DiMichele (1981, 1985) 
1l Diaphorodendron scleroticum (Pannell) DiMichele 
(syn. Lepidodendron scleroticum) 
Pannell (1942); DiMichele (1981, 1985) 
l2 Hizemodendron serratum (Felix) Bateman & DiMichele 
(syn. Lepidodendron serratum) 
Felix (1952); Baxter (1965); Leisman & Rivers (19 
Michele (1991) 
13 Lepidodendron hickii Watson 
(syn. L. aculeatum pro parte, L. obovatum pro parte 
Scott (1906); Seward (1906); Watson (1907); DiMichele (1983); Willard (1989a) 
l4 Lepidophloios harcourtii (Witham) DiMichele 
(syn. Lepidodendron harcourtii) 
Bertrand (1891); Seward (1899); Zalessky (1912); 
15 Lepidophloios johnsonii (Arnold) DiMichele 
Arnold (1940); DiMichele (1979a); Winston (1988) 
16 Lepidophloios hallii (Evers) DiMichele 
Evers (1951); Felix (1952); Andrews & Murdy (1958); 


(1979a) 
S ieee bee ae mmm 


74); DiMichele (1981, 1983), Bateman & Di- 


Koopmans (1928); Calder (1934) 


Brotzman & Schabiliov (1972); DiMichele 


phenotypy (e.g., Bateman & Denholm, 1989c). 
Ironically, this inability to resolve variation at the 
whole-organism level can be considered advanta- 
geous in a cladistic analysis, In most cases, less 
information is discarded when a paleobotanical spe- 


limitation of individual organisms into species using 
morphological discontinuities and hinders attempts 
to distinguish genetic (and thereby taxonomically 

phylogenetically useful) contributions to phe- 
notype from those caused by ontogeny and eco- 
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Correlations of organ-species that constitute the whole-plant OTUs listed in Table 1. Asterisked whole- 


plant species lack recognized autapomorphies (see caption to Table 3 for explanation of abbreviations). 


OTU Rhizomorph Vegetative axes 
PNFR 
OXGR* Oxroadia gracilis 
OXNS* ——— ricndia wm. ior 
ANPU Stigmaria ficoides Anabathra pulcherrima 
CHCO 
SIAP Stigmaria sp. nov. Sigillaria approximata 
SINS* Stigmaria sp. nov. Sigillaria sp. nov. 
SYNS Stigmaria ficoides | Synchysidendron sp. nov. 
SYDI Stigmaria ficoides Synchysidendron dicentricum 
DIPH Stigmaria ficoides Diaphorodendron phillipsii 
DIVA* Stigmaria ficoides Diaphorodendron vasculare 
DISC Stigmaria ficoides Diaphorodendron scleroticum 
HZSE* —————— Hizemodendron serratum 
LNHI Stigmaria ficoides _ Lepidodendron hickii 
LSHC* Stigmaria ficoides Lepidophloios harcourtii 
LSJO* Stigmaria ficoides Lepidophloios johnsonii 
LSHL Stigmaria ficoides _ Lepidophloios hallii 


cies is reduced to a single unvarying (and thus, by 
definition, typological) data set (the normal pro- 
cedure prior to cladistic analysis) than in compa- 
rable studies of extant species. Also, the paleobot- 
anist is effectively constrained to conceptual 
morphospecies and is therefore spared the trauma 
experienced by neobotanists when selecting an ap- 
propriate species concept (cf. de Queiroz & Don- 
oghue, 1988, 1990; Nelson, 1989: Wheeler & 
Nixon, 1990; Nixon & Wheeler, 1990). Although 
we wished primarily to investigate generic rela- 
tionships, our use of species (of whatever kind) as 
basic operational taxonomic units (OTUs) created 
fewer difficulties when typologically eliminating in- 
tra-OTU variation than the more common ap- 
proach of selecting sets of OTUs from higher (and 
often variable) levels in the taxonomic hierarchy 


(cf. Doyle & Donoghue, 1986b). 


SELECTION, PARTITIONING, AND POLARIZATION OF 
CHARACTERS 


The provisional selection of whole-plant species 
preceded the partitioning of the conceptual organ- 
isms into characters and, subsequently, of the char- 
acters into putatively homologous character states. 
Our chosen characters are listed below; they are 
identified by numbers prefixed by the letter C. 


Cones 


——_—————— Pauredendron fraipontii ———Selaginellites crassicinctus ($) 


Flemingites diversus/schopfii ($) 


Chaloneria cormesa ————————————— errr 


Mazocarpon oedipternum (9+8) 
Mazocarpon schorense/cashii (2+6) 
Achlamydocarpon varius (¢+68) 
Achlamydocarpon varius (?+é) 
Achlamydocarpon varius (¢+é) 
Achlamydocarpon varius (+6) 
Achlamydocarpon varius (?+6) 
Achlamydocarpon sp. nov. (2) 
Lepidostrobus minor (ô) 
Achlamydocarpon takhtajanii (2) 
Lepidostrobus cf. oldhamius (8) 
Lepidocarpon lomaxi (2) 
Lepidostrobus oldhamius (ô) 
Lepidocarpon lomaxi (2) 
Lepidostrobus oldhamius (ô) 
Lepidocarpon lomaxi (2) 
Lepidostrobus oldhamius (é) 


Partitioning a representative, conceptual organ- 
ism of a species into morphological characters 6 
based on the assumption that each character rep- 
resents a discrete, recognizable, and homologous 
feature. This is the most subjective and ultimately 
most influential phase of any cladistic analysis; 1t 
is especially unfortunate that the only sources of 
evidence to support assertions of homology © 
circumstantial; consequently, such assertions sed 
not be conclusively verified (e.g., De Beer, 1971; 
Riedl, 1979; Patterson, 1982; Kaplan, 1984; Roth, 
1984, 1988, 1991; Tomlinson et al., 1984; Neff, 
1986; Ridley, 1986; Bryant, 1989; G. P. Wagner, 
1989). Features that do not vary among the chosen 
OTUs provide no information on their phyogs a 
relationships, though such characters are sis 
in characterizing the entire ingroup (they may, ° 
course, have a greater level of universality than 
the ingroup alone). Continuously variable charac- 
ters can be artificially partitioned into binary - 
multiple character states (e.g., by gap coding (AT 
chie, 1985) or segment coding (Chappill, 1989), 
but we believe that such “‘soft”” characters are ar 
appropriately analyzed by phenetic methods (c! 
Bateman, 1990a; Farris, 1990). We therefore ex- 
cluded continuous (metric) and quasicontinuous 
(meristic) characters from our data matrix, eve? 


though they were the only potential source of unique 
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TABLE 2. Continued. 


Megaspore 


Microspore 


Triangulatisporites triangulatus 
Setispora subpalaeocristata 
Setispora pannosa 

Lagenicula rugosa 
Valvisisporites auritus 
Tuberculatisporites reinschii 
Tuberculatisporites mamillarius 
Cystosporites varius 
Cystosporites varius 
Cystosporites varius 
Cystosporites varius 
Cystosporites varius 
Cystosporites giganteus 


Cystosporites giganteus 
Cystosporites giganteus 
Cystosporites giganteus 


Cystosporites giganteus 


Cirritriradiates annulatus 
Auroraspora cf. asperella ‘A’ 
Auroraspora cf. asperella ‘B’ 
Lycospora orbicula 
Endosporites ornatus 
Crassispora kosankei 
Crassispora sp. nov. 
Granasporites medius 
Granasporites medius 
Granasporites medius 
Granasporites medius 
Granasporites medius 
Lycospora cf. pusilla 


Lycospora pusilla 
Lycospora pellucida 
Lycospora sp. nov. 


Lycospora granulata 


SS See. ee ee ee eS EE EEE EE See 


(and thereby distinguishing) characters for eight of 
A l7 species initially selected for study (Table 


We opted for a uniformly bistate data matrix, 
on the grounds that bistate characters are more 
readily analyzed algorithmically and the distribu- 
tions of character states on the resulting trees are 
more easily interpreted. This decision had three 
Potentially deleterious consequences: 

First, the hierarchy of organs that constitute the 
plants introduced a degree of character duplication; 
we often found it necessary to include a character 
Scoring an organ present or absent (e.g., the rhi- 
zomorph: C8) before partitioning additional char- 
acters (C9-C13) to describe its detailed morphol- 
98y in those OTUs that possess that structure. This 
's to be expected, given the strongly hierarchical 
natur eof morphological and anatomical homologies 
(Riedl, 1979; Fortey & Jefferies, 1982; Wimsatt 
& Schank, 1988; Roth, 1991). We adopted a 
Similarly relaxed attitude to the inclusion of poten- 
tially coupled characters “‘correlated for develop- 
mental-genetic reasons” (Doyle & Donoghue, 
1986b: 338); indeed, we hoped that our analysis 
Would reveal such character correlations, which 
are by no means intuitively obvious a priori. 

cond, even more complicated hierarchies of 
related characters, such as stelar (C14—C19) and 


periderm (C39—C43) anatomy, occur in complex- 
es. In such cases, a broad concept of character 
can be selected, allowing coding in multistate rather 
than bistate format. Such multistate characters tend 
to be especially difficult to polarize satisfactorily 
and may have to be input unordered. If they are 
to be ordered, several methods are available for 
coding and polarizing such characters (e.g., O’Gra- 
dy & Deets, 1987). The preferred options are 
additive binary coding (Brooks, 1984) or nonre- 
dundant linear coding (O'Grady & Deets, 1987; 
O’Grady et al., 1989). Both operate via a hypo- 
thetical tree representing transitions between the 
polarized multiple states of the character in ques- 
tion. By reducing the character states to a set of 
bistate subsets, additive binary coding labels every 
node of the character tree and thus generates a 
large number of narrowly defined operational char- 
acters. Nonredundant linear coding avoids this pro- 
liferation of characters, but at the expense of re- 
taining characters in a multistate format and 
arbitrarily designating within the tree a major axis 
that forms the basis for coding the remaining 
branches (minor axes). Once mixed with other (typ- 
ically mostly bistate) characters, multistate char- 
acters complicate the generation and subsequent 
interpretation of cladograms (e.g., Gensel, 1992). 
We therefore preferred the more primitive but also 
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more intuitive additive binary system, despite the 
risk of eliminating a priori some ambiguities in the 
data. 

Third, polarization of spore ornamentation char- 
acters was rendered especially problematic. They 
tend to be evolutionarily displacive rather than 
additive; there appears to be a developmental con- 
straint on the range of features exhibited by the 
exine of any one species, so that a new type of 
ornamentation supplants rather than supplements 
the previous type. We spurned the option of scoring 
these characters as unordered multistate in order 
to maintain a uniformly polarized binary matrix. 
Thus, we treated each type of spore ornamentation 
as a separate character and assumed a hypothetical 
plesiomorphic spore lacking all features. 

Polarization also presented a more general prob- 
lem. All three OTUs (Oxroadia gracilis—Oxroadia 
sp. nov., Paurodendron fraipontii, Chaloneria 
cormosa) originally screened as potential outgroups 
(Maddison et al., 1984) exhibited some character 
states that we were reluctant to regard as plesio- 
morphic. In Oxroadia and Paurodendron, we ini- 
tially believed that such characters were few and 
almost confined to spore ornamentation. However, 
Chaloneria proved too derived to root the tree 
successfully. Consequently, we constructed a hy- 
pothetical ancestor possessing putatively plesio- 
morphic states for all characters; it largely reflected 
character states shared by Oxroadia and Pauro- 
dendron, though for a few problematic characters 
we elected to screen more distantly related lycop- 
sids for presumed plesiomorphic states. 

As in all cladistic analyses, our recognition of 
alternative character states as plesiomorphic and 
apomorphic preceded tree-building. We restrict 
these terms to character states and use primitive 
and derived to describe the relative positions of 
OTUs on the resulting trees. 

Many of the characters considered for inclusion 
were rejected on the grounds that they were known 
for less than two-thirds of the OTUs. The most 
important examples are the detailed histology of 
the rhizomorph and rootlets (which are well known 
for the form-species Stigmaria ficoides but not for 
the different types of this rhizomorph correlated 
with specific whole-plant species), of the leaves 
(well-documented only for Oxroadia: Bateman, 
1988), and of the gametophyte (described for few 
cone-species; e.g., Galtier, 1964, 1970; Brack, 
1970; Phillips, 1979; Stubblefield & Rothwell, 
1981; Pigg & Rothwell, 1983b), and ultrastruc- 
ture of spore walls (e.g., T. Taylor, 1973; W. 
Taylor, 1990). 


One hundred fifteen bistate characters were 


eventually accepted: 69 are vegetative (C1-C4, 
C8-C72) and 46 are reproductive (C5—C7, C73- 
C115). The large number of characters reflects 
the stringent selection criterion of detailed knowl- 
edge that was applied to potential OTUs; inclusion 
of poorly known OTUs would have increased the 
proportion of missing values in some characters 
sufficiently to warrant exclusion of those characters 
from the data matrix. 


THE CHARACTERS 


Characters are apportioned into 11 categories: 
the overall habit of the organism (A) and ten con- 
stituent organs (B-K). For each category, lists of 
characters and character states are preceded by 
discussions of relevant homologies and descriptive 
terms. 


A. Habit (7 characters) 


We perceived habit as an overall property of 
an organism, expressed as a specific bauplan. We 
describe the group of habits colloquially known a 
trees (C1) as “arboreous,” and use “‘arborescent 
strictly to describe the ability to generate secondary 
tissues (C29). Thus, all 16 OTUs are considered 
arborescent, but only 12 are truly arboreous. We 
reject the frequently used term “secondary wood, 
because wood is by definition secondary. Only ae 
ture woody plants greater than 2 m in overall height 
are termed trees; Chaloneria does not qualify a 
a tree on this criterion, despite possessmg W 
(albeit poorly developed) and an elongate, W- 
branched, upright stem (Pigg & Rothwell, 1983a. 
b). Recumbent OTUs generating limited amounts 
of wood are termed pseudoherbs (Bateman, nee 
1992; DiMichele & Bateman, 1989; Bateman 
DiMichele, 1991). bit 

In lepidodendraleans, the main aspects of ha r 
are stem length (C2), frequency of lateral = 
terminal branching (C3-C4, C6), and the seo 
on the bauplan of reproductive structures (C 
C7). Various combinations of these character yer 
generate several distinct architectures (Figs. l; i 
Four of these characters (C3-C4, C6-C7) reflec 
the mode and timing of branching during meget 
Lepidodendralean stems branched isotomously . 
near-isotomously) only during the final cee 
growth (Walton, 1935; Andrews & Murdy. 199 j 
Eggert, 1961; Wnuk, 1985); terminal penne 
is profuse in most OTUs, but infrequent mM i 
bathra, Sigillaria, and Diaphorodendron, anda 3 
sent in Chaloneria (C3). Strongly anisotomous ap 
cal divisions during stem growth result m pa 
branches that were deciduous (C4) in all O 
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except Diaphorodendron scleroticum, where thick 
bark and wood suggest retention (DiMichele, 1981). 
The anisotomies terminate in cones (C6), either 
individually on short peduncles (Sigillaria) or col- 
lectively on large, repeatedly dichotomous pedun- 
cle systems (C7; Anabathra and Diaphoroden- 
dron). In all other OTUs, anisotomous cone-bearing 
branches were borne only on the terminal, isoto- 
mously branched, determinate crown. We regard 
cauline peduncles and cone-bearing lateral branch- 
es as homologous (DiMichele & Bateman, 1989). 
Chaloneria lacked cones (C5), instead bearing spo- 
rophylls directly on fertile zones of the unbranched 
stem (Pigg & Rothwell, 1983a, b). Consequently, 
it is scored as missing for C6 and C7. 

Polarity decisions for several of the habit char- 
acters were taken with considerable reservations 
(especially C1, C2 and C7). Also, in retrospect, the 
addition of a character representing equality of 
“crown” branching among the pseudoherbs would 
have distinguished Oxroadia (dominantly isoto- 
mous) from Paurodendron and Hizemodendron 
(both dominantly anisotomous). 


l. Nonarboreous (0); arboreous (1). 

2. Stem short (relative to any branches), plant 
recumbent (0); stem tall, plant upright (1). 

3. Dichotomy of trunk apex frequent (0); infre- 
quent or absent (1). 

4. ia lateral branches absent (0); present 
1). 

5. Cone present (0); cone absent (1). 

6. Lateral branches and/or cone peduncles borne 
on dichotomous crown (0); excurrent trunk (1). 

T. Number of cones on lateral branches one (0); 
more than one (1). 


B. Rootstock (6 characters) 


Until recently, the stigmarian rhizomorph (C8) 
Was regarded as arguably the most reliable ubiq- 
ultous character state defining the Lepidodendrales 
e.g., Phillips & DiMichele, 1992). Rothwell (1984) 
and Rothwell & Erwin (1985) suggested that the 
stigmarian rhizomorph is a shoot system modified 
lor rooting; we recognize that it is a shootlike de- 
velopmental system, but prefer to regard it as a 
unique organ reflecting limited developmental op- 
"ons within the arborescent lycopsid bauplan. The 
thizomorph is radially symmetrical (C9), undergoes 
“peated isotomous apical dichotomy (C10), is 
Woody (C11), and emits in helical rhizotaxy rigid 
absorbent rootlets (C12), each containing a single 
monarch vascular strand. Rothwell & Pryor (1990, 
1991) concluded that the tracheary elements of 
such rhizomorphs are derived largely from the 


equivalent of a primary thickening meristem. Such 
branched rhizomorphs were considered radically 
different from other unbranched rootstocks, wheth- 
er radial, as in Paurodendron (Rothwell & Erwin, 
1985), or bilateral, as in Protostigmaria-Lepi- 
dodendropsis (Jennings, 1975; Jennings et al., 
1983), Chaloneria (Pigg & Rothwell, 1983a; Pigg 
& Taylor, 1985), and Isoetes (Karrfalt & Eggert, 
1977 et seq.; Karrfalt, 1984). However, other 
studies of rhizomorphs have rendered less profound 
the distinctions between radial and bilateral sym- 
metry (Karrfalt, 1981, 1984; G. W. Rothwell, 
pers. comm. 1989) and branched and unbranched 
vasculature (Bateman, 1988). In our analysis, bi- 
lateral symmetry (C9) is retained as a character 
state; it is scored as an autapomorphy of Chalo- 
neria and thus does not affect tree topologies. 
Scoring the rhizomorph of Oxroadia as branched 
is an oversimplification; it is extremely compact, 
so that the cortex branches shallowly and less fre- 
quently than the vascular system (Long, 1986; 
Bateman, 1988). 

Stigmarian axes exhibit a suite of anatomical 
character states that, with sufficient study, can be 
shown to parallel those of correlated stem genera. 
Given the current paucity of such studies, we used 
only one such character: the ovoid rootlet gaps 
(C13) found in Chaloneria (Pigg & Rothwell, 
1983a) and Sigillaria approximata (Eggert, 
1972). 


8. Rhizomorph absent (0); present (1). 

9. Rhizomorph symmetry radial (0); bilateral (1). 
10. Rhizomorph branched (0); unbranched (1). 
11. Secondary xylem in rhizomorph absent (0); 

present (1). 
12. Rootlets absent (0); present (1). 
13. Rootlet gap in wood fusiform (0); ovoid (1). 


C. Stele (18 characters) 


The morphology and histology of lepidodendra- 
lean vascular systems, described in detail by pre- 
vious authors, are valuable for distinguishing both 
genera and species (Fig. 4a). Unfortunately, much 
less attention has been paid to determining ho- 
mologies and polarizing these complex characters. 
We recognize six distinct forms of stelar medul- 
lation (C14-C19). All genera but Diaphoroden- 
dron and Synchysidendron are regarded as prim- 
itively protostelic (C14). The protosteles of 
Anabathra, Chaloneria, Sigillaria, Hizemoden- 
dron, Lepidodendron, and Lepidophloios are 
medullated (C15). The core of the stele consists of 
unlignified cells whose dimensions are typical of 
tracheids, suggesting that they are procambial de- 
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FIGURE 4. Morphological and anatomical te 
from fig. 9.44 of Gifford & Foster, 1989: 
secondary xylem, ve = vascular cambium, 


Mii 


outer cortex, pr = periderm, pd = phelloderm, pg = phellogen, pm = phellem, It = leaf traces, lc = leaf perge vt 
b. External morphology of leaf base following leaf loss (modified after fig. 11.3B of Stewart, 1983: ls = feld, Ik = 
= vascular trace, fp = foliar parichnos, Ip = ligule pit, uf = upper field, uk = upper keel, If = lower sonnei 
lower keel, p = plications).—c. External morphology of sporophyll (modified after fig. 11.16C of Stewart, 


fig. 2A of Phillips, 1979: p = pedicel, a = alations, h = heel, Im = lamina, 


lg = ligule, s = sporangium, vt = vase 


_ 1A of 
trace).—d. Leaf anatomy in proximal transverse section (modified after fig. 11.9E of Stewart, 1983, and fig. 14 | 


Reed, 1941: vt = vascular trace, ss = sclerench 


ymatous sheath, m = mesophyll, e = epidermis, ab = lateral 


z , ph 
grooves, ad = median adaxial groove).—e. External morphology of spore (pp = proximal pole, dp = distal pani 


= proximal hemisphere, dh = distal hemisphere, 
cf = contact face). 


rivatives that remained parenchymatous (Walton, 
1935; DiMichele, 1979a, b). In Lepidodendron 
and two of the three Lepidophloios species, ran- 
domly oriented filamentous cells apparently infilled 
a central cavity (C16; DiMichele, 1979a). In Dia- 
phorodendron and Synchysidendron, central pa- 
renchyma cells are distinctly smaller than those of 
associated tracheids, suggesting that these genera 


Feb s é i curvaturs, 
e = equator, ts = triradiate (trilete) suture, | = laesura, c = 


had a true pith sensu Beck et al. (1982). ie 
Synchysidendron species share the spent! 
of a solidly parenchymatous pith (C17), an on 
possesses a histological autapomorphy: pre arily 
cells in Synchysidendron sp. nov., and seco 
thickened cells in S. dicentricum. ~ a 
Deep parenchymatous invaginations in í c20) 
centricum and Synchysidendron sp. nov- 
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are not considered homologous with the shallower 
invaginations of Chaloneria (C21); in the Syn- 
chysidendron species, the parenchymatous wedges 
are raylike, many cells wide and high, and are 
often confluent with the pith parenchyma (Di- 
Michele, 1980, 1981). In contrast, the invagina- 
tions of Chaloneria are smaller and do not reach 
the central parenchymatous area of the stele 
(DiMichele et al., 1979; Pigg & Rothwell, 1983a). 

Protoxylem configuration and leaf trace emis- 
sion comprise an integrally linked complex of char- 
acters (C22-C28; Fig. 4). Protoxylem is exarch 
in all OTUs and, with the exception of Oxroadia 
and Paurodendron, forms a continuous sheath en- 
closing the metaxylem (C22). Concentrations of 
protoxylem observed in transverse sections of axes 
of many genera are often described individually as 
“poles” or “points” and collectively as a “corona.” 
This two-dimensional terminology is misleading; 
protoxylem actually occurs as more-or-less longi- 
tudinal strands that are raised to form ridges in 
“coronate” genera (Bateman, 1988). These pro- 
toxylem strands are longitudinal and linear in 
Oxroadia, Paurodendron, and Lepidodendron, 
but reputedly anastomose in Lepidophloios (C28; 
Bertrand, 1891). 

We have coined new terms for four distinct 
modes of leaf trace emission (C23-C25; Fig. 5). 
Leaf traces departing from a longitudinal proto- 
xylem ridge are termed evaginate. Those of 
Oxroadia and Paurodendron are emitted from a 
single uninterrupted ridge and are termed evagi- 
nate-direct (Fig. 5a). Those of Lepidophloios orig- 
inate within a protoxylem ridge at the point where 
it bifurcates and are termed evaginate-internal 
(C25; F ig. 5d). Most of the genera lack discernable 
protoxylem ridges and are said to emit superficial 
leaf traces (C23; Fig. 5b). In Chaloneria, the trace 
originates from a submarginal position in the stele 
and is associated with shallow parenchymatous in- 
Vaginations (C24, Fig. 5c; DiMichele et al., 1979; 
Pigg & Rothwell, 1983a). We used X-coding (Doyle 
& Donoghue, 1986b) to permit evolution of evag- 
iate-internal and invaginate states directly from 


ior a Se es) ee pei 


— 


Ficure 5, Protoxylem morphology and modes of leaf 
race emission. Transverse section above, longitudinal pro- 
jection of surface of xylem bundle below. Protoxylem in 
ack, spots = leaf traces, solid lines = ridges, dashed 

hes = leaf trace orthostiches where these do not coincide 
with ridges, —a, Evaginate direct: Paurodendron and 
roadia. —b. Superficial: Anabathra, Hizemodendron, 
Lepidodendron. —c. Invaginate: Chaloneria. —d. 
vaginate internal: Lepidophloios. 
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the plesiomorphic evaginate-direct state (thus by- 
passing the superficial state), but to inhibit im- 
probable evolutionary routes that treat the evagi- 
nate-internal and invaginate states as intermediate. 

The stele morphology and leaf trace emission of 
Sigillaria are especially difficult to interpret. The 
undulatory outer margin of the continuous primary 
xylem sheath represents tangential variation in the 
amount of primary xylem produced; protoxylem 
thickness is greatest on the flanks of the undula- 
tions. In contrast with other genera possessing dis- 
cernable protoxylem ridges, sigillarian leaf traces 
appear to originate from the intervening troughs, 
leading to suggestions that each trace may be de- 
rived from both of the adjacent protoxylem ridges 
(e.g., Lemoigne, 1961). However, we were unable 
to confirm their putative bipolar origin and there- 
fore scored Sigillaria leaf traces as superficial. 
The alternative option of recognizing Sigillaria 
traces as bipolar would generate an additional ge- 
nus-level autapomorphy. 

Secondary xylem occurs in the stem and at least 
the more proximal crown branches (if present) of 
all OTUs (C29), but extends into the lateral branch- 
es (C30) only in Diaphorodendron vasculare and 
D. scleroticum. Most OTUs possess homogeneous 
rays composed of small-diameter cells, though het- 
erocellular rays characterize Synchysidendron. 
Consistent nonpreservation prevented character- 
ization of the phloem. 

14-19. Solid protostele (000000); medullated 

protostele (010000); medullated proto- 

stele with filamentous core (011000); si- 

phonostele with mixed pith (100000); si- 

phonostele with solidly parenchymous pith 

including secretory cells (100110); si- 

phonostele with solidly parenchymous pith 

including cells with secondary wall thick- 

enings (100101). 

20. Deep parenchymatous invaginations or 
radial partings absent (0); present (1). 

21. Shallow parenchymatous invaginations 
absent (0); present (1). 

22. Exarch protoxylem sheath discontinuous 

(0); continuous (1). 

Leaf trace origin evaginate, direct (000); 

superficial (100); invaginate (X1 0); evag- 

inate, internal (X01). 

26. Longitudinal ridges of protoxylem strands 
discernable (0); indiscernable (1). 

27. Leaf trace originates from one protoxylem 
strand (0); two protoxylem strands (1). 

28. Anastomoses of protoxylem strands ab- 
sent (0); present (1). 


23-25. 


29. Secondary xylem in trunk absent (0); 
present (1). 

30. Secondary xylem in lateral branches and/ 
or peduncles absent (0); present (1). 

31. Rays homogeneous (0); heterogeneous (1). 


D. Cortex (5 characters) 


All OTUs possess a three-zoned cortex (C32; 
Fig. 4a). It consists of a narrow inner cylinder of 
compact, barrel-shaped parenchyma cells, a thick 
middle cylinder of even thinner-walled, more-or- 
less isodiametric parenchyma cells that often decay 
to leave a cavity, and a broader outer zone of 
thicker walled cells that are longitudinally elongate 
(especially in the central portion of the cylinder of 
tissue) and often grade into sclerenchyma (partic- 
ularly in the smallest diameter cells, close to the 
epidermis). This peripheral sclerenchyma 1s €s- 
pecially well developed in Diaphorodendron scle- 
roticum. Leaf traces passing through the middle 
cortex are ensheathed with cells characteristic of, 
and in continuity with, the inner cortex; they are 
secretory in several OTUs (C33) and adaxially 
concentrated in Synchysidendron sp. nov. (C34). 
In Sigillaria sp. nov., Synchysidendron, and Dia- 
phorodendron, each leaf trace is surrounded by a 
broad cylinder of thinner-walled cells when passing 
through the outer cortex (C35), increasing appar- 
ent cellular heterogeneity. No attempt was made 
to divide variations in the angle of passage of the 
leaf traces through the cortex into discrete char- 
acter states. Vertically elongate cavities at the cor- 
tex-periderm transition (C36) characterize D. phil- 
lipsii. 

32. 
33. 


Outer cortex two-zoned (0); three-zoned (1). 

Intracortical leaf-trace sheaths not secretory 

(0); secretory (1). 

34. Intracortical leaf-trace sheaths circumferen- 
tial (0); adaxial (1). 

35. Thin-walled parenchyma surrounding leaf trace 
in thick-walled outer cortex absent (0); present 
(1). 4 

36. Cavities at outer cortex-periderm transition 

absent (0); present (1). 


E. Periderm (14 characters) 


Periderm occurs in all of the OTUs (C37), though 
its distribution (and thereby its protective and sup’ 
portive function) is extremely restricted in the bent: 
plans of Oxroadia and Paurodendron. asa 
references to periderm as “‘secondary cortex ~ 
an anatomical non sequitur; cortex is 4 ya 
an axis (between the stele and the epidermis) rat 


Volume 79, Number 3 
1992 


than a specific tissue type. The periderm of most 
genera cannot be differentiated histologically into 
phellem and phelloderm. The exceptions are Dia- 
phorodendron and Synchysidendron, where clear 
bizonation is strong evidence for a bifacial cambial 
layer that produced much greater quantities of 
centripetal phelloderm than centrifugal phellem 
(C38; DiMichele, 1981). Additional, indirect evi- 
dence indicates bifaciality in Anabathra (Di- 
Michele, 1980) and Lepidophloios (DiMichele, 
1979a). Given the determinate growth and early 
onset of peridermal cambial function that are ev- 
ident in lepidodendraleans, the cambial layer may 
have fully differentiated or been active only near 
the apices of stems and rhizomorph axes. Either 
phenomenon would inhibit preservation of the cam- 
bial layer per se. It is also possible that the cen- 
tripetal and centrifugal products of cell division are 
sufficiently similar to prevent recognition of a po- 
tentially fully differentiated cambium. 

Periderm is the most abundant tissue type pro- 
duced by arboreous lycopsids and often occurs as 
abundant disarticulated fragments in coal-ball as- 
semblages (Phillips & DiMichele, 1981; DiMichele 
et al., 1986). Fortunately, the detailed anatomy 
and histology of the periderm (C39-C43) allow 
identification of five groups of OTUs. Primitive 
genera (Oxroadia, Paurodendron, Anabathra, 
Chaloneria) have a uniform periderm, which is 
modified to include bands of resinous cells in Sig- 
illaria (C40). Diaphorodendron and Synchysi- 
dendron possess bifacial periderm (C41); 
Synchysidendron is distinguished from Diapho- 
rodendron by its uniform (C43) rather than banded 
(C42) phelloderm. Lepidodendron and Lepido- 
Phloios have two- or three-zoned periderm (C39) 
and are X-coded for resinous cell clusters to sup- 
Press improbable evolutionary routes that attain 
histological modifications via the acquisition of zo- 
nation, 

Glandular periderm histology (C47) is shared by 
Lepidodendron and two of the three Lepidophloi- 
°S Species, and resinous sacs (C48) occur in Sig- 
illaria, Lepidodendron, and Lepidophloios. The 
Passage of leaf traces and infrafoliar parichnos 
strands through the periderm also distinguish OTUs. 
Prominent leaf traces, more-or-less perpendicular 
to the length of the axis and surrounded by thin- 
walled parenchyma (C49), are an autapomorphy 
of Anabathra. Similarly, well-developed infrafoliar 
Parichnos strands (C50) characterize Sigillaria 
periderm. 

Many lepidodendraleans retained leaf cushions 
®n the stem surface as periderm production ex- 
panded axial girth. We recognize three retention 
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mechanisms (C44—C46). The first two are, by def- 
inition, mutually exclusive; interareas exhibit either 
a plastic response and expand (C44; Sigillaria and 
Synchysidendron) or a brittle response and fissure 
(C45; Diaphorodendron) (DiMichele, 1981). In 
contrast, Lepidodendron hickii accommodates 
growth by expansion of cells beneath the cushion 
(C46; DiMichele, 1983). This character state could 
have replaced interarea expansions or fissuring, or 
it could have arisen directly from the plesiomorphic 
state; it is X-coded to allow any of these options. 
Evidence for the interarea expansion of Lepido- 
dendron is confined to compression fossils (Tho- 
mas, 1966). 


37. Periderm in stems absent (0); present (1). 

38. Phellem and phelloderm not histologically 

differentiable (0); histologically differen- 

tiable (1). 

Cellular composition of periderm uniform 

(00000); cells form two or three distinct 

zones (1X000); bands of resinous cell clus- 

ters (01000); periderm bifacial, alternat- 

ing bands of thick- and thin-walled cells 

in phelloderm (00110); periderm bifacial, 

phelloderm + uniform (00101). 

Leaf cushion retention mechanism absent 

(000); tangential interarea expansion 

(100); interarea fissuring (010); subcush- 

ion cellular expansion (XX 1). 

47. Periderm nonglandular (0); glandular (1). 

48. Periderm nonresinous (0); at least par- 
tially resinous (1). 

49. Leaf traces in periderm obscure (0); prom- 
inent (1). 

50. Infrafoliar parichnos strands in periderm 
absent (0); present (1). 


39-43. 


44-46. 


F. Leaf bases (15 characters) 


Lepidodendralean leaf base characters are well 
reviewed by DiMichele (1979a, b, 1981, 1983) 
for anatomically preserved species and by Thomas 
(1970b, 1977, 1978) and Thomas & Meyen (1984) 
for adpressed species. 

Lycopsid leaves attenuate bilaterally close to the 
stem, where they are consequently most readily 
detached. The area proximal to the constriction 
persists as a symmetrical structure raised above 
the surface of the axis and is termed a leaf base; 
aggregates of leaf bases preserve the phyllotaxy of 
the axis after leaf loss. In our more primitive OTUs 
(Oxroadia, Paurodendron, Anabathra, Chalo- 
neria), leaf bases are small, ellipsoid in transverse 
section, widely spaced on the axial surface, and 
are fully transitional into the leaf lamina (C51). In 
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the more derived OTUs, the leaf base is more 
elaborate and only a portion emits the leaf; it is 
then termed a leaf cushion (Fig. 4b). Such cushions 
typically exhibit a simple angular outline in tan- 
gential section: hexagonal in Sigillaria, and dia- 
mond-shaped in Diaphorodendron, Synchysiden- 
dron, Hizemodendron, Lepidodendron, and 
Lepidophloios (this character was not coded). 
Cushions are further elaborated by the develop- 
ment of raised upper (C54) and lower (C55) keels, 
and by division into upper and lower fields that are 
usually separated by a lateral line (C58) and are 
independently plicate (C56—C57). Sigillaria is ple- 
siomorphic for all five characters, while Diapho- 
rodendron and Synchysidendron are apomorphic 
for all five. Lepidophloios only possesses keels. 
Hizemodendron possesses plications but lacks an 
upper keel, while Lepidodendron lacks upper field 
plications but possesses a lateral line. 

The leaf cushion can be regarded as an elabo- 
rated leaf base and thus as fundamentally foliar 
and appendicular in nature. This interpretation is 
supported by the leaflike structural and positional 
attributes of leaf cushions; rejection of the foliar 
nature of leaf cushions would require their rec- 
ognition as developmentally distinct organs of a 
kind unknown in other plants. Even the relatively 
simple cushions of Sigillaria are helically ar- 
ranged, closely packed but discrete, well-defined 
features. Leaf bases of all the arborescent lycopsids 
analyzed bear ligules, which are ancestrally foliar 
in the class (Bonamo et al., 1988). All of these 
characters are features of appendicular organs pro- 
duced laterally to the apical meristem through the 
formation of primordia. In order to be axial rather 
than foliar, leaf cushions would have to be epider- 
mal/subepidermal elaborations that enlarged be- 
low, and concurrently with, leaf primordia as the 
leaves expanded. There is no evidence in any ar- 
borescent lycopsid (with or without cushions) for 
such an unintuitive developmental mechanism 
wherein the leaf determines differentiation of the 
axis. Moreover, stomata occur on leaves and leaf 
bases, but are absent from axes; they are exclu- 
sively foliar. 

Together, the above qualitative characters sep- 
arate all of the cushion-bearing genera studied here, 
though morphometric quantification is necessary 
to separate species within each genus (Thomas, 
1970b; DiMichele, 1983; Chaloner & Meyer-Ber- 
thaud, 1983; DiMichele et al., 1984; Wnuk, 1985). 
Although we have deliberately excluded metric 
characters from our analysis, to avoid arbitrary 
division of such characters into states, we made 
an exception for the length: width ratio of the leaf 


cushions (C52—-C53). Since the early nineteenth 
century, the plesiomorphic condition of greater leaf 
cushion length than width has been crucial for 
delimiting Lepidodendron sens. lat. (DiMichele, 
1983). Despite early knowledge that awarding pri- 
macy to this character resulted in the lumping of 
morphologically and anatomically dissimilar species 
(e.g., Scott, 1908; Seward, 1910), only recently 
has Lepidodendron sens. lat. been disaggregated 
into the morphologically distinct segregates Ana- 
bathra (DiMichele, 1980), Diaphorodendron sens. 
str. (DiMichele, 1983), Synchysidendron (Di- 
Michele & Bateman, 1992), Hizemodendron 
(Bateman & DiMichele, 1991), Lepidodendron 
sens. str., and Lepidophloios (DiMichele, 1979, 
1983). Horizontally (i.e., tangentially) elongate leaf 
cushions are, however, a valid generic autapo- 
morphy of Lepidophloios, together with radial 
elongation (C64—C65). Obscure evolutionary re- 
lationships of the arched and the perpendicular 
states of radial elongation among different Lepi- 
dophloios species necessitated X-coding. 

All OTUs possess ligules (C59), and most recess 
the ligule in the cavity that communicates with the 
adaxial surface of the leaf base via a deep pit (C60). 
The plesiomorphic exceptions are Paurodendron. 
where the ligule is fully exposed (Phillips & Leis- 
man, 1966), and Anabathra, where it is afforded 
some protection by the leaf cushions (DiMichele, 
1980). Foliar parichnos (C61) occur in all OTUs 
but Oxroadia and Paurodendron (Bateman (1988) 
was unable to substantiate Long’s (1986) tentative 
identification of parichnos in Oxroadia). In con- 
trast, infrafoliar parichnos (C62) are confined to 
Lepidodendron (Weiss, 1907; DiMichele, 1983) 
and Lepidophloios (W. A. DiMichele, pera 

The presence of a leaf cushion (C51) is positive y 
correlated with deciduousness (C63) in all OTUs 
but Hizemodendron, where retention is probably 
secondary (Bateman & DiMichele, 1991). We de- 
liberately avoid describing leaf loss in lepidoden- 
draleans as abscission (cf. Chaloner & Mopani 
thaud, 1983). Despite the consistent absence 0 
leaf laminae from axes (including almost all twigs) 
an abscission layer has not been detected in y 
OTU at any stage of development. We sist 
that lycopsid leaves atrophy and are mechamca 4 
removed. This occurred more readily in trees SUC 
as Sigillaria than in smaller erect plants such as 
Chaloneria and recumbent pseudoherbs such = 
Oxroadia, due to the higher basal stresses 1m™p°" 
by the long microphylls of the arboreous nee 
Fracture occurs where the lamina constricts © 
form the leaf base; its position relative to the et 
and lamina is more consistent in OTUs possessing 
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leaf cushions, where a sharp structural (and, pre- 
sumably, physiological) boundary is represented as 
a leaf scar. The scar is situated within the leaf 
rather than at the leaf-axis junction. Deciduous 
lateral branches may also have been shed by this 
mechanism; according to Jonker (1976), triangular 
marks below ulodendroid scars indicate that the 
branches were torn away at their junctions with 
the stem. However, R. A. Gastaldo (pers. comm. 
1990) argued that most specimens lack these fea- 
tures, which may be taphonomic overprints. 


51. Leaf is outgrowth of entire leaf base (0); 
portion of leaf base (1). 

52. Length:width ratio of cushions on stems 
and large branches >1:1 (0); <1:1 (1). 

53. Length: width ratio of cushions on small 
branches and twigs >1:1 (0); <1:1 (1). 

54. Upper keel absent (0); present (1). 

55. Lower keel absent (0); present (1). 

56. Upper field nonplicate (0); plicate (1). 

57. Lower field nonplicate (0); plicate (1). 

58. Lateral line separating upper and lower 

fields absent (0); present (1). 

Ligule absent (00); superficial or in shal- 

low depression (10); in deep cavity with 

narrow neck (11). 

61. Foliar parichnos absent (0); present (1). 

62. Infrafoliar parichnos absent (0); present 
(1). 

63. Consistent basal limit to leaf atrophy ab- 

sent (0); present (1). 

Leaf cushion not radially elongate (00); 

elongate, strongly arched (1X); elongate, 

+ perpendicular to axis (X1). 


59-60. 


64-65. 


G. Leaves (7 characters) 


Leaves are the organs most frequently neglected 
when attempts are made to reconstruct lepidoden- 
draleans. Admittedly, their deciduousness and con- 
“equent absence from the upper axes hinders organ 
correlation by organic connection, but even more 
ĉonsistent disarticulation has not prevented indirect 
correlation of cone species with the vegetative axes 
that bore them. Failure to correlate leaves with 
their parent plants is unfortunate, as Graham (1935) 
and Reed (1941) demonstrated the wide range of 
Potentially phylogenetically valuable characters 
Present in isolated lepidodendralean leaves, and 

teman (1988) recorded many characters (in- 
cluding details of the cuticle, epidermis, and sto- 
mata) of leaves attached to Oxroadia axes. We 

arded most of these characters for this analysis, 
because they would have contained unacceptably 
"B® proportions of missing values. 


The presence of two vascular strands (C67) is 
autapomorphic for Sigillaria, which also shares 
V-shaped strands (C68) with Chaloneria. Dorsi- 
ventrally flattened strands (C69) occur in Lepi- 
dodendron, Lepidophloios, and Diaphoroden- 
dron scleroticum; ambiguous evolutionary pathways 
from plesiomorphic terete strands necessitated 
X-coding. All OTUs whose leaves possess nonterete 
strands also possess pronounced lateral abaxial 
grooves (furrows) containing stomata (C70; Fig. 
4d); these are supplemented with a median adaxial 
groove in Sigillaria (C71). The vascular strands 
of most of the more apomorphic OTUs are sur- 
rounded by a sclerenchymatous sheath (C72), 
though this is absent from Hizemodendron. 

Angle of leaf attachment (C66) refers only to 
the angle subtended by the basal portion of the 
mature lamina relative to the distal portion of the 
axis, thus avoiding the effect of recurvation in 
OTUs such as Oxroadia. This character distin- 
guishes genera with hispid, generally short leaves 
(Paurodendron, Anabathra, Chaloneria, Lepi- 
dodendron), but can result from one of several 
developmental mechanisms and is therefore prone 
to homoplasy. 


66. Angle of leaf attachment relative to axial 
apex + horizontal (0); acute (1). 

67. Number of vascular strands per leaf one 

(0); two (1). 

Transverse section of vascular strand te- 

rete (00); dorsiventrally flattened (1X); 

V-shaped (X1). 

70. Lateral abaxial grooves absent (0); present 
at least near base (1). 

71. Median adaxial groove absent (0); present 
at least near base (1). 

72. Sheath of sclerenchyma around trace ab- 
sent (0); present (1). 


68-69. 


H. Cones (4 characters) 


Intensive study of lepidodendralean reproduc- 
tive structures has generated thorough reviews of 
both anatomically preserved (Arber, 1914; Bal- 
bach, 1967; Brack, 1970; Hanes, 1975; Phillips, 
1979; Brack-Hanes & Thomas, 1983; Willard, 
1989a) and adpressed (Lesquereux, 1880; Kid- 
ston, 1923-1925; Willard, 1989b) organ-species. 
Reproductive characters played important roles in 
the delimitation of the genera; not surprisingly, 
many are genus-level autapomorphies. Moreover, 
many of the traits are functionally linked and can 
be used to define reproductive strategies in the 
same manner that vegetative morphology defines 


growth habits. 
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All of the lycopsids included in this analysis are 
heterosporous (C75). Oxroadia, Paurodendron, 
and Anabathra have primitively bisporangiate 
strobili (C76), with microsporangia concentrated 
toward the cone apex. The fertile zones of Chal- 
oneria (arguably a derived condition) are similarly 
bisporangiate; all other OTUs bore monosporan- 
giate cones. 

Characters 73 and 74 describe the relationships 
to the parent stem of the lateral cone-bearing axes, 
whether peduncles or branches (Chaloneria is un- 
branched and was scored as missing for C73). 
Lateral branches are subtended by stelar gaps (C73) 
in Diaphorodendron, Synchysidendron, Lepi- 
dodendron, and Lepidophloios, and are medul- 
lated (C74) in Sigillaria, Diaphorodendron, Lep- 
idodendron, and Lepidophloios. 


73. Stelar vascular gap associated with departure 
of peduncle or lateral branch absent (0); pres- 
ent (1). 

74. Pith in trace of peduncle or lateral branch 
absent (0); present (1). 

75. Plants homosporous (0); heterosporous (1). 

76. Cones/fertile zones bisporangiate (0); mono- 
sporangiate (1), 


I. Sporophylls and sporangia (15 characters) 


Literature review suggests that the terms de- 
scribing most components of the sporangium-spo- 
rophyll complex have become standardized (Fig. 
4c). The sporophyll is divided into a proximal por- 
tion (“pedicel”), perpendicular to the cone axis, 
and a distal portion (“distal lamina”), parallel to 
the cone axis and oriented toward the cone apex. 
The adaxial surface of the pedicel bears the spo- 
rangium and (immediately distal to the sporangium) 
the ligule. The pedicel is triangular in median trans- 
verse section and attenuates abaxially, to a struc- 
ture that has been termed a keel if sufficiently 
prominent (Phillips, 1979), and laterally, to struc- 
tures that have received various names. “Lateral 
laminae” is used most commonly; alternatives are 
“lateral extensions” (e.g., Meyen, 1987), “wings” 
(e.g., Arber, 1914), “flanges” (e.g., Sporne, 1975), 
and “alations”’ (e.g., Phillips, 1979). Some authors 
(e.g., Phillips, 1979) have distinguished the most 
developed state of this character (long and enrolled) 
as “integuments,”” by analogy with the true seeds 
of “spermatophytes.” The apically directed distal 
lamina is much less three-dimensional and usually 
appears as a shallow “V” in transverse section. 
An antapically directed extension from the right- 


angled junction of the pedicel and distal lamina is 
termed the heel. 


Comparing sporophylls with sterile microphylls, 
we believe that the pedicel is homologous with the 
leaf base (including the leaf cushion), and the spo- 
rophyll distal lamina is homologous with the leaf 
lamina. We suggest that the qualifier “distal” should 
be abandoned for the sporophyll lamina (there is 
no proximal lamina), and that the lateral exten- 
sions of the pedicel should be termed alations, ir- 
respective of size and orientations (use of the term 
“‘integuments”’ for extensive enrolled alations mis- 
leadingly implies homology with the integuments 
of true seeds). 

Little attention is paid in the literature to angle 
of pedicel attachment relative to the cone axis 
(C77), which may be prone to ontogenetic change 
as an aid to passive spore dispersal (Bateman, 1988). 
Thus, our identification of Oxroadia as autapo- 
morphic for obtuse sporangia is tentative. With 
this exception, all sporophyll and sporangium char- 
acters (C78-C91) are scored as plesiomorphie for 
the four most primitive OTUs possessing bisporan- 
giate cones (Oxroadia, Paurodendron, Anabath- 
ra, Chaloneria), which differ quantitatively rather 
than qualitatively. A good example is the number 
of megaspores per megasporangium, which also 
separates species of the same genus (e.g-, Ana- 
bathra: Felix, 1954; Brack, 1970; see also Ap- 
pendix 1B). It is tempting to distinguish qualita- 
tively between megasporangia containing four 
spores, derived from one spore mother cell, and 
those containing more than four spores, deriv 
from more than one spore mother cell. However, 
each condition characterizes one of the two spectes 
of Oxroadia (Bateman, 1988), and spore counts 
are complicated by frequent and apparently ran- 
dom abortions. 

Differences among the remaining (monosporal 
giate) genera focus on megasporangiate cones an 
reflect their shared transition in the nature of the 
dispersal unit from isolated megaspores to a meg- 
asporangium-sporophyll complex (C81). The ap% 
morphic state of this character encapsulates a han 
spectrum of morphologies (elaborated in C80 an 
C82-C90) that may reflect parallel (i.e-, hapai 
plastic) responses to similar selective AA 
all monosporangiate genera but Sigillaria, z 
evolutionary trend results in reduction to 4 smg 
functional megaspore (C90) that germinates W! ye 
the sporangium (C80; Phillips, 1979). Probably č 
an aid to dispersal and/or protection, these chang 
are accompanied by lateral expansion of the arp 
to form alations (C82—C84). These are coded 
short and horizontal in Sigillaria, Diaphoroder 
dron, and Synchysidendron, short and erect 
Hizemodendron and Lepidodendron, and long 
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erect (typically enrolled) in Lepidophloios. These 
characterizations require further revision; the lat- 
eral margins of the pedicel can be proportionately 
longer in Anabathra (e.g., fig. 9 of Brack, 1970) 
and more erect in Diaphorodendron and Syn- 
chysidendron (e.g., pl. 8.4 of Phillips, 1979) than 
those of Lepidodendron (e.g., pl. 5.4 of Phillips, 
1979). Moreover, shortness may not be homolo- 
gous between vertical and horizontal alations; hence, 
X-coding was used to allow evolution of short, erect 
alations from either absence of alations or short, 
horizontal alations, and to suppress evolution of 
short, horizontal alations from short, vertical ala- 
tions. 

Megasporangia of Diaphorodendron and Syn- 
chysidendron are strongly dorsiventrally flattened 
(C85-C86) and dehisce proximally (C87-C89), 
while those of Hizemodendron, Lepidodendron, 
and Lepidophloios are strongly bilaterally flat- 
tened and dehisce distally. Sigillaria approximata 
undergoes indehiscent fragmentation (C89), pre- 
sumably an apomorphic character state. We have 
used X-coding to allow its evolution by one step 
from any of the three dehiscence mechanisms. 
Sigillaria approximata possesses another autapo- 
morphy, the enclosure of megaspores with paren- 
chyma (C91), 

Heterocellular sporangium walls (C79) charac- 
terize Diaphorodendron and Synchysidendron. 
Only Lepidodendron possesses a multiseriate spo- 
rangium wall; together with greater sporangium 
size (a quantitative character and therefore not 
coded), this distinguishes Lepidodendron sporan- 
§1um-sporophyll complexes from the otherwise 
identical equivalents of Hizemodendron. 


77. Angle of sporophyll attachment relative 
to cone apex + horizontal (0); obtuse (1). 

78. Sporangium wall uniseriate (0); multise- 
riate (1). 

79. Sporangium wall homocellular (0); het- 
erocellular (1). 

80. Megaspores shed from sporangium prior 
to germination (0); megaspores germinate 
within sporangium (1). 

81. Dispersal unit megaspore (0); megaspo- 

rangium-sporophyll complex (1). 

Alations of megasporophyll pedicel absent 

(000); short, horizontal (100); short, sub- 

erect (X10); long, erect (O11). 

Transverse section of megasporangium + 

cireular (00); strongly bilaterally flattened 

(10); strongly dorsiventrally flattened (01). 

Megasporangium dehisces longitudinally 


82-84, 


85-86. 


87-89. 
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(000); distally (100); proximally (010); 
indehiscent fragmentation (XX 1). 

90. Functional megaspores per megasporan- 
gium more than one (0); one (1). 

91. Parenchyma enclosing megaspores absent 


(0); present (1). 


J. Megaspores (10 characters) 


Figure 4e summarizes terms describing the “‘ge- 
ography” of the exteriors of lycopsid spores. 

Morphological and ultrastructural studies of ly- 
copsid megaspores preserved in situ in cones have 
been undertaken since the earliest applications of 
palynology to biostratigraphic and paleoecological 
problems (Schopf, 1938; Bochenski, 1939; Brack, 
1970; Taylor, 1990; see also Bartram, 1987). 

Polarity decisions for laesural (C94—C96) and 
equatorial (C92-C93) characters were problem- 
atic; they are generally poorly developed in dis- 
tantly related lycopsids, but better developed and 
more complex in putatively more closely related 
outgroups (e.g., Selaginella: Stanier, 1965; Tryon 
& Lugardon, 1978; Minaki, 1984) and in the more 
primitive ingroup members (Oxroadia, Pauroden- 
dron, Chaloneria) than in the more derived OTUs. 
Prominent laesural expansions characterize 
Oxroadia, where they are fimbriate and do not 
extend beyond the curvaturae (C95: Alvin, 1965, 
1966; Bateman, 1988), and Paurodendron, where 
they are plicate and extend to the equatorial flange 
(C96: Guennel, 1952). The laesurae of Diapho- 
rodendron and Synchysidendron megaspores are 
gulate (C94); the spongy, trilobate proximal massa 
is a key taxonomic character. Equatorial expan- 
sions provide autapomorphies for Paurodendron, 
in the form of a perisporial plicate flange (C93: 
Guennel, 1952), and Chaloneria, in the form of 
auriculae (ear-shaped expansions of the exine) op- 
posite laesural rays (C92: Pigg & Rothwell, 1983b). 

Most OTUs lack dispersed proximal and distal 
ornamentation. Contact-face ornamentation is con- 
fined to Oxroadia (sparse, robust, buttressed spines: 
C97), Paurodendron (reticulate: C98), and Chal- 
oneria (rugose: C99). The megaspore of Sigillaria 
sp. nov. could not be scored for this character, but 
its distal surface clearly bears short spines (C100: 
Benson, 1918; Pigg, 1983). Large and more com- 
plex buttressed spines typify Oxroadia mega- 
spores. Paurodendron megaspores bear a striking 
distal reticulum (C101). 


92-93. Equatorial ornamentation absent (00); 
auriculate (10); flanged (01). 
94-96. Laesural ornamentation absent (000); 
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gulate (100); fimbriate (01X); plicate 
(OX1). 
97-99. Contact-face ornamentation absent 
(000); echinate (100); reticulate (010); 
rugose (001). 
100-101. Distal ornamentation absent (00); echi- 
nate (10); reticulate (01). 


K. Microspores (14 characters) 


Several lycopsid microspores commonly en- 
countered in dispersed miospore assemblages (sen- 
su Chaloner, 1970) have been correlated with source 
cones, both anatomically preserved (Brack, 1970; 
Courvoisier & Phillips, 1975; Willard, 1989a) and 
compressed (Thomas, 1970a, 1987; Willard, 
1989b). Classification of lycopsid microspores has 
focused on equatorial elaboration and general sur- 
face ornamentation. The only exception in our list 
of characters, strongly raised laesurae (C105), oc- 
curs in Chaloneria and Lepidophloios harcourtii. 

In Chaloneria, separation of the sexine and 
nexine layers distal to the contact faces has gen- 
erated a pseudosaccus (C102: Brack & Taylor, 
1972). All of the monosporangiate-coned genera 
exhibit some form of equatorial elaboration (C106- 
C108). Anabathra, Hizemodendron, Lepidoden- 
dron, and Lepidophloios microspores possess a 
thickened equatorial band (cingulum: C106); in 
Hizemodendron, Lepidodendron, and some Lep- 
idophloios species, this is supplemented with an 
external membranous flange (zona: C108). Sigil- 
laria, Diaphorodendron, and Synchysidendron 
microspores bear a crassitude equatorial thickening 
(C107) that appears structurally distinct from a 
cingulum. Cingula of some OTUs are further elab- 
orated; those of Hizemodendron, Lepidodendron, 
and Lepidophloios harcourtii are bizonate (C103) 
and those of Anabathra are distally ornamented 
(C104), 

Characterization of microspore general surface 
morphology is increasingly dependent on the great- 
er resolution of scanning electron microscope (SEM) 
studies relative to light microscopy (LM). As in the 
megaspores, ornamentation is described separately 
for contact faces (C109-C1 11) and the distal hemi- 
sphere (C112-C1 15), though Lepidophloios john- 
sonii is insufficiently known to be scored. Contact 
faces of Hizemodendron, Lepidodendron, and 
Lepidophloios microspores are granulate (C109; 
Leisman & Rivers, 1974; Willard, 1989a), those 
of Diaphorodendron and Synchysidendron grade 
into a more foveolate texture (C110; Courvoisier 
& Phillips, 1975), and those of Oxroadia and 
Paurodendron are echinate (C111; Schlanker & 


Leisman, 1969; Bateman, 1988). Several OTUs 
lack contact-face ornamentation, but only Chalo- 
neria and Lepidophloios harcourtii lack distal 
ornamentation. In Oxroadia and Paurodendron, 
the echinate contact faces are paralleled by the 
distal ornamentation (C115). Diaphorodendron and 
Synchysidendron microspores are papillate (C114), 
those of Sigillaria are characterized by a mixture 
of spines and cones (C113), and those of Hize- 
modendron, Lepidodendron, and Lepidophloios 
hallii bear dense grana (C112). 


102. 
103-104. 


Pseudosaccus absent (0); present (1). 
Equatorial expansion absent (00); 
unornamented bizonate cingulum com- 
plex (10); distally ornamented cingu- 
lum complex (01). 

Laesurae subdued (0); strongly raised 
(1). 

Equatorial ornamentation absent (000); 
cingulum complex (100); crassitude 
(010); zona (001). 

Contact-face ornamentation absent 
(000); granulate (100); granulo-foveo- 
late (010); echinate (001). 

Distal ornamentation absent (0000); 
densely granulate (1000); echino-co- 
nate (0100); papillate (0010); echinate 
(0001). 


105. 


106-108. 
109-111. 


112-115. 


THE DATA MATRIX 


After all 115 characters had been scored, two 
whole-plant species (Oxroadia gracilis and 
Oxroadia sp. nov.) possessed identical data sets, 
demonstrating that they differ in quantitative but 
not qualitative characters. The duplicate data set 
provided no useful information and was therefore 
omitted, reducing the original 17 whole-plant spe 
cies to 16 OTUs that, in our opinion, represent 
10 genera. This gave a ratio of characters: OTUs 
of 115:16 (7.2: 1). 

The resulting 1,840-byte data matrix (Table 3, 
excluding HYAN) contained 94 missing values that 
each reflected one of four factors. Two of these 
factors result from ignorance of whether the O 
possesses the feature or what state the feature 
exhibits, one from unwillingness to specify the pre- 
cursor state of the character (the X-coding pt” 
cedure discussed in detail by Doyle & Donoghue 
1986b: 344-350; see also Appendix 2), and one 
from absence of the relevant feature, which there- 
fore cannot be scored. The first three categorie 
indicate varying degrees of ignorance vanci 
the nature of the character and can be replat 
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by a posteriori optimization, substituting values most 
parsimonious with the cladogram in question. The 
fourth type of missing value (coded # in Table 3) 
is also replaced during optimization, but here the 
result is merely an operational necessity. It does 
not generate a potentially biologically meaningful 
hypothesis; there is no obvious meaning in scoring 
a character state as present for an OTU that com- 
pletely lacks the feature in question (for example, 
in Fig. 6 Chaloneria is optimized as possessing 
peduncle/branch gaps (C73) but actually lacks both 
types of organ). This situation is more frequent in 
data matrices that include OTUs of highly diver- 
gent morphology, where the probability of obtain- 
ing fully compatible sets of homologous features is 
less. The most gap-ridden (12% missing values) of 
our data sets is Chaloneria cormosa, where 11 of 
the 14 missing values reflect lack of the coded 
feature (i.e., absence rather than ignorance). Com- 
pared with OTU selection, we were less rigorous 
when excluding gap-ridden characters from the 
matrix; the worst examples (C34, C40, C69) pos- 
sess 31% missing values, most of these X-coded 
to represent ambiguity of precursor states. Nev- 
ertheless, the overall proportion of missing values 
n our data matrix (5%) compares favorably with 
those of other studies (e.g., 24% in Table 2 of 
Doyle & Donoghue, 1986b). 

We believe that every cladistic data matrix should 
routinely carry character: OTU ratios to indicate 
the average strength of empirical support for nodes 
and percentage missing values to indicate the com- 
pleteness of the data matrix, just as the resulting 
trees now routinely carry consistency indices to 
Summarize levels of homoplasy (Kluge & Farris, 
1969; Brooks et al., 1986). 


OPTIMIZATION, CHARACTER STATES, AND 
MISSING VALUES 


The distribution of an apomorphic character 
state among all the OTUs can be assigned to one 
of three categories: 

The apomorphic state is confined to a single 
OTU and is thus an autapomorphy at the least 
inclusive (species) level in the taxonomic hierarchy. 
It is important to note that the autapomorphic 
condition is a relative concept; a synapomorphy 
(shared derived character) at the species level in 
Sur analysis can be an autapomorphy (unshared 
derived character) at the more inclusive level of 
genus. Autapomorphies at the least inclusive level 
analyzed distinguish OTUs, but are phylogeneti- 
cally uninformative, Consequently, they are omit- 
ted from algorithmic analysis to avoid artificially 


increasing length and decreasing perceived levels 
of homoplasy by including characters that are, by 
definition, nonhomoplastic (e.g., Brooks et al., 1986; 
Kluge, 1989; Sanderson & Donoghue, 1989). 

2. The apomorphic state is ubiquitous among 
the OTUs (including the outgroups if used), thereby 
justifying their status as a potential clade. Although 
such character states are conventionally described 
as “basal synapomorphies” or “invariant charac- 
ters” (e.g., Sanderson & Donoghue, 1989), we 
prefer to coin the more parsimonious term “‘holapo- 
morphy.” It could be argued that such character 
states are merely plesiomorphies, but they cannot 
be defined as such in the absence of an equivalent 
apomorphic state. Our definition of the term syn- 
apomorphy is also unconventional in implicitly ex- 
cluding holapomorphic character states. Holapo- 
morphy, like autapomorphy, is a relative concept; 
addition to the suite of taxa analyzed of an OTU 
lacking the apomorphic character state transforms 
a formerly ubiquitous holapomorphy into a non- 
ubiquitous synapomorphy. Holapomorphies also re- 
semble autapomorphies in being phylogenetically 
uninformative within the confines of a particular 
data matrix. Hence, like autapomorphies, holapo- 
morphies should be (and in this study were) omitted 
from tree length calculations. 

3. When the apomorphic state occurs in more 
than one but less than all of the OTUs, it is deemed 
phylogenetically informative and included in the 
algorithmic computation of tree length. Most such 
character states are synapomorphic, although cat- 
egory (3) also encompasses homoplasies (these are 
generally regarded as refuting the initial hypothesis 
of homology between the plesiomorphic and apo- 
morphic states: e.g., Wiley, 1981; Funk & Brooks, 
1990). The perception of synapomorphies as ho- 
moplastic (i.e., as parallelisms and/or reversals) or 
nonhomoplastic (evolving only once and persisting 
throughout the derived portion of the clade) is the 
least stable aspect of cladistic analysis, since these 
conditions are a property only of the interaction 
of a specific tree topology with a specific optimi- 
zation algorithm (see below). Also, character states 
can be both homoplastic and partially synapo- 
morphic; for example, a particular state may be a 
meaningful synapomorphy of the OTUs forming 
clade A+B but be represented as a parallelism in 
their non-sister OTU D (J. I. Davis, pers. comm. 
1990). Surprisingly, the terminologically rich dis- 
cipline of cladistics does not appear to have gen- 
erated unique terms to describe the important (al- 
beit ad hoc) distinction between synapomorphies 
sens. str. (i.e., holapomorphies excluded) that are 
homoplastic and those that are not; only the latter 
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TABLE 3. Cladistic data matrix. Operational taxonomic units: Hypothetical ancestor (HYAN), Paurodendron 
fraipontii (PNFR), Oxroadia gracilis-sp. nov. (OXGR), Anabathra (Paralycopodites) pulcherrima (ANPU), Chal- 
oneria cormosa (CHCO), Sigillaria approximata (SIAP), Sigillaria sp. nov. (SINS), Synchysidendron (Diapho- 
rodendron) sp. nov. (SYNS), Synchysidendron dicentricum (SYD), Diaphorodendron phillipsii (DIPH), Diapho- 
rodendron vasculare (DIVA), Diaphorodendron scleroticum (DISC), Hizemodendron (Lepidodendron) serratum 
(HZSE), Lepidodendron hickii (LNHI), Lepidophloios harcourtii (LSHC), Lepidophloios johnsonii (LSJO), Lepi- 
dophloios hallii (LSHL). Known values: primitive (0), derived (1); missing values (coded 9 in PAUP; functional states 
in preferred most parsimonious tree are indicated by subscripts): not known whether OTU possesses relevant feature 
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reliably characterize an entire monophyletic por- _ teriori procedure performed using one of a range 
tion of a clade. of algorithms that are designed to apply specific 
We are also surprised at the paucity of literature logical precepts to specifying the nature (e.g-+ 7 
concerning optimization, as it proved tobe a crucial versal vs. parallelism) and location of each char- 
aspect of our analysis. Optimization is an a pos- acter transition on a tree whose topology and | 
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TABLE 3. Continued. 
(—), OTU possesses relevant feature but character state unknown (?), OTU lacks relevant feature (#), precursor 


state ambiguous (X; X-coded sensu Doyle & Donoghue, 1986b). Functional level of generality of derived character 
state in preferred most parsimonious tree: + = species level autapomorphy (character state restricted to a single 
OTU; these provide no information on historical relationships of species or genera), * = genus level autapomorphy 
(character state restricted to a single genus but occurring in more than one species; these provide no information on 
historical relationships of genera), @ = holapomorphy (= basal synapomorphy; these provide no information on 
historical relationships of species or genera). 


Stele (14-31) 
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Periderm (cont’d.) Leaf bases (51-65) 
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have already been fixed. In theory, optimization had the advantage of maximizing the number of 


cannot alter the length of a tree (though see below optimization algorithms that could legitimately be 
for a critical reappraisal of this conventional wis- applied (Swofford, 1985: 3.10). We tested all five 
dom). optimization algorithms available in PAUP 2.4: 


Rooting the tree using the hypothetical ancestor ACCTRAN maximizes reversals, as in practice does 
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TABLE 3. Continued. 


Leaf bases (cont’d.) 
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Leaves (66-72) 
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Megaspores (92-101) 


eS ee en 
$ + S + + + + + + 
HYAN 0 0 0 0 0 0 0 0 0 0 
PNFR 0 l 0 Xe l 0 l 0 0 l 
OXGR 0 0 0 l X, l 0 0 1 0 
ANPU 0 0 0 0 0 0 0 l 0 0 
CHCO l 0 0 0 0 0 0 0 0 0 
SIAP 0 0 0 0 0 0 0 0 0 0 
SINS 0 0 0 0 0 2, % 2, l 0 
SYNS 0 0 l 0 0 0 0 0 0 0 
SYDI 0 0 1 0 0 0 0 0 0 0 
DIVA 0 0 1 0 0 0 0 0 0 0 
DISC 0 0 l 0 0 0 0 0 0 0 
HZSE 0 0 0 0 0 0 0 0 0 0 
LNHI 0 0 0 0 0 0 0 0 0 0 
LSJO 0 0 0 0 0 0 0 0 0 0 
LSHL 0 0 0 0 0 0 0 0 0 0 


FARRIS; in contrast, DELTRAN maximizes par- 
allelisms. MINF concentrates character state tran- 
sitions toward the terminal branches, MINRES 
concentrates all possible transitions toward the root, 
then concentrates the remainder toward the ter- 
minal branches. For our preferred most parsimo- 
nious tree (PMPT), FARRIS yielded similar (though 
not identical) results to ACCTRAN, and MINRES 
yielded similar results to DELTRAN and, to a lesser 


degree, MINF. We preferred MINRES and MINF, 
as they minimized perceived homoplasy in the lower 
(genus-level and above) branches that were ° 
greatest interest to us; they also tended to ye 
the greatest number of intuitively satisfactory =A 
lections when (1) choosing between reversals a 
parallelisms, and (2) substituting 0 or 1 for missing 
values of specific characters (see below). eer 
vations expressed by Swofford (1985: 3.9-3.1 


= a SS 
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TABLE 3. Continued. 
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Microspores (102-115) 
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concerning MINRES encouraged our consistent result are different for every tree of every analysis; 
use of MINF. for example, those presented as subscripts to miss- 

Dur ing optimization, each missing value (coded ing values in our data matrix (Table 3) refer only 
9 in the PAUP data matrix) is replaced with O or to the preferred MPT of analysis A following ap- 

in accordance with (1) the topology of the tree plication of the MINF optimization algorithm. Re- 
and (2) the intended effect on patterns of character placement of missing values is achieved parsimo- 
State transformation of the chosen optimization al- niously in accordance with the topology of the tree, 
gorithm (this is an a priori opportunity to modify so that further homoplasy will not be introduced. 
evolutionary interpretations). The substitutions that Hence, we presume that this procedure artificially 
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reduces levels of homoplasy relative to those that 
would have been determined from the same matrix 
if complete (i.e., gap-free). 

Doyle & Donoghue (1986b: 352) imposed dual 
origin (parallelism) where a structure was simpli- 
fied, and an origin and a loss (reversal) where a 
structure became more complex, on the assumption 
that “‘it is easier to reduce or lose a complex struc- 
ture than to elaborate one from a simple structure.” 
We did not consider this generalization sufficiently 
reliable to warrant a posteriori modification of our 
algorithmically optimized character state distribu- 
tions, though we did emphasize complex, appar- 
ently conservative characters in subsequent evo- 
lutionary interpretations. 

According to Bateman (in prep.), when all OTUs 
are scored as possessing either the plesiomorphic 
or apomorphic state for a character, the assignment 
of that character to one of these three ostensibly 
exclusively categories (autapomorphy, holapomor- 
phy, synapomorphy sens. str.) is definitive; its con- 
dition is fixed for that data matrix. However, this 
principle does not apply if the character column 
contains at least one missing value (9), when the 
results of optimization determine the status of the 
character as informative or uninformative. For ex- 
ample, in our data matrix, a character scored as 
fifteen ls and one 9 (e.g., C29) will be holapo- 
morphic if the 9 is replaced a posteriori with a 1, 
but synapomorphic if the 9 is replaced with a 0 
(in practice, this will occur only if the OTU scored 
9 is placed at the base of the cladogram). Similarly, 
a character scored as fourteen Os, one 1, and one 
9 (e.g., C94, C98) will be autapomorphic if the 9 
is replaced with a 0 but nonautapomorphic if the 
9 is replaced with a 1. Hence, there is a need for 
the concepts of transient autapomorphy, transient 
holapomorphy, and transient synapomorphy, to 
accommodate characters that contain missing val- 
ues. These concepts are relative, even within a 
single data matrix; they characterize only a single 
combination of a specific topology and a specific 
optimization algorithm. In our novel terminology, 
definitive and transient autapomorphies together 
constitute operational autapomorphies (likewise for 
holapomorphies and synapomorphies). 

In contrast with definitive autapomorphies and 
holapomorphies, transient autapomorphies and hol- 
apomorphies cannot be screened out of an analysis 
a priori. Consequently, they contribute to tree length 
as calculated during cladistic analysis; they intro- 
duce spurious extra steps in an unpredictable man- 
ner, often generating alternative trees that are 
incorrectly considered of equal length by the par- 
simony algorithm. For example, a data matrix gen- 
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erates two topologies, A and B, both X steps in 
length. A posteriori screening for transient autapo- 
morphies/holapomorphies reveals four in topology 
A and two in topology B. A is then preferred over 
B as its true length is X — 4, relative to X — 2 
in topology B. Unfortunately, topology C, per- 
ceived by the algorithm as X + 1 steps long but 
containing six transient autapomorphies/holapo- 
morphies, has a true length of X — 5 steps and is 
actually the most parsimonious tree. Thus, the true 
lengths of trees generated from a data matrix con- 
taining missing values can only be calculated a 
posteriori. Algorithmically determined tree lengths 
are unreliable, and trees other than those that are 
ostensibly the most parsimonious must also be 
screened individually via apomorphy lists for spu- 
rious additional steps. The alternative option of 
omitting all potential transient autapomorphies and 
holapomorphies a priori (appendix 1 of Sanderson 
& Donoghue, 1989) deleteriously discards poten- 
tial synapomorphies merely because their frequen- 
cies among the OTUs approach zero or unity. 


PARSIMONY ANALYSIS 


Cladograms were generated from the data ma- 
trix using Version 2.4 of PAUP (Swofford, 1985), 
which employs unrestricted parsimony via the 
Wagner method (Kluge & Farris, 1969; Farris, 
1970; Felsenstein, 1982; Swofford & Maddison, 
1987; Wiley et al., 1991). Some of the compu- 
tational difficulties encountered by us and discussed 
below have been at least partially surmounted by 
more recent software (see Appendix 2). Despite 
the long run-times incurred on our IBM-PS2/# 80, 
the branch-and-bound option (a modification of the 
algorithm devised by Hendy & Penny, 1982) was 
used routinely to obtain the definitive shortest trees. 

Once character scoring had been finalized, five 
different configurations of the data matrix were 
analyzed: 

A. All 16 OTUs and all 115 characters inched 
This basic analysis provided a yardstick by whic 
to measure the remaining analyses. 

B. All OTUs included, “habit” characters (Cl- 
C7) excluded. We wished to reassess the phylogeny 
without these characters for two reasons. as. 
they describe the most generalized aspects of plan 
morphology and are thus most prone to ace 
Second, we wished to map the distribution of E 
copsid bauplans onto a phylogeny constructed i? 
dependently of such characters (see also gee 
& DiMichele, 1991; Phillips & DiMichele, 199°) 

C. All characters included, but Chaloneria cor 
mosa excluded. Survey of tree topologies rom 
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analyses A and B demonstrated that Chaloneria 
is the most unstable OTU; it is supported by the 
least robust node, characterized only by homo- 
plastic and autapomorphic characters, and pos- 
sesses more autapomorphies (six) than any other 
OTU (Fig. 6). This most awkward OTU was omitted 
in order to determine how the topologies of the 
more parsimonious trees would be altered and 
whether homoplasy would decrease significantly. 

D. All OTUs included, analysis restricted to veg- 
etative characters (C1—C4, C8-C72). 

E. All OTUs included, analysis restricted to re- 
productive characters (C5-C7, C73-—C1 15). Anal- 
yses D and E were performed to determine the 
relative contributions of vegetative and reproduc- 
tive characters to the whole-plant phylogeny (cf. 
Bateman & DiMichele, 1991), and to assess the 
likely accuracy of phylogenies based on the partial 
plants that constitute most paleobotanical “‘spe- 
cies.” Data matrices for organ-species phylogenies 
are much easier to construct than those for whole- 
plant phylogenies, given the difficulty of correlating 
vegetative and reproductive organs. 

After some experimentation, we developed an 
analytical routine that was applied to each of our 
main groups of analyses (A-D above; analysis E 
generated more equally most parsimonious trees 
than PAUP 2.4 can store). In each case, an initial 
run used the BANDB command to find all equally 
Most parsimonious topologies by branch-and-bound, 
and the combination of the OPT=MINF optimi- 
zation algorithm and APOLIST print command to 
identify the putative location and direction of each 
character state transition. Having thus determined 
the length of the shortest tree(s) (Lan), we then 
reanalyzed the data matrix by replacing the BANDB 
command with BB = ‘X’, where ‘X’ was one step 
longer than the shortest tree (i.e., L,,,,,). This 
Second run found and saved trees of lengths Ln 
and Lans; in order to determine the number of 
topologies of length L,,,,,, the total number of 
topologies found at BB = L,,, was subtracted from 
the total number of topologies found at BB = L,,..1- 

Procedure was repeated up to lengths of about 
Lant: (depending on the particular submatrix under 
Scrutiny). Tree number increases more-or-less ex- 
Penentially with increase in length; the maximum 
Capacity of PAUP 2.4 to store trees (N = 100) is 
Soon exceeded, so that it becomes untenable to 
routinely scrutinize topologies much longer than 


We found such scrutiny desirable for two rea- 
Sons. First, we wished to know how many genus- 
level topologies occurred at each length, as opposed 
to species-(OTU-)level topologies routinely detect- 
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ed by the tree-building algorithm. Second, we wished 
to use lists of apomorphies following optimization 
to assess each tree for transient holapomorphies 
and autapomorphies, so that they could be sub- 
tracted to obtain its true length (determined entirely 
by synapomorphies sens. str.). In practice, this 
time-consuming screening procedure was not ap- 
plied to trees longer than L,,,,, (a new algorithm 
is required for this purpose), rendering optimization 
and the retention of apomorphy lists redundant 
from L,,,4. onward (analyses A, B, D) or Lunts 
onward (analysis C). 

Having surveyed all optimally (L,,,,) and subop- 
timally (Lain <3) parsimonious trees, we focused on 
particular trees of interest, including all most par- 
simonious trees (MPT, i.e., tree of length Lun) 
These were reprinted with APOLIST (a list of node 
by node character state transitions) and CHGLIST 
(a list of character by character state transitions) 
for full interpretation (e.g., Figs. 6, 7). Topologies 
of potential interest longer than those routinely 
surveyed (i.e., longer than ca. L,,,,,; Fig. 7) were 
specified using the TOPOLOGY command in “user 
tree” mode (Swofford, 1985: 2.20-2.22). 

In summary, our cladistic analyses were exper- 
imental sensu Doyle & Donoghue (1986b). The 
basic philosophy of this approach was well sum- 
marized by Johnson (1982) and Bryant (1989: 
218): ‘‘Parsimony determines the order by which 
viable hypotheses should be tested; one starts with 
the simplest” (our italics). Alternative hypotheses 
are then considered in order of increasing com- 
plexity until a self-imposed threshold is reached. 
In contrast, nonexperimental cladistic studies both 
start and finish with the simplest. 


SPECIES-LEVEL RELATIONSHIPS 


Although this experimental cladistic study was 
aimed primarily at elucidating genus-level relation- 
ships, we chose to perform our analyses at the 
species level. This decision partly reflected a sub- 
sidiary interest in species-leve] relationships, but 
was taken primarily because species-level OTUs 
provide a test of the presumed monophyly of gen- 
era. Genera can then be re-delimited if necessary. 
The following discussion is based primarily on anal- 
ysis A, but also applies to analyses B-D (analysis 
E produced an untenably large number of equally 
most parsimonious trees). 

Only four of the ten genera in the data matrix 
are represented by more than one OTU (Table 1): 
Sigillaria (two species), Diaphorodendron (three 
species), Synchysidendron (two species), and Lep- 
idophloios (three species, but see Appendix 1D). 
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i05 Paurodendron fraipontii 


i 
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. ¢ “a i (19) 
56 
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= aia 1 66 
rr E eee = 58 Lepidodendron hickii 
a (46 78) 
oa, Lepidophloios harcourtl 
26 39 47 48 53 62 68 70 105 (*) 
25 28 52 -57 64 84 z Lepidophloios johnsonii 
z (*) 
-47 -103 Lepidophloios hallii 
(65) 


FIGURE 6. Preferred most parsimonious cladogram for analysis A (all OTUs and all characters). Holapomorphies 


are placed in square brackets below the hypothetical ancestor, and autapomorphies are placed in parentheses 


below 


OTUs; all character state transitions on terminal branches are therefore homoplastic. Characters that experience 
reversal are underlined (with a minus sign where a reversal occurs), parallelisms are overlined. Asterisks indicate 


absence of character state transitions at specific nodes. 


Each of these four genera proved very robust (i.e., 
dismantling each genus in any way resulted in much 
longer trees), but relationships of the species within 
at least two of the genera are less clear (the two 
species of Sigillaria and Synchysidendron re- 
spectively do not allow multiple topologies). 

The three-taxon problem presented by Diapho- 
rodendron phillipsii, D. vasculare, and D. sclero- 


ticum is very poorly resolved (Fig. 6), resting ©": 
tirely on the nonhomoplastic synapomorphy i 
secondary xylem in lateral branches (C30) that 
unites D. scleroticum and D. vasculare (Fig. 8d). 
Treating C30 as a synapomorphy of the clade 

as a secondary loss in D. phillipsii (a tena?’ 
hypothesis more consistent with stratigraphic i 
dence; Fig. 3) costs only one extra step and cok 


—_—— 
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jas. Paurodendron fraipontii 
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FIGURE 7. Fully annotated cladogram for analysis A showing alternative generic topo 


logies of potential interest. 


t differs from Figure 6 in that (1) Paurodendron is a sister group of Oxroadia (cost = two steps), (2) Anabathra 
and Chaloneria form a clade (cost = nil). (3) Sigillaria and Diaphorodendron form a clade (cost = one step), and 
(4) Hizemodendron and Lepidodendron sens. str. form a clade (cost = five steps). Character notation follows Figure 6. 


lapses the relationship into an uninformative tri- 
chotomy (Fig. 8e). This polychotomy was a 
Persistent cause of trivially multiple topologies in 
our analyses (see Appendix 2). 

The phylogenetic relationship of the three pu- 
vars Lepidophloios species is obscured by exten- 
sive homoplasy that is compounded by ambiguities 
caused by missing values for some characters. Lep- 


idophloios hallii and L. harcourtii are subtended 
only by homoplasies, while L. johnsonii is not 
subtended by any characters (it even lacks quali- 
tative autapomorphies; Fig. 6). We obtained three 
equally most parsimonious solutions to this three- 
taxon problem (Fig. 8a-c). The first (Fig. 8a) links 
L. johnsonii and L. hallii by the homoplastic syn- 
apomorphy of a filamentous core to the protostele 


(65) (a) (x) 


JO HC HL 
(*) (*) (65) 


25 28 52 -57 64 84 


PH VA SC 
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VA PH sc 
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(C16). The second (Fig. 8b) links L. johnsonii and 
L. harcourtii by the homoplastic synapomorphy 
of a glandular periderm (C47) and by the ostensibly 
nonhomoplastic synapomorphy of arched leaf cush- 
ions (C64). Unfortunately, this character is scored 
as missing for L. hallii as a result of X-coding 
(Table 3) and can therefore be treated as either 
apomorphic or plesiomorphic during optimization. 
As a result of inconsistent replacement during op- 
timization, the apomorphic state is depicted as 
characterizing all three Lepidophloios species in 
Figure 8a but only two species in Figure 8b. The 
third topology (Fig. 8c) treats C16 as a loss in L. 
harcourtii and C64 as present in all three species, 
which consequently collapse to a trichotomy. 
Given that all three solutions are equally most 
parsimonious, every topology that differed in the 
positions of OTUs other than Lepidophloios spe- 
cies was repeated three times by the tree-building 
algorithm to accommodate the multiple solutions 
to the Lepidophloios problem (hence our division 
by three of the algorithmically determined numbers 
of species-level topologies to yield the smaller, more 
meaningful values listed in Table 4). The arrange- 
ment of Lepidophloios species shown in Figure 8a 
best fits their reported sequence of relative ap- 
pearance in the fossil record (Fig. 3); on the basis 
of this weak evidence, it was preferred when sê- 
lecting the trees shown in Figures 6 and 7. 


EXPERIMENTAL CLADISTICS: A SURVEY OF 
GENUS-LEVEL TOPOLOGIES 


In all analyses, Synchysidendron and Diapho- 
rodendron sens. str. consistently remained united 
as a monophyletic group and are not distinguished 
in Figures 9 and 10. Also, the following discussion 
occasionally refers to Synchysidendron as derived 
relative to Diaphorodendron and Lepidophloios 
as derived relative to Lepidodendron. As these 
pairs of genera are sister groups (Fig. 6), these 
assertions of derivation are subjectively impo 
by us, based on comparison of the number an 
inferred evolutionary significance of the character 
state transitions supporting each genus. 


eS eee 


n= 


FicurE 8. Poorly resolved relationships between T 
cies of the same genus (see caption to Table 3 for key 5 
abbreviations). a-c show three equally parsimonious ( 
step) solutions to the three taxon problem present ú 
the Lepidophloios species. d and e show two — n- 
the three taxon problem presented by the Diaphor : 
dron species (d = 7 steps, e = 8 steps). Notation ca 
follows Figure 6, though lines below character n 
emphasizing reversals are omitted. 


TABLE 4. Comparison of cladistic parameters for five permutations (A-E) of the data matrix in Table 3. A = all taxa and all characters included; B = all taxa 
included, habit characters (C1—C7) omitted; C = Chaloneria cormosa omitted, all characters included (the total number of characters declines from 115 to 109 due 
to the elimination of Chaloneria’s six autapomorphies; also one former synapomorphy becomes an autapomorphy of Paurodendron); D = vegetative characters (C1- 
C4, C8-C72) only; E = reproductive characters (C5-C7, C73-C115) only. For numbers of trees, the initial figure is the number of species-level topologies of that 
length, the figure in parentheses is the number of genus-level topologies of that length, and the figure in square brackets is the cumulative total of genus-level tree 
topologies (ND = not determined). Figures for topologies longer than L,,,,,, may be inaccurate, as it was not feasible to screen by eye large numbers of trees for 
transient holapomorphies and transient autapomorphies (see “Optimization, character states, and missing values”). Strictly, each species-level value should be multiplied 
by three to reflect three equally parsimonious solutions to the relationships of the three species of Lepidophioios (this triangular relationship became unreliable in all 
analyses at Lana). The figure for analysis E remains greatly inflated by repeated dichotomies. Note that Hizemodendron serratum was treated as a new genus distinct 


from Lepidodendron sens. str. when calculating the genus-level values. Segregating Synchysidendron from Diaphorodendron sens. str. did not affect calculation of 
tree numbers, as the two genera consistently behaved as a monophyletic group. 


Analysis 
Parameter A B C D E 
Number of OTUs (excluding HYAN) 16 16 15 16 16 
Total number of characters 115 108 109 69 46 
Number of operational holapomorphies 8 8 8 7 1 
Number of operational species-level autapomor- 
phies 27 25 22 il 16 
Number of informative characters remaining 
(synapomorphies sens. str.) 80 75 79 51 29 
Minimum length (L,,,; steps) 128 115 121 we 43 
Consistency index at Lu» 0.625 0.652 0.653 0.662 0.674 
Number of trees of minimum and | a: 3) [3] 2 ANS] 1 ely) tty 2, 2y [2] 270 (ND) [ND] 
near-minimum length ks 14 (11) [11] ae Cs es T a. (2y fi 13 (11) [11] — 
| a 31 (19) [22] 28 (18) [18] 5 (4) [5] 44 (ND) [ND] _ 
Lss 62 (ND) [ND] 70 (ND) [ND] 14 (10) [11] 129 (ND) [ND] — 
| ae 134 (ND) [ND] 133 (ND) [ND] 32 (ND) [ND] — — 
Lains = = 51 (ND) [ND] = = 
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Generic topologies of analysis A at L,,, (a-c) and Lanti (d-k) (see caption to Table 3 for key t 
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abbreviations). For analysis B (habit characters omitted), a and b occur at La, and c-f, h at L,.,41 (the age nes 
between Chaloneria and Anabathra became unresolved in c). For analysis C (Chaloneria omitted), only sing kis wala 
topologies occur at L,,, (that seen in a-c less Chaloneria) and at L,,,,,, (that seen in e-g). 1 is a see on 
tree for analysis A at L,..,., analysis B at L,,,,,, and analysis C at L,,,,, (less Chaloneria). Synchyside fe eee 
Diaphorodendron are not distinguished in Figures a-g as they consistently behave as sister groups. All c 


information is omitted. 


Analysis A (all OTUs and all characters includ. 
ed). Analysis A yielded three equally most parsi- 
monious trees (MPTs; Fig. 9a-c) that unite 
Oxroadia and Paurodendron as a basal clade. The 
preferred most parsimonious tree (PMPT) depicts 
the remaining genera as a perfectly nested se- 
quence of increasingly apomorphic OTUs (Figs. 6, 
9a). The three MPTs differ in the position of Chal- 
oneria, which occurs immediately above Anabath- 
ra in the preferred MPT (Fig. 9a), immediately 
below Anabathra in the second MPT (Fig. 9c), 


and is united with Anabathra to form a mone 
phyletic clade in the third MPT (Fig. 9b). af 
Decreasing the level of parsimony by adding D 
step (L,,,.,) yields another eight generic topologi h 
(Fig. 9d-k). Three of these unite Sigillaria al 
Diaphorodendron-Synchysidendron (Fig. ah 
otherwise repeating the three possible relations 
between Chaloneria and Anabathra seen n ee 
9a~c. Four others dissociate the Oxroadia- 
rodendron clade, so that each arises directly Fie. 
the major axis (Fig. 9d, h-j); three of these (P£ 
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h-j) also involve changes in the relative positions 
of Anabathra and Chaloneria. It is then equally 
parsimonious to have Oxroadia (Fig. 9d, j) or 
Paurodendron (Fig. 9h, i) as the basal OTU of the 
tree. The eighth topology unites Sigillaria and 
Chaloneria, linking both to Anabathra to form a 
substantially different topology (Fig. 9k). 

At Lantz Some topologies unite Chaloneria and 
Sigillaria, others allow the exchange of Hizemo- 
dendron and Sigillaria across Diaphorodendron— 
Synchysidendron (cf. Fig. 9a). The putative Sig- 
illaria—Diaphorodendron-Synchysidendron clade 
can be placed above Hizemodendron. Alterna- 
tively, placing Chaloneria immediately below An- 
abathra allows the Sigillaria—Diaphorodendron- 
Synchysidendron clade to be situated between or 
immediately below these genera. Together, the five 
genera can form a sister clade to Hizemodendron— 
Lepidodendron—Lepidophloios, either with Chal- 
oneria~Anabathra and Sigillaria—Diaphoroden- 
dron—Synchysidendron as sister groups or as a 
nested clade ((((Diaphorodendron—Synchysiden- 
dron) Sigillaria) Anabathra) Chaloneria). 


Analysis B (all OTUs, “habit” characters C1- 
C7 omitted). The two MPTs of analysis B are 
identical to two of the three MPTs of analysis A 
(Fig. 9a, b). Five additional topologies occur at 
Lin+13 One fails to resolve the relationship between 
Chaloneria and Anabathra (dashed line on Fig. 
9c), two unite Sigillaria and Diaphorodendron- 
Synchysidendron (Fig. 9e, f), and two allow dis- 
sociation of the Oxroadia—Paurodendron clade 
(Fig. 9d, h). At Lan+2, Chaloneria can be united 
with Sigillaria, Hizemodendron and Sigillaria 
can be transposed across Diaphorodendron-Syn- 
chysidendron, and the putative Diaphoroden- 
dr on-Synchysidendron-Sigillaria clade shown in 
Figure 9e can be placed above Hizemodendron or 
below Anabathra and Chaloneria. All topologies 
found in analysis B at L,,,,. were also found in 
analysis A at Long. 


Analysis C (Chaloneria omitted, all characters 
included) . The only MPT from analysis C merely 
deletes Chaloneria from the preferred MPT of 
analysis A (Fig. 9a). Increasing the number of steps 
allow a Sigillaria—Diaphorodendron—Synchyst- 
dendron clade, either below (L, 1) or above (Lan+2) 
Hizemodendron, and disaggregation of the 
Oxroadia~Paurodendron clade (L,,,,.,). The range 
of topologies substantially increases at L,,,,.3. Once 
again, Sigillaria and Hizemodendron can be ex- 
changed across Diaphorodendron—Synchysiden- 
dron. A putative Sigillaria—Diaphorodendron— 
Synchy sidendron clade can be placed low in the 


tree, between Oxroadia—Paurodendron and An- 
abathra, or it can be appended to Anabathra to 
yield a more innovative topology. 

A strict consensus tree (Nelson, 1979, 1983) 
at L,,,.42 (analyses A, B) or L,,,,, (analysis C), has 
only four nodes (Fig. 91); only Oxroadia and Pau- 
rodendron (at the base of the tree), Diaphoro- 
dendron and Synchysidendron, and Lepidoden- 
dron and Lepidophloios (at the apex) are 
consistently conjoined and/or juxtaposed. 


Analysis D (all OTUs, vegetative characters 
only). The preferred MPT for analysis D (Figs. 
10a, 11) is pectinate (each genus is connected 
directly to the major axis) and differs substantially 
from the preferred MPT of analyses A-C (see also 
Bateman & DiMichele, 1991). Lepidodendron and 
Lepidophloios remain linked at the top of the tree, 
but two pairs of adjacent branches are transposed 
(Sigillaria and Diaphorodendron-Synchysiden- 
dron, Anabathra and Chaloneria). The Oxroad- 
ia—Paurodendron clade is split into its constituent 
genera, each of which forms an equally parsimo- 
nious sister group to the rest of the ingroup (Fig. 
10a, b). Hizemodendron is sister group to the most 
derived genera Lepidodendron and Lepidophloios 
in the preferred MPT of analysis A (Fig. 9a), but 
is derived relative only to Oxroadia and Pauro- 
dendron in the preferred MPT of analysis D. 

All of the variation among the eight additional 
generic topologies at L,,,,,, occurs below Diapho- 
rodendron—Synchysidendron in the tree, indicat- 
ing that the more apomorphic portion of the tree 
is the most robust. An unresolved trichotomy re- 
places the Oxroadia—Paurodendron clade (Fig. 
10c). The remaining topologies place Anabathra 
below Chaloneria (Fig. 10d-f) or unite Anabathra 
and Chaloneria as a separate clade (Fig. 10g-)). 
Hizemodendron is the least stable genus; it can 
occur below (Fig. 10a-e, g, i) or above (Fig. 10f, 
h, j) Anabathra and Chaloneria. In the most rad- 
ical topology, Hizemodendron is the sister group 
of Anabathra and Chaloneria, together forming 
a separate clade (Fig. 10k). 


Analysis E (all OTUs, reproductive characters 
only) . In contrast with the other analyses, it was 
not possible using PAUP 2.4 to store and thereby 
screen all 810 MPTs for analysis E (270 trees, if 
the equally most parsimonious solutions to the Lep- 
idophloios species relationships are ignored). 
Available evidence suggests that the number of 
MPTs was grossly exaggerated by repeated poly- 
chotomies and conceals a much smaller number of 
substantially different topologies. A representative 
and fully annotated MPT is shown in Figure 12a. 
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FicureE 10. Generic topologies of analysis D (vegetative characters only) at Lan (a, b) and Lanyi (C me 
caption to Table 3 for key to abbreviations.) Only those portions of trees that differ from the preferred most pars 
tree (a) are shown, and all character information is omitted, 


Unlike analyses A-D, the relationship between the 
four most primitive OTUs is unresolved, as are the 
relationships between (1) the three species of Dia- 
Phorodendron sens. str., and (2) Hizemodendron 
and Lepidodendron. The most radical innovation 
is the depiction of Lepidophloios as polyphyletic; 
the relatively primitive L. harcourtii is separated 
from L. johnsonii and L. hallii by Lepidodendron 
and Hizemodendron. Restoring Lepidophloios to 
monophyly costs one additional step (Fig. 12b). 


Methodological conclusions. For any cladistic ma- 
trix, progressive one-step increases in length rel- 
ative to the MPT result in a rapid increase in the 
number of topologies obtained (Table 4). The gra- 


dient of this increase provides an estimate pha 
crude and dependent on matrix size) of the ae 
resolution of different data matrices; more Be 
dence can be placed on a most parsimonious ‘ 
from a matrix that yields few alternative a 
optimal or near-optimal length (more rigorous; © 
tistical methods are now available for paS 
confidence limits of specific topologies; €-8- 89) 
senstein, 1985; Archie, 1989b; Sanderson, 19 Ua 
For our data, the complete data matrix (analysis 
A) provides a yardstick by which to measure y 
relative resolution of analyses based on se 
reduced permutations of the matrix present a 
Table 3 (i.e., analyses B-E). Omitting the yield: 
synapomorphic habit characters (analysis B) 
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FICURE 11. Preferred most parsimonious cladogram for analysis D (vegetative characters only). Character notation 


follows F igure 6, OTU notation follows Table 3. 


ed fewer topologies of length L,,,,—L,,,,,, and slightly 
increased the consistency index (Table 4). Omitting 
Chaloneria, the least stable OTU (analysis C), in- 
creased the consistency index by a similar amount 
to that of analysis B and generated an even more 
highly resolved set of trees that included only one 
most parsimonious topology (Fig. 9a). 
Substantially reducing the size of the data matrix 
analyzed by including only vegetative characters 
(analysis D) or only reproductive characters (anal- 
se E) also increased consistency index values rel- 
ative to those of analysis A. However, in contrast 
with analyses B and C, analyses D and E yielded 
less resolved sets of topologies (Table 4). Analysis 
D Provided acceptable results (Fig. 10), but analysis 


E yielded 270 equally most parsimonious trees, all 
containing at least one polychotomy. Thus, it is 
tempting to argue that vegetative characters are 
more phylogenetically informative (i.e., less ho- 
moplastic) than reproductive characters, but levels 
of homoplasy are very similar in the MPTs of the 
vegetative and the reproductive submatrices (con- 
sistency index values = 0.66 and 0.67 respec- 
tively; Table 4). Rather, the crucial difference ap- 
parently lies in the different sizes of the submatrices, 
which are reflected in different values for the av- 
erage number of synapomorphic character states 
per OTU: 3.2 in analysis D and 1.8 in analysis E. 
For the preferred MPT of our complete matrix, 
the number of steps per synapomorphy (1.6) and 
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FIGURE 12. a. Preferred most parsimonious cladogram for analysis E (reproductive characters only). os 
notation follows Figure 6, OTU notation follows Table 3.—b. Terminal portion of an alternative topology at baint 


the resulting consistency index (0.63) are average 
relative to those of other cladistic data matrices 
containing similar numbers of OTUs (the main 
variable influencing consistency index values: cf. 
figs. 2a and 3 of Archie, 1989a; fig. 1b of San- 
derson & Donoghue, 1989). Thus, we suggest that 
for a data matrix of average homoplasy (as here), 
there is a threshold of minimum empirical support 
(2-3 synapomorphies per OTU) below which well- 
resolved sets of fully dichotomous trees cannot be 
expected (see also Felsenstein, 1985; Guyer & 
Slowinski, 1991), 


ES 
EMPIRICAL SUPPORT FOR ALTERNATIVE HYPOTHES 
OF GENERIC RELATIONSHIPS 


The relative merits of the topologies described 
above are best assessed by examining the op 
distributions of character state transitions er 
the trees and subjectively estimating the pro j 
ities of alternative evolutionary scenarios for par 
ticular suites of characters. To this end, many T 
the near-optimally parsimonious pairings of gen®" 
obtained during the topological survey of ont 
A are summarized in a single fully annotated 


Volume 79, Number 3 
1992 


Bateman et al. 
Arborescent Lycopsid Phylogeny 


535 


(Fig. 7), to enable comparison with the preferred 
most parsimonious tree (Fig. 6). 


Oxroadia and Paurodendron. Oxroadia and Pau- 
rodendron share only two synapomorphies (Fig. 
6). Both are nonhomoplastic and describe echinate 
distal and proximal surfaces of microspore exines 
(C111, C115). Disaggregating this clade induces 
parallelism in these characters, at the cost of one 
step (Fig. 7). Given that echinate microspore exines 
could be regarded as a single character, the status 
of these two genera as a monophyletic group may 
be even less well supported than parsimony sug- 
gests. A phylogenetic study focusing specifically on 
these and other similarly primitive genera is re- 
quired to resolve this ambiguity. 


Anabathra and Chaloneria. Topologies treating 
this group as paraphyletic (Fig. 6) or as a mono- 
phyletic clade (Fig. 7) are equally parsimonious, 
differing only in the distributions of a few homo- 
plasies. When united, the genera possess only two 
synapomorphies: weakly branched or unbranched 
trunk apex (C3) and acutely angled leaf attachment 
(C66). Both characters are parallelisms: trunks with 
little or no branching are depicted as homoplastic, 
also occurring in Sigillaria and Diaphoroden- 
dron. Angle of leaf attachment, a character prone 
to ontogenetically related variation and ecopheno- 
typic modification, is represented as a parallel trait 
in Paurodendron and Hizemodendron. Thus, ev- 
idence for an Anabathra—Chaloneria clade is weak. 
Uniting Anabathra and Chaloneria affects op- 
timized character state transitions elsewhere in the 
tree. In Figure 6, unbranched trunk apex (C3) is 
perceived as reversed in Synchysidendron and in 
the Hizemodendron—Lepidodendron—Lepido- 
Phloios clade, but in Figure 7 this character state 
is depicted as a parallel acquisition in Anabathra— 
Chaloneria, Diaphorodendron, and Sigillaria. The 
absence of a ligule cavity (C60) in Anabathra is 
a plesiomorphy in Figure 6 but a reversal in Figure 
7. Missing values for C33 (secretory intracortical 
leaf sheaths) and C73 (branch gap associated with 
peduncle) in Chaloneria allow these characters to 
be optimized as apomorphic in Figure 6 (despite 
the absence of a cone in Chaloneria) but plesiomor- 
phic in Figure 7. 


Anabathra, Chaloneria, and Sigillaria. In anal- 
ysis A, these genera formed a clade at L,,.+1 (Fig- 
9k). Chaloneria and Sigillaria are united by two 
nonhomoplastic character states, fusiform rootlet 
gaps (C13) and V-shaped leaf traces (C69), and 

the homoplasy of abaxial grooves in the leaf 
(C70). The apparent synapomorphy for C69 may 


be false if the V-shaped leaf trace of Sigillaria 
reflects origination from two protoxylem strands, 
and the abaxial grooves may be developmentally 
related to leaf trace morphology. Anabathra is 
united with Chaloneria and Sigillaria on the basis 
of two homoplasies reflecting habit: the trunk pos- 
sesses an apex that shows little if any branching 
(C3) and bears lateral branches and/or cones (C6; 
a character that could not be coded for the branch- 
less and coneless Chaloneria). Both states also 
characterize Diaphorodendron (Fig. 6). 


Diaphorodendron and Synchysidendron. The 
distinction between the segregate Synchysiden- 
dron (S. dicentricum and Synchysidendron sp. 
nov.) and the three species of Diaphorodendron 
sens. str. is well supported (Fig. 6), the former by 
seven synapomorphies (four nonhomoplastic) and 
the latter by four synapomorphies (two nonhomo- 
plastic). These characters, which represent habit 
sens. lat. (C3, C6-C7, C74), leaf base retention 
(C44-C45), stele histology (C17, C20, C31), and 
periderm histology (C42-C43), reflect two sub- 
stantially different growth habits (see Evolutionary 
Patterns). Nesting Synchysidendron within Dia- 
phorodendron requires a minimum of four extra 
steps in analyses A, C, and D, and three extra 
steps when habit characters C1-C7 are omitted 
(analysis B). Placing the genera on separate ter- 
minal branches is even less parsimonious. Although 
there is little doubt that Diaphorodendron sens. 
lat. (i.e., sensu DiMichele, 1985) is monophyletic, 
we believe that the differences between the two 
monophyletic groups of species that it contains are 
sufficiently profound to warrant segregation of the 
new genus, Synchysidendron (Appendix 1C). 


Sigillaria and Diaphorodendron-Synchysiden- 
dron. Sigillaria and Diaphorodendron-Synchy- 
sidendron are empirically the best supported of all 
the clades under scrutiny; Sigillaria is supported 
by 14 characters (six nonhomoplastic), and Dia- 
phorodendron-Synchysidendron by 15 charac- 
ters (nine nonhomoplastic). Uniting Sigillaria with 
Diaphorodendron-Synchysidendron adds only 
one step and transforms two reproductive char- 
acter states from parallelisms into nonhomoplastic 
synapomorphies: short, horizontal alations on the 
megasporophyll pedicel (C82), and crassitude mi- 
crospore laesurae (C107). Balancing these gains, 
homoplasy is induced by losses in three character 
states: the presence of a lower keel on the leaf 
cushion (C55), germination of megaspores within 
the sporangium (C80), and reduction of viable 
megaspores to one per sporangium (C90, a logical 
functional correlate of C80). Although Figure 7 
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portrays these character state transitions as re- 
versals in Sigillaria, they are more likely to rep- 
resent parallel acquisitions in Diaphorodendron— 
Synchysidendron and in the more derived Hize- 
modendron—Lepidodendron—Lepidophloios clade 
(loss of these traits in Sigillaria would probably 
confer a severe competitive disadvantage). As all 
five of the above characters are considered poten- 
tially homoplastic, the possible monophyly of Sig- 
illaria—Diaphorodendron—Synchysidendron re- 
mains equivocal. 

Differential optimization alters perceptions of 
evolutionary patterns of three other characters in 
Figure 7. The first is stem apical branching (C3, 
already discussed under Chaloneria—Anabathra), 
and the remaining two describe leaf cushion mor- 
phology. In Figure 6, the plicate lower field (C57) 
evolves below Diaphorodendron—Synchysiden- 
dron and is subsequently lost in Lepidophloios, 
whereas in Figure 7 it is represented as a parallel 
acquisition in Diaphorodendron—Synchysiden- 
dron, Hizemodendron, and Lepidodendron. Sim- 
ilarly, in Figure 6, the upper keel evolves below 
Diaphorodendron-Synchysidendron and is lost in 
Hizemodendron, whereas in F igure 7 it evolves 
below the Oxroadia—Paurodendron clade (despite 
the absence of leaf cushions in these highly ple- 
siomorphic genera) and is independently lost in 
Sigillaria and Hizemodendron. For both char- 


acters, the optimizations in Figure 6 are more 
intuitive. 


Hizemodendron and Lepidodendron. As depicted 
in Figure 6, the branches immediately subtending 
Hizemodendron and Lepidodendron share only 
one character state transition: the loss of secretory 
intracortical leaf-trace sheaths (C33). The other 
four characters that support the Hizemodendron— 
Lepidodendron clade in Figure 7 are all parallel- 
isms. Three describe microspore equatorial (C103) 
and contact face (C109, C112) ornamentation and 
are homoplastic among Lepidophloios species. 
Their frequencies and distributions differ between 
topologies, due to inconsistent optimization of miss- 
ing values. These characters are more appropri- 
ately treated as reversals within Lepidophloios 
(Fig. 6) than as parallelisms in Lepidophloios and 
Hizemodendron—Lepidodendron (Fig. 7). In Fig- 
ure 6, the fourth character state, plication of the 
lower field of the leaf cushion (C57), originates 
below Diaphorodendron—Synchysidendron and is 
reversed in Lepidophloios. In Figure 7, this char- 
acter state originates twice, in Diaphorodendron— 
Synchysidendron and in Hizemodendron (a less 
probable scenario). Missing values allow demotion 


of infrafoliar parichnos from a nonhomoplastic syn- 
apomorphy of Lepidodendron and all Lepido- 
phloios species (Fig. 6) to a parallelism of Lepi- 
dodendron and Lepidophloios hallii only (Fig. 7). 
Similarly, a missing value in Hizemodendron for 
multizoned periderm (C39) allows the genus to be 
plesiomorphic for this character in Figure 6 but 
apomorphic in Figure 7. 

However, the most deleterious consequence of 
uniting Hizemodendron and Lepidodendron is the 
generation between Lepidodendron and Lepido- 
phloios of five vegetative parallelisms: discernable 
protoxylem ridges on the stele (C26), resinous peri- 
derm (C48), tangentially elongate leaf cushions on 
twigs (C53), and a leaf with a dorsiventrally flat- 
tened vascular strand (C68) and lateral abaxial 
grooves (C70). The overall cost to parsimony (five 
steps) appears sufficient to reject the hypothesis of 
monophyly. Nevertheless, treating these hypoth- 
esized character state transitions as reversals in 
Hizemodendron suggests a heterochronic evolu- 
tionary mechanism that could allow monophyly (see 
Evolutionary Patterns). 


Lepidodendron sens. lat. Prior to the studies of 
DiMichele (1985), Bateman & DiMichele (1991), 
and DiMichele & Bateman (1992), ‘Lepido 
dendron’ sens. lat. encompassed four of 
anatomically preserved genera analyzed by u5 
Lepidodendron sens. str., Hizemodendron, Dia- 
phorodendron, and Synchysidendron. Forcing 
these four genera into a single clade repor 
the traditional concept of Lepidodendron (not t- 
lustrated) cost nine additional steps and reveal 
only one synapomorphy uniting the clade: a plicate 
lower field to the leaf cushion, which is reve 6) 
in Lepidophloios in the preferred MPT (Fig. 0- 
In contrast, five nonhomoplastic synapomorphis 
in Figure 6 are rendered homoplastic: short, pit 
megasporophyll alations (C83) and bilaterally fa i 
tened (C85), distally dehiscent (C87) megasport 
gia are lost in Diaphorodendron—Synchys! = 
dron, and zoned periderm (C39) and — d 
parichnos (C62) become parallelisms in Lep i 
dendron and some Lepidophloios species wie 
result of ambiguous missing values in the noe 
Moreover, homoplasy is increased in other $ 
acters already depicted as homoplastic m bag 
6. We conclude that ‘Lepidodendron’ sens. 1 
clearly a paraphyletic group. 


MULTIVARIATE ANALYSIS 


ical 
In order to examine patterns of morpholog! 
variation free from the rigid constraints ee 
by cladistic nested hierarchies, we subject 
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cladistic data matrix (Table 3) to multivariate anal- 
ysis. A value of zero or unity was substituted a 
priori for each missing value according to optimized 
distributions of character state transitions in the 
preferred most parsimonious cladogram (Table 3 
subscripts, Fig. 6). The resulting uniformly binary 
matrix allowed generation of a symmetrical matrix 
comparing OTUs without a priori standardization, 
simply using the number of character state conflicts 
(ie., O vs. 1) as a direct measure of dissimilarity 
between pairs of OTUs. The dissimilarity values 
were used to construct an unrooted minimum span- 
ning tree (Gower & Ross, 1969); links in the tree 
represent specific sets of character state transitions, 
thus contradicting frequent assertions that phenetic 
trees inevitably lack such information. Also, prin- 
cipal coordinates (Gower, 1966) were calculated 
from the data matrix via Manhattan distances, 
using unpublished software written by J. Alroy. 
Holapomorphies (which are invariant) were ex- 
cluded (as in cladistic analysis), but autapomorphies 
contributed to both the unrooted tree and the or- 
dination. 

Links between genera on the minimum spanning 
tree (Fig. 13a) represent at least 10 character 
conflicts, those within genera represent no more 
than five. The 21 conflicts between Hizemoden- 
dron and Lepidodendron, and 14 conflicts be- 
tween Synchysidendron and Diaphorodendron, 
emphasize the need to segregate these new genera 
(Appendix 1C). 

The minimum spanning tree resembles the pre- 
ferred most parsimonious cladogram (Fig. 6) in 
depicting a progression from Paurodendron and 

roadia through Anabathra, Hizemodendron, 
and Lepidodendron to Lepidophloios (though L. 
johnsonii is shown as ancestral to the two remain- 
'ng species). However, Anabathra is also depicted 
as the ancestor of a second lineage, consisting of 
Chaloneria, Sigillaria, Diaphorodendron, and 
Synchysidendron, that is not represented among 
the range of cladograms shown in Figure 9. This 
Second lineage is held together by the weakest links 
ont the tree (Chaloneria—Sigillaria = 29 conflicts, 
Sigillaria~Diaphorodendron = 33 conflicts) and 
consequently can be dissociated at the cost of very 
few additional steps (Fig. 13b), demonstrating that 
these genera are the most problematic in both the 
cladistically and phenetically generated phyloge- 
mes, Attaching Sigillaria directly to Anabathra, 
rather than via Chaloneria, creates an intuitively 
more credible evolutionary hypothesis at the ex- 
Pense of two steps. Diaphorodendron and Syn- 
chysidendron can be attached to Hizemodendron 
at the cost of only one step, but this results in the 
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FicureE 13. a. Minimum spanning tree (unrooted phe- 
netic tree) as based on cladistic data matrix (Table 3, q.v. 
for key to OTUs). Links represent numbers of character 
state conflicts. —b. Alternative links that are almost max- 
imally parsimonious; additional steps relative to a are 


marked. 


putatively more derived Synchysidendron giving 
rise to the more primitive Diaphorodendron. Pau- 
rodendron can replace Oxroadia as the closer 
relative of Anabathra at the cost of one step. 

In summary, the phenetic trees serve primarily 
to emphasize the potentially pivotal role of Ana- 
bathra as the most primitive arboreous lycopsid 
analyzed. 

The first three principal coordinates (Fig. 14) 
account for an unusually large proportion (91%) 
of the total variance. The first coordinate separates 
Diaphorodendron—Synchysidendron from the 
other genera, the second coordinate separates Lep- 
idodendron—Lepidophloios from the bisporan- 
giate-coned group of Paurodendron—Oxroadia— 
Anabathra—Chaloneria (with Hizemodendron in- 
termediate), and the appreciably weaker third co- 
ordinate separates Sigillaria (and, to a lesser ex- 
tent, Chaloneria) from the remainder. The resulting 
tetrahedral arrangement of four clusters (excluding 
Hizemodendron), separated by broad morpholog- 
ical discontinuities, underlines the distinctiveness 
of the three groups of monosporangiate-coned trees 
(Sigillaria, Diaphorodendron—Synchysiden- 
dron, Lepidodendron-Lepidophloios) and the 
consequent difficulties of resolving their phyloge- 
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FIGURE 14. Principal coordinates ordination of the 
OTUs, based on Table 3 (q.v. for OTUs). Only genera 
are labeled. Two Diaphorodendron species are indistin- 
guishable on the first three coordinates. 


netic relationships relative to each other and to 
Anabathra, their most likely sister group. 


EVOLUTIONARY PATTERNS AND PROCESSES 


Overall trends. Much of the variation among ar- 
borescent lycopsids can be resolved into a vege- 
tative trend, reflecting in particular morphological 
and anatomical expression of different growth ar- 
chitectures, and a reproductive trend, representing 
increasingly sophisticated reproductive strategies. 
The phylogenetic analyses show that the two trends 
are not entirely concordant; the preferred MPTs 
using vegetative characters only (analysis D; Fig. 
11) and reproductive characters only (analysis E; 
Fig. 12) have substantially different topologies. The 
following discussion of these trends emphasizes 
characters for which we developed strong (often a 
priori) hypotheses of high burden and pays partic- 
ular attention to the relative temporal order of 
appearance of apomorphic states of different char- 
acters (e.g., Donoghue, 1989), 

Many (possibly all) of the MPTs of analysis E 
(reproductive characters only) distinguish the four 
bisporangiate-coned genera (Paurodendron, 
Oxroadia, Anabathra, Chaloneria) from the five 
relatively derived genera that possess a suite of 
characters reflecting the developmental partition- 
ing of mega- and microsporangia into monosporan- 
giate cones (Fig. 12). In contrast, the preferred 
MPT of analysis D (vegetative characters only) 
distinguishes four primitive pseudoherbs/ shrubs 
(Paurodendron, Oxroadia, Hizemodendron, 
Chaloneria) from five derived arboreous genera 
(the derived clade is supported by the tree habit 
only; Fig. 11) and depicts Sigillaria as derived 


relative to Diaphorodendron—Synchysidendron, 
The genera whose positions differ most between 
the reproductive and vegetative cladograms are 
the bisporangiate-coned tree Anabathra, which is 
reproductively plesiomorphic and vegetatively apo- 
morphic, and the monosporangiate-coned pseu- 
doherb Hizemodendron, which is reproductively 
apomorphic and vegetatively plesiomorphic (see 
also Bateman & DiMichele, 1991). 

Thus, monosporangiate cones (Fig. 12) and the 
tree habit (Fig. 11) cannot both be nonhomoplastic, 
though a full analysis using vegetative and repro- 
ductive characters together could have yielded a 
compromise solution involving homoplasy in both 
suites of characters. In fact, although its topology 
differs in detail from those of analyses D and E, 
the preferred MPT for analysis A (all characters 
and OTUs; Fig. 6) more closely resembles the 
exclusively reproductive cladogram (Fig. 12) than 
the exclusively vegetative cladogram (Fig. 11); in 
particular, the monosporangiate-coned clade is re- 
tained at the expense of depicting the tree habit 
as homoplastic. 

The two subsections that follow discuss in greater 
detail the reproductive and vegetative trends, fo- 
cusing on the functional morphology and adaptive 
(or nonadaptive) significance of specific character 
states, before returning to the evolutionary impli- 
cations (and limitations) of the preferred whole- 
organism cladogram. 


Reproductive morphology. Reproductive char- 
acters proved to be of little value in elucidating 
phylogenetic relationships within the plesiomorphic 
group of bisporangiate-coned OTUs, which differ 
primarily in autapomorphic spore character states 
(Fig. 12). Careful revision of bisporangiate cone 
is desirable, to search for potential synapomorphies 
less inclusive than the entire group. It is particu- 
larly important to understand the ontogeny ® 
reproductive biology of the bisporangiate cones 
order to determine how they could have given ms 
to monosporangiate-coned descendants. All ei 
rangiate cones have apically concentrated ai 
sporangia and basally concentrated megasporang™ 
irrespective of presumed geotropic orientation, oe 
gesting a shared developmental control of sp? 
genesis. k 
In contrast, the functional morphology of ce 
productive characters within the monoeparne a 
coned portion of the lepidodendralean are: 
prompted much discussion (e.g., Thomas, l Phil- 
Phillips, 1979; DiMichele & Phillips, 1985; ‘i 
lips & DiMichele, 1992; see Appendix 1D for 
taxonomic implications). In Figures 6 and 12, 
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appearance of monosporangiate cones (immediate- 
ly below Sigillaria) is accompanied by lateral ex- 
pansion of the sporophyll pedicel to form alations 
and by the functionally important transition in basic 
dispersal unit from isolated megaspores to the 
megasporophyll-megasporangium complex. All 
three character states persist without reversal 
throughout the derived clade, indicating strong 
functional linkage that may be evolutionarily tied 
to elaboration of the leaf bases. 

Current evidence suggests that the leaf lamina 
and leaf base are derived from the same primor- 
dium; the leaf cushion, an elaborated leaf base, is 
also fundamentally foliar. Furthermore, we believe 
that the sporophyll lamina is homologous with the 
leaf lamina, and that the sporangium-bearing ped- 
icel is homologous with the leaf base, including an 
elaborated cushion if present. The strongest evi- 
dence supporting these homologies is provided by 
the ligule (e.g., Phillips, 1979; Bateman, 1988), 
which occurs adaxially on (Paurodendron) or with- 
in (all other OTUs) the bases of leaves (F ig. 4b) 
but on sporophylls occurs close to the distal end 
of the pedicel, between the sporangium and the 
more-or-less perpendicular junction of the pedicel 
and lamina (Fig. 4d). Regarding the attachment of 
the ligule as a homologous point implies that the 
pedicel is indeed homologous with the leaf base. 
Both the leaf base (and thereby cushion) and lamina 
originate from the same primordium, as do the 
sporophyll pedicel and lamina. Leaf and sporophyll 
both bear a ligule and both are fundamentally ap- 
pendicular in origin. Moreover, there is a strong 
Positive correlation between the complexity of the 
pedicel-sporangium unit and that of the leaf base; 
definable leaf cushions appear at the same node of 
the cladogram as monosporangiate cones, and both 
Structures progressively increase in complexity 
through the remainder of the clade, culminating in 
the large, elaborate leaf cushions and equally large, 
seedlike megasporophyll of Lepidophloios (e.g., 
Reed, 1941; Phillips, 1979). Thus, the evolution 
of the leaf cushion may have been developmentally 
linked to that of the sporophyll (T. L. Phillips, pers. 
comm. 1989). It is not clear whether elaboration 
of the leaf prompted modification of the sporophyll, 
Door whether increase in size of the appendicular 
Primordial meristems allowed simultaneous expan- 
“on and elaboration of both leaf bases and spo- 
rophyll pedicels. 

hatever its driving mechanism, the transition 
tom bisporangiate to monosporangiate cones (im- 
mediately below Sigillaria on Figs. 6 and 12) 
represents a crucial release from developmental 
Constraints (cf, Endress, 1987, on angiosperms). 


In particular, it allowed spatially independent de- 
velopment of the mega- and microsporophylls, 
thereby permitting modifications of the megaspo- 
rophyll-megasporangium units that could have im- 
paired the function of microsporophyll units if sim- 
ilarly modified (a likely consequence in bisporangiate 
cones, where mega- and microsporophylls form a 
developmental continuum). Thus, free megaspores 
were superseded as the basic dispersal unit by the 
megasporophyll-megasporangium complex. 

The remaining reproductive modifications that 
delimit increasingly exclusive portions of the mono- 
sporangiate-coned clade can be envisioned as a 
progressive evolutionary trend toward K-selection 
(sensu Pianka, 1970). Megaspores decrease in 
number and increase in size, and the pedicellate 
tissues surrounding the megasporangium become 
adapted for increasingly specialized modes of mi- 
crospore/microgametophyte capture and diaspore 
dispersal (e.g., Phillips, 1979). However, argu- 
ments that the most derived product of this evo- 
lutionary trend, Lepidophloios—Lepidocarpon, 
possesses true seeds (Zhang et al., 1986) are phy- 
logenetically unhelpful; the megasporophylls are 
clearly analogs rather than homologs of gymno- 
sperm ovules. 

Reduction in megaspore number to one per spo- 
rangium, and concomitant germination of mega- 
spores within the sporangium, distinguish the re- 
mainder of the monosporangiate-coned clade from 
Sigillaria (Figs. 6, 12). The dichotomy immedi- 
ately above Sigillaria results in two clades well 
supported by reproductive characters: Diaphoro- 
dendron-Synchysidendron (proximally dehiscent, 
dorsiventrally flattened, heterocellular megaspo- 
rangium containing gulate megaspores; granulate- 
foveolate Granasporites microspores), and Hize- 
modendron—Lepidodendron—Lepidophloios (dis- 
tally dehiscent, cylindrical sporangium subtended 
by suberect alations; cingulate Lycospora micro- 
spores). Enveloping, integumentlike alations delimit 
Lepidophloios (Fig. 6; in the analysis of repro- 
ductive characters only (Fig. 12a), this single char- 
acter is insufficient to override differences in mi- 
crospore ornamentation between Lepidophloios 
species, resulting in depiction of the genus as poly- 
phyletic; inducing monophyly in Lepidophloios 
costs one extra step (Fig. 12b)). Interestingly, the 
vegetatively well-differentiated species of Diapho- 
rodendron and Synchysidendron (Fig. 11) are 
effectively uniform in reproductive characters (the 
medullated stele attributed to Diaphorodendron 
(C74) strictly applies to the cone-bearing lateral 
branches rather than the cones per se; Appendix 


1D). 
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Distinguishing between monoecious and dioe- 
cious strategies is especially valuable in interpreting 
the phylogeny and the functional morphology of 
extant plants (e.g., Bawa, 1980; Givnish, 1980, 
1982; Donoghue, 1989). Unfortunately, the al- 
most inevitable disarticulation of cones from veg- 
etative axes undermines attempts to identify the 
reproductive strategies adopted by members of the 
monosporangiate-coned lycopsid clade. Were in- 
dividual plants monoecious or dioecious? If mon- 
oecious (as seems more likely), did megasporangiate 
and microsporangiate cones mature synchronously 
or sequentially? Or were arborescent lycopsids ca- 
pable of even more complex strategies, such as 
gynodioecy? 

Overall, spore morphology proved less phylo- 
genetically informative than cone morphology. 
Many of the spore character states are species- 
level autapomorphies, and the remainder exhibit 
significantly greater homoplasy than cone char- 
acters (in Fig. 6, the proportion of holapomorphies 
plus autapomorphies and consistency index for cone 
and sporophyll characters are 26% and 0.78 re- 
spectively, contrasting with values of 42% and 
0.64 for mega- and microspore characters). In 
retrospect, some cases of mistaken homology are 
clearly evident among spore characters. For ex- 
ample, the distal spines (C100) of Oxroadia and 
Sigillaria sp. nov. megaspores differ in detail; those 
of the former are long and buttressed, those of the 
latter are short and almost papillate. Also, polar- 
ization of spore characters was especially problem- 
atic, as character states (especially those repre- 
senting ornamentation) tend to replace each other 
in entirety rather than accumulating as sequential 
elaborations of form (i.e., they are displacive rather 
than additive). Consequently, patterns of increasing 
complexity cannot be expected. Indeed, many spore 
character states (including the more elaborate forms 
of ornamentation) are confined to the primitive, 
bisporangiate OTUs. Reduction to a single func- 
tional megaspore resulted in the loss of all types 
of dispersed ornamentation (Table 3), suggesting 
that they were redundant once megaspores had 
ceased to be the basic units of dispersal. 

Relatively little attention has been paid to the 
functional morphology of lycopsid spores. Promi- 
nent equatorial (Paurodendron, Chaloneria) and 
laesural (Paurodendron, Oxroadia, Diaphoroden- 
dron, Synchysidendron) megaspore ornamenta- 
tion has been invoked as an aid to flotation and 
thereby dispersal. Phillips (1979) argued that open- 
ings in the massa of Diaphorodendron and Syn- 
chysidendron megaspores trapped microspores to 
facilitate fertilization, and a similar function was 
attributed to the anastomosing fimbriate laesural 
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ornamentation of Oxroadia megaspores by Bate- 
man (1988). The pseudosaccus of Chaloneria mi- 
crospores and cingulum of Hizemodendron, Lep- 
idodendron, and Lepidophloios microspores 
probably acted as buoyancy aids. 


Vegetative morphology. Vegetative characters 
uniting the arborescent lycopsids are those asso- 
ciated with the production of rhizomorphs and wood, 
together with ligules and exarch xylem maturation 
(inferred centrifugal maturation in stigmarian rhi- 
zomorphs (Frankenberg & Eggert, 1969) is suspect 
(Phillips & DiMichele, 1992)). Beginning at the 
root of the preferred MPT and passing along its 
major axis (Fig. 6), the node above Paurodendron- 
Oxroadia is characterized by the appearance of 
the tree habit (habit is discussed more fully in 
subsequent subsections) and by modifications to 
vascular tissue that probably reflect greatly in- 
creased body size: medullation of the stele, contin- 
uous protoxylem sheath, loss of protoxylem ridges 
(together causing superficial leaf trace emission: 
Fig. 5b), and the advent of foliar parichnos in the 
leaf bases, now much more distant from the axial 
vasculature following acquisition of the arboreous 
habit. Beyond Anabathra, the ligule pit is ubiq- 
uitous, though it also occurs in the ostensibly prim- 
itive genus Oxroadia. 

Beyond Chaloneria, the evolution of discrete 
leaf cushions provided a consistent basal limit to 
leaf atrophy. Many of the characters that support 
nodes higher in the cladogram represent elabora- 
tions of leaf cushion morphology, notably in cushion 
complexity below Diaphorodendron-Synchyst 
dendron, and in overall shape below Lepidoden- 
dron—Lepidophloios. In contrast, cushion and leaf- 
trace simplification is evident in Hizemodendron. 
The analysis of vegetative characters only (Fig. 
11) places Hizemodendron much lower in the tree 
eliminating the many character losses shown ” 
Figure 6 but depicting the leaf cushion per 5° us 
iterative in (a) Hizemodendron and (b) all of the 
arboreous genera, shown clustered above Cha 
neria in Figure 11 (also, arbitrary optimization Ks 
missing values representing the absence of cod 
structures led to the nonsensical apparent evoh 
of specific leaf cushion features below Paur ane 
dron, prior to the evolution of the cushions them 
selves). Further modifications of cushion eee 
together with the appearance of infrafoliar eo 
nos and a return to longitudinal protoxylem pel 
characterize the highly derived Lepidodendr ci 
Lepidophloios clade (Fig. 6). f 

Stelar characters, together with peridermal e 
tures, play important roles in the more deri Ps: 
portion of the arborescent lycopsid clade, partic 
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ularly in delimiting the three main cushion-bearing 
clades: Sigillaria, Diaphorodendron—Synchysi- 
dendron, and Hizemodendron—Lepidodendron— 
Lepidophloios. The solid protosteles of Oxroadia 
and Paurodendron give way to medullation by 
parenchymatous vitalization in the remaining gen- 
era. In most of these, the central portion of the 
stele consists of parenchyma cells that have lengths 
and diameters similar to those of the innermost 
metaxylem tracheids; these appear to be procam- 
bial cells that remained unlignified and thus met- 
abolically active throughout the life of the plant. 
In the largest axes of Lepidodendron and Lepi- 
dophloios (less L. harcourtii), the stelar core con- 
sists of filamentous cells that probably proliferated 
into a central void. In contrast with the other 
lycopsid genera but in parallel with ferns, the Dia- 
phorodendron—Synchysidendron clade evolved a 
true siphonostele. Pith parenchyma cells are much 
shorter and narrower than the adjacent metaxylem 
tracheids, suggesting different developmental ori- 
gins for these tissues. Diaphorodendron has a mixed 
pith of parenchyma and tracheids, with parenchy- 
ma increasing in relative abundance toward the 
centers of larger axes. Synchysidendron has a pith 
region sharply delineated from the tracheary cells, 
and wood with heterogeneous rays and deep pa- 
renchymatous invaginations. 

Most characters of the stelar margin constitute 
Senus-level autapomorphies, notably the distinctly 
different modes of leaf trace emission observed in 
Chaloneria, Sigillaria, and Lepidophloios (Fig. 
5). Of greater interest is the apparent switch from 
distinct longitudinal protoxylem ridges (yielding 
_Coronate” cross sections) to no discernible ridges 
immediately above the primitive Paurodendron— 
Oxroadia clade, followed by a return to similar 
(but not identical) coronate morphology in the most 
derived Lepidodendron—Lepidophloios clade; the 
Protoxylem ridges are further modified in Lepi- 
dophloios, where they anastomose (Fig. 6). Deri- 
vation of Hizemodendron from Lepidodendron 
(contra Fig. 12; see Heterochrony) would imply 
loss of protoxylem ridges. We suggest that the loss 
of ridges does not reflect complete absence of dis- 
crete protoxylem strands, even though the strands 
are no longer discernible. 

Cortical characters of the arborescent lycopsids 
are surprisingly conservative compared with the 
other axial tissues. The persistent inner cortex may 

ave provided a barrier of live cells along the outer 
margin of the phloem, protecting this delicate tissue 
fom exposure to the central void created by the 
Presumed in vivo disintegration of the thin-walled 
Parenchyma of the middle cortex. The medullated 
Steles of Diaphorodendron and Synchysidendron 
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emit leaf traces that are sheathed by parenchyma 
when passing through the outer cortex, a character 
state also found in Sigillaria sp. nov. Leaf traces 
are secretory in most of the arboreous species 
analyzed. 

Periderm is arguably the most unusual and de- 
velopmentally intriguing vegetative tissue in ar- 
borescent lycopsids. As with stelar morphology, the 
Diaphorodendron-Synchysidendron clade is dis- 
tinct from the other cushion-bearing arboreous gen- 
era. Bifaciality in the former group is evident in 
the clear histological distinction between the thin 
phellem and much thicker phelloderm. The rela- 
tively homogeneous periderm of the latter group 
may conceal cryptic cambial bifaciality (for ex- 
ample, this may be manifested in the peridermal 
trizonation of Lepidodendron—Lepidophloios), es- 
pecially if the phellem is very weakly developed or 
the phellem and phelloderm are histologically iden- 
tical. Details of periderm histology tend to be ho- 
moplastic or species-level autapomorphies, and 
therefore of limited phylogenetic value. In the anal- 
ysis of vegetative characters only (Fig. 11), pro- 
motion of Sigillaria to sister group of the Lepi- 
dodendron-Lepidophloios clade united the three 
genera that possess resinous periderm. However, 
the other two characters supporting this node 
(cushions on twigs wider than long, leaves with 
lateral abaxial grooves) are almost certainly mis- 
coded as homologs shared by Sigillaria and Lep- 
idodendron—Lepidophloios. 

Periderm, the main support tissue of the arbo- 
reous lycopsids, reached thicknesses of at least 20 
cm in some species (DiMichele, 1979a, b, 1981; 
Phillips & DiMichele, 1981). This considerable 
increase in trunk girth over that of the primary 
tissues is difficult to reconcile with the persistence 
of primary leaf cushions, which probably remained 
photosynthetic after leaf loss; they are covered in 
stomata (Thomas, 1970b, 1977; DiMichele, 1979a, 
b), and both leaf traces and parichnos connections 
with leaf cushions were maintained through the 
periderm (Delevoryas, 1957; DiMichele, 1980). 
Several specialized mechanisms for accommodat- 
ing girth increase evolved: tangential interarea ex- 
pansion in Sigillaria and Synchysidendron, in- 
terarea fissuring in Diaphorodendron, and 
subcushion cellular expansion in Lepidodendron 
(DiMichele, 1981, 1983). In arboreous genera with 
well-developed crowns (Synchysidendron, Lepi- 
dodendron, Lepidophloios), periderm thickness 
diminished through the branching systems and the 
cortex was probably a major support tissue. 

Given the determinate growth of arborescent 
lycopsids (Andrews & Murdy, 1958; Eggert, 1961), 
most of the periderm probably formed and differ- 
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entiated near the stem apex, during (though prob- 
ably continuing after) differentiation of the primary 
cortex, leaf cushions, and traces. This coordinated 
development is indicated by the persistence of the 
vascular linkages between the stele and leaf cush- 
ions through the periderm, and the occurrence of 
arboreous genera lacking specialized cushion-re- 
tention mechanisms (Anabathra, Lepidophloios). 
This mode of growth implies a stem apex analogous 
to the primary thickening meristem of some mono- 
cotyledonous angiosperms, an assertion previously 
made for stigmarian apices (Rothwell & Pryor, 
1990, 1991). Unfortunately, there have been few 
discoveries of anatomically preserved stigmarian 
apices (Rothwell, 1984; Rothwell & Pryor, 1991) 
and none of stem apices, despite detailed and pro- 
longed studies of coal balls. Their rarity may be 
at least partly explained by preferential decay of 
the apex, as observed in Oxroadia branches (Bate- 
man, 1988). Rapid trunk elongation is suspected 
in the arboreous species (Phillips & DiMichele, 
1992). 

The most significant transitions in leaf base char- 
acters are the evolution of leaf cushions immedi- 
ately below Sigillaria and their elaboration im- 
mediately above, involving the advent of upper and 
lower keels and lower field plications. The upper 
keel was subsequently lost in Hizemodendron, as 
were the lower field plications in Lepidophloios. 
Moreover, iteration is evident in several cushion 
characters: upper field plication in Diaphoroden- 
dron—Synchysidendron and Hizemodendron, the 
lateral line in Diaphorodendron—Synchysiden- 
dron and Lepidodendron, and the undoubtedly 
nonhomologous broader-than-long leaf cushions on 
the twigs of Sigillaria and Lepidophloios. Such 
convergences extend beyond the coded characters; 
for example, the similarity between Synchysiden- 
dron and Lepidodendron leaf cushions is reflected 
even in patterns of shape change from large to 
small branches. Also, cushions of both genera are 
protuberant and therefore have deep ligule pits. 
Only the infrafoliar parichnos of Lepidodendron 
allow distinction of its axial surface from that of 
Synchysidendron. Thus, many of the characters 
detailing leaf cushion morphology that have figured 
so prominently in previous classifications of Paleo- 
zoic lycopsids are shown to be homoplastic. 

We regard the leaf cushions of derived genera 
as elaborated leaf bases of their more primitive 
antecedents such as Anabathra, which lacks a 
clear lamina-cushion distinction. Several features 
of the leaf base, notably the position of the ligule, 
suggest development from a leaf primordium. If 
so, the so-called “tabscission”’ of leaves, which leaves 


a scar immediately external to the ligule pit ap- 
erture, occurred within the leaf rather than at the 
leaf-axis junction. Moreover, there is no evidence 
of a discrete abscission layer. Thus, we suspect 
that leaf laminae merely withered and sheared off 
at the physically weakest point, where the leaf 
constricts and is perforated by the ligule pit and, 
in the more derived genera, by the foliar parichnos. 
Acute leaf posture is strongly homoplastic, char- 
acterizing unrelated genera with relatively short, 
broad laminae (Paurodendron, Anabathra, Chal- 
oneria, Hizemodendron). The derived arboreous 
genera (above Chaloneria) all possess leaves with 
sclerenchymatous sheaths, which presumably sup- 
ported the long, narrow laminae. Interestingly, the 
sheaths were lost in the short-leaved pseudoherb 
Hizemodendron (Fig. 6). Expansion and invagi- 
nation of traces may also have aided structural 
support; dorsiventrally flattened traces character- 
ize the Lepidodendron—Lepidophloios clade and 
Diaphorodendron scleroticum, V-shaped traces 
occur in Chaloneria and Sigillaria. Such vascular 
elaboration is invariably manifested externally as 
lateral abaxial grooves. The V-shaped trace of Sig- 
illaria is an inevitable consequence of the vascu- 
larization of each leaf by two adjacent protoxylem 
strands of the parent axis and results in a median 
abaxial groove in addition to the lateral grooves: 
Postmortem accentuation of the grooves may have 
occurred (cf. Rex, 1986). 
Arborescent lycopsids as integrated homeostatic | 
organisms. Analogizing arborescent lycopsids with 
woody seed plants or tree ferns is unsatisfactory 
when considered in detail; in particular, perider™ 
rather than wood constitutes the main physica 
support system, and the primary function of 
branches is reproduction rather than light capture | 
(DiMichele & Phillips, 1985; Phillips & DiMichele, 
1992). Perhaps the most profound character dis- 
tinguishing the arborescent lycopsids is their de- | 
terminate growth (Walton, 1935; Andrews & Mur- 
dy, 1958; Eggert, 1961), which allows analogy 
with the ontogeny of vertebrates and legitimate U* | 
of the terms body and bauplan (Bateman 
DiMichele, 1991). : r 
The body of each arborescent lycopsid is a | 
posed of at least two of four major detenn 
structural units (modules): rhizomorph, stem 3 
fined as the length of the axis from the je ie 
root-shoot divergence to the first isotomy © 
apical meristem), crown branches (resulting ink 
isotomy of the apical meristem), and lateral brat 
es/cauline peduncles (resulting from strong p 
omy of the apical meristem) (DiMichele & | 
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man, 1989; Bateman & DiMichele, 1991). 
Rhizomorph and stem are ubiquitous modules, 
though both occur in a severely reduced form in 
some OTUs. Lateral branches/cauline peduncles 
and crown branches are each rare or absent from 
some OTUs, thus defining three basic bauplans: (1) 
neither type of branch present, (2) crown branching 
frequent, lateral branching generates peduncles 
only, and (3) crown branching infrequent, lateral 
branching dominant. Variations in the sizes and 
secondary tissue contents of all four modules, and 
in the frequency of dichotomy of lateral and crown 
branches, generate a range of growth habits (e.g., 
Hallé & Oldeman, 1970; Hallé et al., 1978; White, 
1979; Tomlinson, 1982, 1983). We suspect that 
growth of the arborescent lycopsids in general, and 
arboreous taxa in particular, was largely deter- 
ministic (genetically induced), offering much less 
potential than angiosperms for opportunistic mod- 
ification of growth architecture by environmental 
influences or chance factors (e.g., Tomlinson, 
1982). Consequently, the conceptual architectural 
model was unusually faithfully reproduced in the 
actual habit (Bateman & DiMichele, 1991). 

The three basic bauplans outlined above broadly 
correspond to three of the growth models reviewed 
by Hallé et al. (1978); two can be subdivided using 
growth habit. 

Bauplan (1) corresponds to Corner’s architec- 
tural model. In our analysis, this exclusively mono- 
axial growth is confined to Chaloneria, though 
several penecontemporaneous lycopsids apparently 
Possessed the same habit; these include Spencerites 
(Leisman & Stidd, 1967), Sporangiostrobus 
(Wagner & Spinner, 1976; R. H. Wagner, 1989), 
Porostrobus (Leary & Mickle, 1989), and the as 
yet unnamed, almost fully articulated compression 
from the Upper Devonian Cleveland Shale of Ohio 
(Chitaley, 1982, 1988; S. Chitaley & K. B. Pigg, 
in prep.). Chaloneria is the only OTU in our anal- 
ysis that consists only of the two ubiquitous mod- 
ules, a rhizomorph and an unbranched stem (Pigg 
& Rothwell, 1983a). We have classified its re- 
peated zones of cauline sporophylls as lateral rather 
than terminal fructifications, as their production 
did not necessarily result in cessation of stem growth. 

© wood cylinder is narrow and the stem erect 
but much shorter than those of the truly arboreous 
OTUs (cf. Fig. 2b with Fig. 2a-g). The low wood 
fentent and lack of branches and cones in Chal- 
oneria demonstrate highly economical construc- 
tion, implying rapid growth. Distributional evidence 
Suggests that Chaloneria was an ecological dom- 


mant in marshlike associations (DiMichele et al., 
1979), 


Bauplan (2), which corresponds to Schoute’s 
architectural model, is the most widespread among 
the OTUs. Stems divide acrotonously (distally) by 
equal division of the apical meristem, yielding mod- 
ular (determinate), orthotropic (three-dimensional) 
branches. The bauplan encompasses two distinct 
subgroups categorized by the arboreous and pseu- 
doherbaceous habits respectively. 

The first subgroup includes the classic arboreous 
genera Lepidodendron and Lepidophloios, to- 
gether with Synchysidendron (Fig. 1). Throughout 
much of their life history, these trees consist of a 
rhizomorph and telegraph polelike stem capped by 
a massive primary body, undergoing frequent di- 
chotomous branching to form a determinate crown 
only during the final phase of growth and subse- 
quent monocarpic reproduction (DiMichele & Phil- 
lips, 1985). The cones were borne on stout pe- 
duncles that we regard as the homologs of more 
elaborate lateral branches found in bauplan (3); 
although evolutionarily significant, the peduncles 
do not define the architecture of bauplan (2) trees. 
These trees were cheaply constructed. Secondary 
thickening ceased well before termination of growth; 
the resulting poor development of wood in both the 
crowns and the trunks meant that they relied pri- 
marily on periderm for structural support (Di- 
Michele, 1979a, b, 1983; Wnuk, 1985). This 
probably allowed channeling of more photosynthate 
into reproduction. Rapid generation times and an 
opportunistic life strategy were postulated for the 
subgroup by DiMichele & Phillips (1985). 

The second subgroup, consisting of Pauroden- 
dron, Oxroadia, and Hizemodendron, is delimited 
by the pseudoherbaceous growth habit (Bateman, 
1988, 1989, 1992; Bateman & DiMichele, 1991). 
These genera possess the same modules as Lepi- 
dodendron and Lepidophloios, but differ in the 
relative sizes and shapes of the modules (Fig. 2a, 
c, d). Also, overall body size of mature individuals 
is one to two orders of magnitude less than those 
of comparable arboreous species (cf. Fig. 1 with 
Fig. 2). In particular, the ubiquitous modules of 
rhizomorph and stem are greatly reduced (most 
drastically in Paurodendron) relative to the crown, 
which develops much earlier in the life histories of 
these genera. Their minimal stems result in a re- 
cumbent growth habit (Baxter, 1965; Schlanker 
& Leisman, 1969; Bateman, 1988), though de- 
velopmental constraints preclude adaptations typ- 
ical of truly prostrate growth (Bateman & Di- 
Michele, 1991). Nonetheless, these plants 
superficially resemble the basally branched archi- 
tectural model of Tomlinson (Hallé et al., 1978). 
Restriction of wood to the rhizomorph and highly 
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reduced stem defines the pseudoherbaceous habit 
sensu Bateman (1988). Rapid determinate growth 
and a strongly r-selected life strategy are inferred 
for these species, though whether their reproduc- 
tion was extended, or monocarpic upon cessation 
of growth, remains equivocal (Bateman, 1988; 
Bateman & DiMichele, 1991). 

Bauplan (3) also encompasses two subgroups. It 
characterizes Anabathra, Diaphorodendron sens. 
str., and Sigillaria, as well as the reconstructed 
adpression Bothrodendron punctatum (Wnuk, 
1989). All share Stone’s architectural model and 
the arboreous growth habit. Stone’s model resem- 
bles Schoute’s model in many parameters but ex- 
hibits clear differentiation between stem and lateral 
branches, which were produced throughout much 
of the life of the individual. Although a degree of 
dorsiventral flattening has been inferred for the 
lateral branches (Wnuk, 1985), their stable mor- 
phology and spiral phyllotaxy suggest persistent 
orthotropy sens. lat. rather than a transition to 
plagiotropy (cf. Hallé et al., 1978, table 7). De- 
viation from Stone’s model occurs during the final 
phase of determinate growth, occasional isotomous 
divisions (Hirmer, 1927; DiMichele & Phillips, 
1985; Wnuk, 1985), presumably heralding ex- 
haustion of the apical meristem. However, the ef- 
fect of such divisions is much less profound in 
bauplan (3) than in bauplan (2). Our assignment 
of Sigillaria to Stone’s model, which contradicts 
Hallé et al.’s (1978, fig. 71) assertion that the 
genus conforms to Schoute’s model, reflects our 
view that the stout cauline peduncles of Sigillaria 
(and the bauplan (2) genera) are homologous with 
entire lateral cone-bearing branches of Anabathra 
and Diaphorodendron; thus, by definition, Sigil- 
laria possesses lateral branches, though we distin- 
guish it as a separate architectural subgroup. 

This group of polycarpic plants possessed ex- 
current trunks and deciduous lateral branches. Wide 
cylinders of wood and periderm occur in the trunks 
of all the members of the group; they extend into 
the lateral branches of D. scleroticum, suggesting 
greater persistence (DiMichele, 1980, 1981, 1985). 
Growth and reproduction were both prolonged and 
sustained, conferring greater tolerance to extrinsic 
stress and allowing these species to occupy suitable 
habitats for considerable periods (DiMichele & 
Phillips, 1985; DiMichele et al., 1987). Relatively 
sporadic reproduction and apomixis (presumably 
facultative) have been reported in Sigillaria ap- 
proximata (Schopf, 1941; Phillips, 1979; Di- 
Michele & Phillips, 1985). 

DiMichele & Phillips (1985) argued that growth 
architecture and mode of reproduction largely de- 
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termine the ecological niches of specific arboreous 
lycopsids, a hypothesis that can be broadened to 
encompass nonarboreous arborescent species 
(Bateman, 1988; Bateman & DiMichele, 1991). 
Bateman & DiMichele (1991) further suggested 
that transitions in many of the characters describ- 
ing more detailed aspects of vegetative morphology 
reflect evolutionary changes in growth architec- 
ture. If so, it is especially important to assess the 
frequency and polarity of architectural changes 
during the history of the arborescent lycopsids. 

Our character analysis (Table 3) focused on 
specific homologous structures, whereas the five 
growth habits outlined above are polythetic sum- 
maries of several individual character states, some 
apomorphic and some plesiomorphic (C1-C7). For 
example, the arboreous habit corresponds with the 
apomorphic state of C1, the pseudoherbaceous habit 
with the plesiomorphic state of C2, and the pe- 
dunculate habit with the plesiomorphic state of C7; 
Schoute’s model with the apomorphic state of C3, 
Stone’s model with the apomorphic state of C6, 
and Corner’s model with the apomorphic state of 
C5. Despite their heterogeneity and partial depen- 
dence on the cladistic characters, there is consid- 
erable interpretative value in mapping the distri- 
butions of the growth habits across the preferred 
MPTs. The procedure is not wholly tautologous, 
as the preferred MPT for analysis B (habit char- 
acters omitted) is identical to that for the complete 
analysis (A; Fig. 9a), suggesting that habit char- 
acters had little direct effect on the topology of the 
latter. 

Two of the five growth habits (Corner’s model 
in Chaloneria, pedunculate Stone’s model in Sig- 
illaria) are autapomorphic at the generic level 
preventing assessment of their phylogenetic Sig- 
nificance. The preferred MPT for analysis A (Fig- 
6) requires homoplasy in at least two of the re 
maining three growth habits (pseudoherbaceou 
Schoute, arboreous Schoute, laterally branch 
Stone). It is equally parsimonious to assume a pseu- 
doherbaceous or arboreous hypothetical ancestor. 
In the first case, pseudoherbaceousness is replac 
by arboreousness immediately below Anabathr s 
with a reversal to pseudoherbaceousness in ee 
modendron. The transition to Stone’s model 5 
occurs immediately below Anabathra, with -i 
sumption of Schoute’s model in Diaphorodendron 
and the Hizemodendron—Lepidodendron-Lepr 
dophloios clade. The only sister groups that a 
equivocally possess the same habit are Paurode 
dron and Oxroadia at the base of the oe 
and Lepidodendron and Lepidophloios at the ape*- 
Moreover, there is no clear evolutionary t° 
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through the clade; Schoute’s model characterizes 
both the most primitive and most derived genera. 

We expected analysis D, based only on vege- 
tative characters, to provide a less homoplastic 
distribution of the major vegetative architectures. 
The most profound difference is the unification of 
the pseudoherbs (all attributed to Schoute’s model) 
as a basal paraphyletic group, thus depicting Hize- 
modendron as much more primitive than it appears 
in the full analysis (cf. Fig. 11 with Fig. 6). The 
topology for Figure 11, unlike that for Figure 6, 
requires recognition of the pseudoherbaceous habit 
as plesiomorphic. Although Anabathra is promoted 
to sister group of Diaphorodendron—Synchysi- 
dendron, a transition from Schoute’s to Stone’s 
model immediately above the pseudoherbs, fol- 
lowed by independent re-acquisition of Schoute’s 
model in the derived Synchysidendron and Lep- 
idodendron—Lepidophloios clades, remains the 
most parsimonious distribution of major architec- 
tures. 


Heterochrony. The strongly iterative occurrence 
of the tree habit among the arborescent lycopsids 
is mirrored in extant tree ferns and seed plants 
(Hallé & Oldeman, 1970; Hallé et al., 1978; White, 
1979; Funk, 1982; Tomlinson, 1983). We believe 
that radical changes in growth architecture are 
more readily achieved instantaneously, by mutation 
of genes involved in the production of the mor- 
Phogens that control early development (D-genes 
of Arthur, 1984, 1988), than gradually, by selec- 
tively driven adaptation. The resulting morpholog- 
ical change would appear instantaneous on a geo- 
logical time-scale. In this case, we believe that such 
saltational events were expressed as heterochrony 
Sens. lat.: a change in the timing of the appearance 
of a trait between ancestor and descendant (e.g., 
Gould, 1977; Alberch et al., 1979; Fink, 1982; 
McNamara, 1982; Rothwell, 1987; DiMichele & 
Bateman, 1989; Bateman & DiMichele, 1991). 
Many other correlated morphological and anatom- 
ical changes probably occurred immediately, as a 
result of epigenetic changes within the new bauplan 
and habit, or subsequently, as a result of adaptive 
honing by natural selection (Arthur, 1984; Bate- 
man & DiMichele, 1991). This scenario predicts 
= such macromutants very rarely generated evo- 
lutionary lineages, requiring a competition-free 
niche to allow establishment of the new population 
Prior to adaptive honing (Valentine, 1980; Arthur, 
1984, 1988; DiMichele et al., 1987). Hence, het- 
*rochronic anomalies must arise frequently to over- 
“ome their almost inevitable failure. 

Bateman & DiMichele (1991) argued that Hize- 
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modendron became pseudoherbaceous by hetero- 
chronic reduction from an arboreous ancestor. Pre- 
cocious division of the primary apical meristem 
minimized the length of the stem and prompted 
many subsequent character changes to accom- 
modate the new growth habit. Reduction in size 
and change in shape of the stem of Hizemodendron 
imply progenesis, a form of paedomorphosis (re- 
tention of ancestral characters in the descendant 
adult). 

In this paper, we are concerned less with the 
details of the postulated mechanism of vegetative 
reduction than its potential consequences for phy- 
logenetic reconstruction. The preferred MPT for 
all characters (Fig. 6) depicts Hizemodendron ser- 
ratum as primitive relative to its former congener, 
Lepidodendron hickii; together, the two genera 
constitute a paraphyletic sister group of Lepido- 
phloios. The preferred MPT for reproductive char- 
acters only (Fig. 12) shows an unresolved trichot- 
omy, thus allowing monophyly of Hizemodendron 
and Lepidodendron. In contrast, the preferred 
MPT for vegetative characters only (Fig. 11) de- 
picts a very distant relationship between these two 
genera. 

Although the phylogeny based on the greatest 
number of characters (Fig. 6) has the highest over- 
all probability of accuracy, it may be misleading 
in this case. If H. serratum evolved by progenesis 
directly from L. hickii (admittedly an improbable 
event, but a useful working hypothesis) and thereby 
lost all autapomorphies of L. hickii (H. serratum 
and L. hickii have no exclusive synapomorphies; 
Table 3), H. serratum would be depicted as prim- 
itive sister group of L. hickii. Losses of character 
states that the ancestor (L. hickii) shared with its 
former sister group (Diaphorodendron—Synchysi- 
dendron) are often depicted as reversals on the 
branch of the descendant (H. serratum) (this is 
only guaranteed if ACCTRAN optimization is used), 
whereas losses of character states not present in 
Diaphorodendron—Synchysidendron are not rec- 
ognized as character transitions at all; it is more 
parsimonious to assume that the ancestor of Hize- 
modendron never possessed them. Thus, hetero- 
chrony in general and progenesis in particular can 
undermine phylogenetic reconstruction, and mono- 
phyly of Hizemodendron and Lepidodendron re- 
mains a viable hypothesis. 

Bateman (1988, 1992) and Bateman & Di- 
Michele (1991) postulated a similar progenetic or- 
igin for the other two pseudoherbs, Oxroadia and 
Paurodendron. They are shown as the most prim- 
itive OTUs in all analyses (cf. Figs. 6, 11, 12), 
though they form a polychotomy with the remain- 
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ing bisporangiate-coned genera (Chaloneria and 
Anabathra) in Figure 12. As in the case of Hize- 
modendron, their primitiveness relative to OTUs 
with similar reproductive morphology is determined 
by vegetative characters. Anabathra is the only 
bisporangiate-coned tree included in our analysis 
and therefore provides the only potential arboreous 
ancestor for Oxroadia and Paurodendron. How- 
ever, the three OTUs differ in many spore char- 
acters, and the unbranched rhizomorph and su- 
perficial ligules of Paurodendron label the genus 
as relatively primitive or relatively derived, de- 
pending upon near-arbitrary polarization decisions. 
Other bisporangiate-coned trees, once reconstruct- 
ed, will provide more credible ancestors. 

A heterochronic origin for Oxroadia and Pau- 
rodendron would weaken our analysis, as an a 
priori assumption of their primitiveness was used 
to polarize most of the characters (i.e., they were 
used as partial outgroups). Inclusion in the data 
matrix of even more primitive OTUs may support 
our original assumption that Paurodendron and 
Oxroadia are sister groups to the remainder of the 
arborescent lycopsid clade. Our concern is largely 
driven by our opinion that the first arborescent 
lycopsid would have generated secondary tissues 
throughout its bauplan (the most simple develop- 
mental transition from inability to generate sec- 
ondary tissues), and that restriction of wood to 
certain modules reflects subsequent developmental 
modifications. Moreover, determinate growth and 
a centralized rhizomorphic rootstock are characters 
shared by all the OTUs, suggesting that they have 
exceptionally high burden (i.e., they play pivotal 
roles in the development and function of the or- 
ganism and influence other dependent characters: 
Riedl, 1979; Fortey & Jefferies, 1982; Donoghue, 
1989). These high-burden characters represent se- 
rious ontogenetic constraints to a truly prostrate 
growth habit (Bateman & DiMichele, 1991), and 
the bauplan appears much better adapted to up- 
right growth. Increase in body size to arboreous 
proportions may have occurred subsequently rath- 
er than concomitantly with acquisition of wood. 

Although Chaloneria is not a pseudoherb and 
is erect, the cladograms for all characters (F ig. 6) 
and vegetative characters only (Fig. 11) show that 
its unbranched, bilaterally symmetrical rhizome, 
unbranched stem, and cauline sporophylls (all aut- 
apomorphies in our analysis) are derived. This con- 
clusion is tempered by the possibility that Chalo- 
neria is not a genuine member of the ingroup, 
which would explain its numerous autapomorphies 
and its role as the greatest cause of topological 


instability in each analysis (Figs. 9, 10). The only 


well-supported conclusion from our study is that 
Chaloneria is more primitive than the most prim- 
itive member of the clade delimited by reduction 
to a single functional megaspore per megasporan- 
gium (i.e., than Diaphorodendron—Synchysiden- 
dron). 

The phylogenetic position of Chaloneria is es- 
pecially significant because it is the oldest recon- 
structed genus currently assigned to the Isoetales 
(Pigg & Rothwell, 1983a; Rothwell & Erwin, 1984). 
A sister-group relationship with a widely recognized 
lepidodendralean genus such as Sigillaria (e.g. 
Fig. 9k) would imply paraphyly of the Lepidoden- 
drales and support Meyen’s (1987: 70-81) deci- 
sion to synonymize the Lepidodendrales into the 
Isoetales. Further resolution of these problems re- 
quires a broader cladistic analysis that includes 
other bona fide isoetaleans (including Isoetes), po- 
tential arboreous ancestors (e.g., Lepidodendrop- 
sis—Protostigmaria: Jennings, 1975; Jennings et 
al., 1983), and ostensibly more primitive OTUs 
(e.g., Selaginella) (see Bateman, 1992). 


Outgroups and ancestors. Thus, we return to the 
fundamental questions that prompted this study. 
What character states delimit the Lepidoden- 
drales? Is the group monophyletic? If so, what i$ 
the most appropriate outgroup? The chosen answer 
to this question leads to an even more loaded ques- 
tion: What is the most probable ancestor of the 
ingroup? It also largely determines percepton of 
the phylogenetic relationships among the ingroup 
members. 

We believe that the greatest weakness of our 
analysis is the narrow temporal and ecological range 
represented by our OTUs; most of the species arè 
restricted to at most the ca. 10 Ma of the West 
phalian (Fig. 3) and to the coal swamps of Euramer- 
ica. However, the main phylogenetic groups wis 
the Lepidodendrales (or at least species possessing 
many of their diagnostic character states; whole- 
plant reconstructions have not yet been ee 
for pre-Westphalian arboreous lycopsids) can 
traced back at least another 20 Ma, to the pers 
Moreover, reproductive organs consistent with 
most apomorphic genus, Lepidophloios, have ee 
recovered from Ivorian strata, a further 15 s 
older (Fig. 3, inset; Long, 1968). This implies th 
all of the sister groups of this genus had se 
by the Ivorian; unfortunately, only one of our OT A 
(Oxroadia gracilis-Oxroadia sp. nov.) was ay 
constructed from such early assemblages. AS Y n 
incompletely reconstructed arboreous lycopsids es 
necessarily bona fide lepidodendraleans) were Wi i 
spread and at least locally ecologically dominan 
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by the latest Devonian (e.g., Scheckler, 1986a, b; 
DiMichele et al., 1992). 

Thus, the combination of the stratigraphic rec- 
ord and our phylogeny suggests that the three main 
groups of monosporangiate-coned genera recog- 
nized in our cladistic analysis (Sigillaria, Diapho- 
rodendron—Synchysidendron, Hizemodendron— 
Lepidodendron—Lepidophloios) diverged at least 
35 Ma prior to the Westphalian coal-swamp lager- 
statten that provided most of our OTUs. This would 
explain why these groups show similarly large de- 
grees of divergence from their putative bisporan- 
giate-coned ancestor(s) (Fig. 14) and are supported 
by many character-state transitions (Fig. 6). Al- 
ternatively, the saltational evolutionary scenario 
erected for major vegetative changes may be ex- 
tended to encompass reproductive innovations, 
eliminating the need for intermediate taxa during 
the early radiation of the group. 

Current evidence suggests that our morpholog- 
ically divergent OTUs together exhibit most of the 
character states possessed by the arborescent ly- 
copsids as a whole, including other Pennsylvanian 
Species and their Mississippian and Devonian an- 
tecedents. However, the paucity of genera in our 
analysis, and the fact that most represent only the 
final period of the history of the group, implies that 
we have sampled only a restricted range of the 
combinations of character states that existed. This 
would explain the large number of character states 
that occur as genus-level autapomorphies in our 
cladogram (Fig. 6), leaving few character states to 
Support the consequently weak links that constitute 
the main axis of the cladogram and determine 
perceived genus-level relationships. Inclusion of 
older OTUs, dating back to the main radiation of 
the group, would probably alleviate this problem 
by transforming genus-level autapomorphies into 
fenus-level synapomorphies. In an alternative less 
Sradualistic scenario, the large number of genus- 
level autapomorphies may reflect evolutionary dy- 
namics, particularly the simultaneous origin of 
blocks of characters linked by pleiotropic or epi- 
Benetic factors (e.g., Levinton, 1988). This mode 
of evolution may be difficult to resolve cladistically 
for a variety of methodological reasons (these will 

discussed in a future paper; see also Lemen & 
Freeman, 1989), 

Older OTUs are also needed to determine con- 
vincingly whether the Lepidodendrales are mono- 
Phyletic and in particular to provide more satis- 
factory outgroups. However, before these questions 
ĉan be addressed, the character states that sup- 
Posedly delimit the Lepidodendrales should be re- 


"ewed. Four are most commonly cited: the pos- 


session of rootlet-bearing rhizomorphs, secondary 
tissues (wood and periderm), ligules, and hetero- 
spory (Chaloner, 1967; Stewart, 1983). These 
structures provided the holapomorphies that unite 
all the OTUs included in our analysis (Fig. 6), but 
most (possibly all) have a greater level of univer- 
sality. For example, a wide range of enigmatic latest 
Devonian and earliest Mississippian lycopsids pos- 
sessed wood (Meyer-Berthaud, 1981, 1984; 
Scheckler, 1986a, b; Matten, 1989; Roy & Mat- 
ten, 1989). Rhizomorph-like, rootstocks, ligules, 
and heterospory all characterize homophyllous Se- 
laginella, the type genus of the Selaginellales (e.g., 
Bierhorst, 1971; Bold et al., 1980). Moreover, 
Paurodendron (and therefore, by implication, its 
sister genus Oxroadia) was assigned by Schlanker 
& Leisman (1969) to Selaginella, and it is widely 
accepted as a member of the Selaginellales (e.g., 
Taylor, 1981; Stewart, 1983; Meyen, 1987). On 
these criteria, the Lepidodendrales could be cir- 
cumscribed to include both Paurodendron and 
homophyllous Selaginella. 

Moving progressively up the clade, the next 
OTU encountered is the most primitive tree, An- 
abathra. Regarding this OTU as the most primitive 
lepidodendralean would allow delimitation of the 
order using the arboreous habit and associated 
modifications of stelar anatomy, together with foliar 
parichnos. Unfortunately, many of these character 
states are homoplastic as a result of loss during the 
hypothesized progenetic evolution of pseudoherbs 
from trees: Hizemodendron from a Lepidoden- 
dron-like ancestor, and possibly Oxroadia and/or 
Paurodendron from Anabathra-like ancestor(s). 
Although phylogenetically valuable, these char- 
acter states are not ubiquitous within the clade. 
Many workers would argue that the clade is delim- 
ited primarily by possession of a stigmarian rhi- 
zomorph. However, we were unable to identify any 
profound characters that distinguish the stigmarian 
rhizomorph of Anabathra from the supposedly 
nonstigmarian rhizomorph of Oxroadia, which is 
much smaller and more compact but otherwise very 
similar. Also, this clade contains Chaloneria, an- 
other nonarboreous OTU. Chaloneria possesses 
several autapomorphies, notably an unbranched 
stem lacking cones and a bilaterally symmetrical 
rhizomorph, that suggest affinities with the extant 
genus Isoetes. If Pigg & Rothwell (1983a) cor- 
rectly ascribed Chaloneria to the Isoetales, and if 
the genus is correctly positioned in our phylogeny 
(which is by no means certain; Fig. 9), inclusion 
of Anabathra in the Lepidodendrales and contin- 
ued recognition of the Isoetales would render the 
former order paraphyletic (Appendix 1A). 
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Perhaps the most cohesive clade includes Sig- 
illaria as its most primitive genus. It is delimited 
by the nonhomoplastic synapomorphies of leaf 
cushions and monosporangiate cones that generate 
megasporangium-megasporophyll disseminules. We 
are confident of the monophyly of the three main 
groups of OTUs that constitute the clade (Sigil- 
laria, Diaphorodendron—Synchysidendron, Hi- 
zemodendron—Lepidodendron—Lepidophloios), 
and our parsimony analysis strongly supports 
monophyly of the clade as a whole (nevertheless, 
we note that such an adaptively valuable suite of 
character states could reflect parallel responses in 
two or more lineages to similar selective regimes, 
thus confounding parsimony: cf. Coddington, 1988). 
Given that our assumption of homology among the 
three lineages in monosporangiate cones and the 
megasporangium-megasporophyll complex as dis- 
seminule is the crux of the preferred MPT (Fig. 
6), these characters merit even more careful scru- 
tiny. For now, we refer to this clade, more narrowly 
defined than most perceptions of the Lepidoden- 
drales, as the “‘Segregationists” (referring to the 
segregation of megasporangia and microsporangia 
in different cones). Members of the less inclusive 
clade that excludes Sigillaria and is delimited pri- 
marily by reduction to a single functional mega- 
spore that germinates within the sporangium are 
the “‘Isolationists.” 


CONCLUSIONS 


Empirical observations. We are confident that 
each of the 10 genera analyzed by us is mono- 
phyletic. This conclusion is not especially profound, 
as six of the genera are here represented by only 
one species and thus not cladistically testable (Fig. 
8). The analysis prompted segregation of two new 
genera: Hizemodendron from Lepidodendron 
(Bateman & DiMichele, 1991), and Synchysiden- 
dron from Diaphorodendron (DiMichele & Bate- 
man, 1992). However, these decisions were taken 
primarily on the grounds of differences in several 
characters (many directly or indirectly reflecting 
different growth architectures) rather than as at- 
tempts to disaggregate para- or polyphyletic groups. 
Diaphorodendron-—Synchysidendron is undoubt- 
edly monophyletic; Hizemodendron—Lepidoden- 
dron is depicted in Figure 6 as paraphyletic but, 
as we have argued, may nonetheless be monophy- 
letic. Derived (monosporangiate-coned) genera 
constitute three distinct monophyletic clades that 
are most appropriately regarded as families: Sig- 
illariaceae (Sigillaria), Diaphorodendraceae (Dia- 


phorodendron—Synchysidendron), and Lepido- 
dendraceae (Lepidophloios-Lepidodendron- 
Hizemodendron). Each family is supported by sev- 
eral character state transitions (Fig. 6). 

Evidence for the postulated monophyly of the 
monosporangiate-coned clade, and for the primi- 
tiveness of the Sigillariaceae relative to the Dia- 
phorodendraceae and Lepidodendraceae, is more 
equivocal. This largely reflects our inability to make 
confident statements concerning phylogenetic re- 
lationships among the four primitive, bisporangiate- 
coned genera. They form a highly heterogeneous, 
paraphyletic (or possibly polyphyletic) plexus of 
disparate morphologies that share a free-sporing 
mode of reproduction. Oxroadia and Pauroden- 
dron differ in many characters and are united in 
Figure 6 by arguably only one synapomorphy; we 
doubt their apparent monophyletic status. Only 
Anabathra is a tree; the pseudoherbs Oxroadia 
and Paurodendron and possibly even the shrub- 
sized phallos Chaloneria are potential progenetic 
descendants of trees broadly similar to, but prob- 
ably distinct from, Anabathra. 

Determining the origin(s) of the monosporan- 
giate-coned clade will require inclusion of pre- 
Pennsylvanian monosporangiate-coned species and 
a broader selection of bisporangiate-coned trees; 
several potential candidates, all requiring further 
reconstruction before they can be used with con- 
fidence in cladistic analyses, are listed in Table 5. 
Determining the origin(s) of the arboreous lycop- 
sids, and of putatively progenetic bisporangial® 
coned genera such as Oxroadia and Pauroden- 
dron, also necessitates inclusion of more distantly 
related nonarborescent lycopsids (e.g-, Selaginella 
sens. lat., Leclercgia) to reassess character suis 
polarities. Given these observations, we envisage 
eventual redelimitation, or possibly amalgamatıon, 
of three widely recognized lycopsid orders (Appen 
dix 1A; see also Bateman, 1992). 

Absence from the present analysis of nonarbo- 
rescent species and of any credible ancestor of the 
most primitive arboreous genus (Anabathra) t% 
gether prevent determination of whether arbores- 
cence (secondary thickening) and arboreousness 
(large body and upright growth) evolved ee 
neously or sequentially. Other especially gai 
and potentially linked innovations that posse i 
higher level of universality than our ingroup a 
determinate growth and the centralized ae 
morphic rootstock. In general, the most significan 
evolutionary advances within the bisporangia'® 
coned plexus appear to have involved vegetatiy 
rather than reproductive organs, indicating 
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TABLE 5. Selected genera that are potentially phylogenetically informative but are currently insufficiently known 
to provide satisfactory cladistic data sets (listed in order of appearance in the stratigraphic record). We note that 
‘Lepidodendron’ calamopsoides is not closely related to Lepidodendron sens. str. (R. M. Bateman, unpublished 
obs.) and that Bothrodendron sens. lat. is probably a polyphyletic aggregate of several disparate species (Scott, 1920; 
Thomas & Meyen, 1984). See also reviews by Chaloner (1967), Meyer-Berthaud (1981, 1984), and Matten (1989). 


Lepidosigillaria whitei Krausel & Weyland 


White (1907); Arnold (1947); Grierson & Banks (1963) 


[Late Givetian-Early Famennian: New York State] 
Protolepidodendropsis spp. Gothan & Zimmermann 
Høeg (1942); Schweitzer (1965) 
[Famennian: Euramerica] 
Trabicaulis spp. Meyer-Berthaud 


Meyer-Berthaud (1981, 1984); Roy & Matten (1989) 


[Famennian-Early Tournaisian: south-central France; New York State] 


Cyclostigma kiltorkense Haughton 


Johnson (1913); Chaloner (1967, 1968, 1984); Chaloner & Meyer-Berthaud (1983) 


[Strunian: southwestern Ireland] 
Landeyrodendron spp. Meyer-Berthaud 
Meyer-Berthaud (1981, 1984) 
[Early Tournaisian: south-central France] 


Lepidodendropsis spp. Lutz~Protostigmaria eggertiana Jennings 
Lutz (1933); Turina & Lemoigne (1975); Jennings (1975); Jennings et al. (1983) 


[Tournaisian: Euramerica] 
Valmeyerodendron triangularifolium Jennings 
Jennings (1972) 
[?Tournaisian: Illinois] 
Bothrodendron spp. Lindley & Hutton 


Scott (1908); Weiss (1908); Calder (1933b); Stubblefield & Rothwell (1981); Wnuk (1989) 


[Tournaisian- Westphalian: Euramerica] 
‘Lepidodendron’ (?Anabathra) calamopsoides Long 
Long (1964, 1971, 1986); Scott & Galtier (1988) 
[Late Tournaisian: northern Britain] 
Levicaulis arranensis Beck 
Beck (1958); Pant & Walton (1961) 
[Mid-Visean: southwestern Scotland] 
‘Lepidodendron’ (?Phytokneme) brownii Unger 


Chodat (1911); Calder (1933a); Meyer-Berthaud (1981) 


[Visean: southern Scotland] 


they are economic adaptations (sensu Eldredge, 
1989) employed continuously in competition for 
resources. The cladograms imply that the well- 
documented sequence of reproductive innovations 
in the monosporangiate-coned clade, which ulti- 
mately led to seed analogs in Lepidophloios (e.g. 
Phillips, 1979), occurred later, though they may 
have been developmentally linked to additional veg- 
*tative modifications (this hypothesis requires fur- 
ther study). Interestingly, the equally well-docu- 
mented, climatically driven end-Westphalian 
extinctions of specific elements of the coal-swamp 
floras (Phillips et al., 1977, 1985; Phillips & Pep- 
pers, 1984) most seriously affected the most de- 
"ved portion of the arborescent lycopsid clade, 
eliminating the “‘Isolationist’’ families Diaphoro- 


dendraceae and Lepidodendraceae (Fig. 3). Such 
extinctions of major monophyletic groups are un- 
common (Smith & Patterson, 1988) and require 
a causal rather than a purely stochastic explana- 
tion. 

These analyses provide useful (if circumstantial) 
evidence for the relative burden of particular types 
of character, in the guise of amounts of homoplasy. 
Interestingly, this partly reflects the physical scale 
(dimensions) of the feature represented by the char- 
acter relative to that of the plant body. Characters 
of largest scale (notably overall growth habit) and 
smallest scale (e.g., various details of cellular his- 
tology and spore ornamentation) are generally more 
homoplastic than those of intermediate scale, such 
as stelar and associated trace morphology, the basic 
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structure of leaf bases and sporophylls, and the 
nature of the dispersal unit. As a broad general- 
ization, large- and small-scale features delimit spe- 
cies and genera, intermediate scale features delimit 
families or still higher taxa. 

The results of our study will not encourage pro- 
ponents of organ phylogenies. Merely bisecting our 
data matrix into submatrices representing vege- 
tative (Fig. 11) and reproductive (Fig. 12) organs 
generated substantially different preferred MPTs 
that were clearly inferior to the preferred MPT of 
the full matrix (Fig. 6). The analysis of reproductive 
organs could not satisfactorily resolve the relation- 
ships among the four most primitive and three most 
derived genera. The analysis of vegetative organs 
misplaced Hizemodendron as unduly primitive and 
Sigillaria as unduly derived and could only dis- 
tinguish Synchysidendron from Diaphoroden- 
dron using cone axis characters that mirror those 
of ultimate vegetative axes. Nonetheless, the ar- 
boreous members of the three most derived families 
(Sigillariaceae, Diaphorodendraceae, Lepidoden- 
draceae) persist as clades in the analyses of both 
submatrices. This shows that the families as cur- 
rently known can be approximately delimited using 
either vegetative or reproductive characters alone, 
even if their relationships cannot be determined 
accurately. 

Our results are even less encouraging for clas- 
sifications based on even more reduced suites of 
characters. We have identified homoplasy in many 
supposedly diagnostic character states, including 
some of the leaf-base details that are traditionally 
used to classify adpressed lycopsid axes. On the 
basis of these observations, we support in principle 
the hierarchical system of well-known core taxa 
and less well-known satellite taxa proposed for the 
Lycopsida by Thomas & Brack-Hanes (1984), but 
are convinced that reconstructed, anatomically 
preserved whole plants provide better core taxa 
than the reproductive organs favored by Thomas 
& Brack-Hanes. Certainly, whole-plant reconstruc- 
tions are essential prerequisites for convincing phy- 
logenetic and ecomorphic interpretations. 


Methodological observations. We do not regard 
phylogenetic reconstruction as an isolated, objec- 
tive procedure divorced from hypotheses of evo- 
lutionary mechanisms; rather, it is positive feed- 
back between the two sets of paradigms that leads 
to greater understanding. The evolutionary history 
of the arborescent lycopsids is not a simple story 
of progressively increasing complexity expressed 
throughout the bauplan. In particular, high-burden 
characters such as determinate growth and a cen- 


tralized rhizomorphic rootstock, together with the 
small number of module types that constitute the 
bauplan, predisposed the plants to profound het- 
erochronic changes in body size and body plan; we 
believe that these were manifested as geologically 
instantaneous events resulting from changes in de- 
velopmental regulation. This saltational evolution- 
ary scenario has considerable predictive value, par- 
ticularly if considered in tandem with advances in 
understanding of the ontogeny, functional mor- 
phology and physiology of these remarkable or- 
ganisms (e.g., Phillips & DiMichele, 1992). More- 
over, saltational scenarios can be falsified (or at 
least highly modified) by cladograms, if the pre- 
dicted positively correlated suite of character state 
transitions is dissociated (e.g., by the interpolation 
of additional OTUs onto the internode in question: 
R. M. Bateman, in prep.). 

We have been unable to discern any substantive 
differences between reconstructing the morpholog- 
ical phylogenies of extinct and extant species (cf. 
Stein, 1987; Gauthier et al., 1988; Donoghue et 
al., 1989; Boy, 1990). The inevitable typological 
nature of conceptual whole-plant fossils is not det- 
rimental in the essentially typological realm of cla- 
distics. Our 16 OTUs undoubtedly represent a 
highly rarified sample of all the arboreous lycopsid 
species that ever existed. This contributed to wn 
eral problems, notably the broad morphological 
discontinuities separating some clades (Fig. 14) and 
the absence of satisfactory outgroups. However, 
cladistic analyses based exclusively on extant sp®- 
cies are even more selective; unique character com 
binations found only in the fossil record, especially 
during the initial radiation of a major clade, are 
deliberately excluded. Similarly, opportunities to 
use stratigraphic-temporal evidence to assist po 
larization and characters, and (more importantly) 
to select among alternative topologies generat 
from the same data matrix, are squandered. The 
question of excluding fossils does not arise mM the 
case of the monosporangiate-coned lycopsid clade, 
which apparently lacks extant descendants. Despite 
the serious problems posed by incomplete prese™ 
vation in the fossil record, we were able to score 
a large number of characters representing all or- 
gans of our OTUs and generated a large date 
matrix containing only a small proportion of H 
values. Although technically feasible, inclusion 3 
less well-known, partially reconstructed ha 
should be postponed pending further investigat! ' 
of the effects of missing values on tree-le 
culation and character state optimization. all 

This study indicates that well-understood fi 
are as valuable for phylogenetic studies as 3Y 
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extant organisms and do not require any special 
methodological concessions. We note that profound 
gaps in the fossil record caused by nonpreservation 
are matched by equally profound gaps in extant 
floras, namely those resulting from extinction of 
their precursors. Access to the time dimension, and 
thereby to character states (and combinations of 
character states) no longer in existence, justifies a 
pivotal role for paleontological data in phylogenetic 
reconstruction. 
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APPENDIX l. NOMENCLATURAL AND TAXONOMIC NOTES 


A. Higher taxa. Bateman (1990b) recommended the 
supra-ordinal classification of Knoll & Rothwell (1981) 
and ordinal classification of Stewart (1983); together, 
these classifications imply monophyly of the Division Tra- 
cheophyta and Class Lycopsida, but present a provisional, 
egalitarian (grade) arrangement of six orders within the 
Lycopsida (Drepanophycales, Protolepidodendrales, Ly- 
copodiales, Selaginellales, Lepidodendrales, Isoetales). Ac- 
quisition of further phylogenetic information will require 
re-delimitation of these orders and their rearrangement 
into a more hierarchical classification (Bateman, 1990b, 
1992; Hueber, 1992). Colloquial (informal) names are 
consistently rooted in their formal counterparts; thus, 
“lycopsid” is used for the Class Lycopsida and “‘lepidoden- 
dralean” for the Order Lepidodendrales (Bateman, 1990b). 

At present, the Lycopsida are perceived as being de- 
limited by exarch protostele, scalariform metaxylem with 
Williamson fimbrils, vascularized “microphylls,” and fo- 
liar /axillary eusporangia (e.g., Stewart, 1983). With the 
possible exception of stigmarian rhizomorphs, none of the 
characters traditionally used to delimit the Lepidoden- 


drales (secondary thickening, periderm, ligules, hetero- 
spory) is confined to the order. The present analysis shows 
that the Lepidodendrales as currently delimited may not 
be monophyletic, and even if monophyletic may not be 
the most meaningful clade to use as the basis of an ordinal 
level classification (see also Matten, 1989). Together with 
the Isoetales and Selaginellales, the Lepidodendrales re- 
quire redelimitation or amalgamation. 


B. Anabathra versus Paralycopodites. We accept Pear- 
son’s (1986) arguments that Paralycopodites brevifolius 
(Williamson) DiMichele (formerly ‘Lepidodendron’ brevi- 
folium Williamson) is identical to, and a junior synonym 
of, Anabathra pulcherrima Witham (cf. Witham, 1833; 
Williamson, 1872). We also recognize that the adpression 
genus Ulodendron Lindley & Hutton (1831) both resem- 
bles and nomenclaturally pre-dates Anabathra, but sus- 
pect that Ulodendron is considerably more inclusive (i.e., 
more broadly delimited) than Anabathra. 

Unfortunately, A. pulcherrima is a form-species (sensu 
Bateman & Rothwell, 1990), having been correlated with 
several cone-species of the plesiomorphic bisporangiate 
genus Flemingites (DiMichele, 1980; Pearson, 1986). 
The type material of 4. pulcherrima, from the Late 
Tournaisian of Allanbank, southeastern Scotland, co-0c- 
curred with Lagenicula horrida megaspores (these were 
incorrectly referred to L. subpilosa by Pearson, 1986: 
K. M. Bartram, pers. comm. 1987; H. L. Pearson, pers. 
comm. 1987), which have been found in Flemingites 
gracilis cones (Chaloner, 1953; Brack-Hanes & Thomas, 
1983). Elsewhere, A. pulcherrima co-occurs with F. scotti 
at the Late Visean locality of Pettycur, southeastern Scot- 
land (e.g., Williamson, 1872; Jongmans, 1930; Di- 
Michele, 1980; Scott et al., 1984; Pearson, 1986), and 
with either F. diversus or F. schopfii at many Westphalian 
coal-ball localities in Euramerica (DiMichele, 1980). The 
differences between these cone-species are subtle and bed 
readily resolved cladistically. In particular, our scoring 
A. pulcherrima focused on Pennsylvanian rather i 
Mississippian assemblages; F. diversus (Wesyihan a 
Felix, 1954) and F. schopfii (Westphalian B: Brach 
1970) differ primarily in the mean number of pac we BP 
per megasporangium and thus were not differentiat 
our data matrix (Table 3, Fig. 3). aca 

We anticipate taxonomic revision of Mere 
reproductive organs of Anabathra as increased know RA 
(particularly of Mississippian forms) allows further kres 
plant reconstructions; for example, ‘Lepidode age 
calamopsoides Long (1964), which co-occurs wit = 
marian rootstocks, Flemingites-like strobili and pea $ 
cula crassiaculeata megaspores (Scott & Mer 
thaud, 1985; Scott & Galtier, 1988; Scott, 1990), may 
be a species of Anabathra (cf. Table 2). Such plants 
form the basis of a new arborescent lycopsid family, 
Anabathraceae. 


C. Erection of new genera and species. We have mni 
erately avoided formal reclassification of the 17 W a 
plant species included in our study. Thus, the pai Syn- 
species (one each of Oxroadia, Sigillaria, 992: 
chysidendron) remain unnamed (but see Bateman, 
DiMichele & Bateman, 1992). us to 
However, the results of this study encouraged 4 by 
restrict further the range of variation encomp 
Lepidodendron sens. str., a process that was begun de- 
DiMichele (1981, 1983, 1985) when clarifying hloios, 
limiting parameters of Lepidodendron, Lepidop 
and Diaphorodendron sens. lat. In order to © 
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Lepidodendron from an apparently paraphyletic (Fig. 6) 
to a monophyletic entity, we retained only the anatomi- 
cally preserved equivalent (L. hickii) of the type species 
(L. aculeatum) and erected ‘Lepidodendron’ serratum 
as the type species of a new genus, Hizemodendron 
(Bateman & DiMichele, 1991). The two genera share 
similar reproductive organs, but differ in many vegetative 
characters; at least most of these differences may reflect 
the imposition of radically different growth habits on a 
shared bauplan. 

Although our study strongly supports monophyly for 
Diaphorodendron sens. lat. (i.e., sensu DiMichele, 1985; 
see Fig. 6), the precedent of generic distinction of species 
sharing similar reproductive organs but exhibiting major 
differences in growth habit and ontogeny requires the 
recognition of the two most apomorphic species, ‘D.’ 
dicentricum and ‘Diaphorodendron’ sp. nov., as a new 
genus of arboreous lycopsid, Synchysidendron. In con- 
trast with the more plesiomorphic Diaphorodendron (epit- 
omized by the type species, D. vasculare), S. dicentricum 
and Synchysidendron sp. nov. lack lateral branches and 
were probably monocarpic (DiMichele, 1981, 1985; Bate- 
man & DiMichele, 1991; DiMichele & Bateman, 1992), 
thus possessing the same bauplan, growth habit, and re- 
productive strategy as Lepidodendron sens. str. (Fig. 1). 


D. Revision of the cone-genera. As currently delimited, 
cone form-genera serve as shorthand for co-occurring 
complexes of character states. A few cone-genera are 
assignable to single stem-genera (Table 2), notably Ma- 
zocarpon to Sigillaria (e.g., Schopf, 1941; Feng & Roth- 
well, 1989) and Lepidocarpon to Lepidophloios (Di- 
Michele, 1983). The microsporangiate genus 
Lepidostrobus characterizes Hizemodendron, Lepido- 
dendron, and Lepidophloios; cone-species of each of 
these stem-genera can only be distinguished by continuous 
quantitative characters and microspore morphology (e.g., 
Willard, 1989a). In contrast, Achlamydocarpon is a 
greatly inflated form-genus (e.g., Leisman & Phillips, 1979) 
encompassing three main morphological groups: (1) mega- 


sporangiate A. takhtajanii-type, (2) megasporangiate A. 
varius-type, and (3) microsporangiate A. varius-type (Ta- 
ble 2). Groups (2) and (3) are restricted to Diaphoro- 
dendron—Synchysidendron; cones of all three Diapho- 
rodendron species and both Synchysidendron species can 
only be distinguished by the medullated steles of the latter, 
despite major differences in growth architecture and veg- 
etative anatomy between the two genera (DiMichele, 1981). 
Group (1) cones characterize Hizemodendron and Lep- 
idodendron (DiMichele, 1983; Bateman & DiMichele, 
1991). We believe that megasporangiate A. takhtajanii 
and A. varius are sufficiently distinct to merit generic 
distinction (the latter would require a new organ-genus). 
Assignment of megasporangiate and microsporangiate A. 
varius to different cone-genera would be more consistent 
with the systematic treatment of Hizemodendron, Lep- 
idodendron, and Lepidophloios cones in Table 2, though 
Phillips (1979: 256, 259) presented several arguments 
against this option. 


APPENDIX 2. ANALYTICAL ADVANCES 


Given the appropriate microcomputers and software, 
more elegant solutions are now available to some of the 
difficulties that we encountered when performing these 
analyses in 1989. For example, the problem of storage 
of only superficially different topologies, resulting from 
polychotomies, has been solved in Version 3.0 of PAUP 
(Swofford, 1991). There is much to commend an ana- 
lytical approach that entails initial parsimony analysis in 
PAUP 3.0, subsequent comparison of MPTs with those 
generated by using Hennig86 Version 1.6 (Farris, 1989), 
and printing of interesting topologies and character state 
distributions using MacClade Version 3.0 (Maddison & 
Maddison, 1991). Estimation of degrees of support for 
particular nodes using bootstrapping (Efron, 1982; Fel- 
senstein, 1985; Sanderson, 1989) is gaining in popularity, 
though there is no statistical substitute for detailed ex- 
amination of suboptimal-length topologies. 


COMPARATIVE ECOLOGY AND Tom L. Phillips? and William A. DiMichele’ 
LIFE-HISTORY BIOLOGY OF 

ARBORESCENT LYCOPSIDS 

IN LATE CARBONIFEROUS 

SWAMPS OF EURAMERICA! 


ABSTRACT 


The comparative ecologies of Diaphorodendron, Lepidodendron, Lepidophloios, Paralycopodites (= Anabathra), 
and Sigillaria in Late Carboniferous coal swamps serve as a context for assessing life cycles and exploring possible 
structure-function relations. The distinctive aspects of the “‘lycopsid tree habit” in lepidodendrids are emphasized as 
part of the arborescent reproductive architecture of relatively short-lived (10-15 years) plants. These include: 
determinate apical growth of the aerial shoot system and the anchoring stigmarian system, which have marked 
homologies; limited amount of secondary xylem, lack of secondary phloem, and periderm as the major support tissue; 
differentiated gas-diffusion system associated with appendages, including lacunae and parichnos; abscission of ap- 
pendages and lateral branches; retention of leaf cushions; and close relationships between mode and timing of branching, 
to cone-bearing and heterosporous reproductive biologies. Vegetative structure-functions explored include the possibility 
that lacunae and parichnos were involved in internal mediation of gas diffusion as opposed to aeration functions. The 
possibility exists that parts of the stigmarian system were involved in CO, acquisition from substrates; some may have 
been photosynthetic. These functions are considered in the context of the light sharing and diffuse photosynthesis 
evident in the pole architecture. The combination of such possibilities is related, in part, to the xeromorphie char- 
acteristics of the arborescent habit, raising the question about a modified kind of C, photosynthesis such as CAM 
(Crassulacean Acid Metabolism). Stigmarian lycopsids dominated tropical Westphalian coal swamps as an array of 
genera with relative distributional abundances reflecting habitat partitioning according to edaphic conditions, including 
temporal disturbance patterns. Species appear to exhibit different levels of tolerance to disturbances and range from 
colonizers to site occupiers. Paralycopodites, with prolific, free sporing, bisporangiate cones, was most abundant in 
frequently disturbed, partially exposed, peat- to mineral-rich habitats (ecotonal). Monosporangiate Lepidophlotos and 
Lepidodendron were associated typically with deeper, standing-water habitats, and, in association with te 
branching, monocarpically produced specialized monosporic megasporangium-sporophyll units, termed aquacarps 
Diaphorodendron species were monosporangiate with aquacarps and range from a typically persistent, low level 
reproductive output on deciduous lateral branches, to monocarpy (D. dicentricum) with terminal branching. Sigillaria 
was less closely associated with peat swamps, as a sporadic occupant associated with major disturbances, such as 
flood/dry down cycles. Whorls of monosporangiate cones were produced intermittently, perhaps in seasonally wel 
dry conditions, with megaspore-sporangial dispersal units derived from cone fragmentation. In North America Sigillaria 
was the principal lepidodendrid survivor of the swamp extinctions near the Middle-Upper Pennsylvania boundary. 
The ecological roles of stigmarian lycopsids in Westphalian coal swamps are distinctive as an array of heterosporous 
reproductive architectures that were collectively dominant. Plants were characteristically much taller than other trees, 
yet did not shade out lower vegetation. They also were a major stabilizing influence on substrates with their extensive, 
baffling and anchoring systems in the high disturbance and abiotically stressed environments of peat swamps: © 
environmental circumstances of the first major coal age appear to have selected against long-lived or slow-growing 
trees in most coal swamps. Lepidodendrids constituted the most important of the arborescent genera both because © 
their unusual array of reproductive biology in such large structural and yet short-lived growth habits, as well as many 
physiological attributes that are only partially known or conjectured. 
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Trees referred to as lepidodendrids or stigmarian ically from arborescent forms of other non-lycopod 
lycopsids encompass lycopsid zeniths in size, struc- lineages as to merit distinction of a “lycopsid tee 
tural complexity, and heterosporous diversification habit.” The principal lepidodendrid genera of the 
as well as ecological dominance in Late Carbonif- Late Carboniferous coal swamps, known from an- 
erous tropical swamps. Lepidodendrid trees were atomical preservation, include Diaphorodendro™ 
so different anatomically, functionally, and ecolog- Lepidodendron, Lepidophloios, Paralycopodites 
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(= Anabathra), and Sigillaria. Anatomical and 
developmental studies drawn from these genera 
provide a composite picture of the lycopsid tree 
habit that emphasizes some of the differences from 
conventional trees. 

Lepidodendrids were apically determinate in both 
the aerial shoot system (Andrews & Murdy, 1958; 
Eggert, 1961; Lemoigne, 1966; Chaloner & Mey- 
er-Berthaud, 1983) and apparently also in the an- 
choring stigmarian system (Frankenberg & Eggert, 
1969); the anatomical transition from stem to rhi- 
zomorph occurred in the unbranched base of the 
trunk. Marked homologies between stem and rhi- 
zomorph systems, both in axial organization and in 
appendage anatomy and arrangement, indicate that 
the stigmarian rhizomorph was, in part, function- 
ally but not morphologically a “root system.” Its 
anchorage and support of the pole-type trunk was 
provided by dichotomous branchings and the ex- 
tensive appendages they bore. There were no root 
caps or root hairs. The apex of the growing stig- 
marian axes tapered abruptly to a rimlike apical 
roove associated with appendage origins; this was 
terminated with a protective plug of apparently 
parenchymatous tissue (Rothwell, 1984). The cy- 
lindrical appendages were largely air filled, appar- 
ently buoyant, with a large external surface to 
biomass ratio. They radiated for lengths of 0.5 m 
or more; some are known to have been once di- 
chotomous in the sigillarian stigmarias (Lemoigne, 
1963). 

The aerial shoot system exhibited structural sup- 
port mainly by nonwoody tissues, especially peri- 
derm, formed beneath the leaf cushions or bases 
in larger stem parts. The relatively modest amount 
of secondary xylem in the trunk was high in con- 
ductance (Cichan, 1986). Both aerial and stig- 
marian systems had wood and periderm, but neither 
had secondary phloem (Eggert, 1972; Eggert & 
Kanemoto, 1977). In the Lepidodendron and Lep- 
idophloios-type trees with dichotomous terminal 
branching, the primary body showed progressive 

inution of all tissue zones in the branching sys- 
tem. There were progressively fewer and smaller 
leaf cushions and lamina, and finally maturation of 
apices. Secondary tissues diminished upward from 
the trunk base with periderm extending further up 
than wood. For trees with divided lateral branches 
(Paralycopodites, most Diaphorodendron spe- 
ĉies) the same determinate changes occurred in 
the lateral branches, and presumably in the main 
trunks. Genera that ultimately terminated their 
Pole architecture with at most a fork or two (Para- 


lycopodites, Sigillaria) also are thought to have 
had primary body diminution with apical dichoto- 
my. 

Pole branching was coordinated apparently with 
reproduction (DiMichele & Phillips, 1985). Lateral 
branches were either ramified, bearing multiple 
cones, or pedunculate with a terminal cone. Ter- 
minal “crown type” branching was simply an ex- 
panded scaffolding for cone production and dis- 
persal. Deciduous lateral or terminal branching 
represented minimal interference with light pene- 
tration and wind dispersal. In short, the lepidoden- 
drids represent arborescent reproductive architec- 
ture with determinate growth and fairly short life 
spans (Stewart, 1983: 104; DiMichele & Phillips, 
1985) despite their large statures. 

Lepidodendrids were indeed tree sized with a 
pole design; however, their morphology suggests a 
short life span, perhaps up to 10-15 years as a 
conservatively high estimate. Their developmental 
design indicates a different way to be a big pole 
tree, rapidly and temporarily. Some were also small 
trees, especially in the coal swamps. Most recon- 
structions show lepidodendrids in late to terminal 
reproductive or “death” stages. As unconventional 
as these trees were structurally, and, by inferences, 
functionally, we have good reason to query what 
else is known or likely about their biology, their 
earlier history, and the origin of genera. In turn, 
what do these aspects suggest about the wetland 
environments where they dominated? 

Lepidodendrids were confined to the tropical 
wetlands of the late Paleozoic, apparently associ- 
ated with soft substrates (peat as well as mostly 
inorganic). Such habitats are accordingly referred 
to as coal (peat-forming) or clastic (inorganic sub- 
strates) swamps. There were additional lacustrine 
and fluvial environments associated with the better- 
drained parts of the wet lowlands. Because the 
habitats of lepidodendrids were also depositional 
environments, there is an unparalleled fossil record 
of compressions, casts/molds (see Gastaldo, 1986), 
permineralized anatomy, and spore floras. Rec- 
ognition and interpretation of this record is partic- 
ularly enhanced by the unusually distinctive struc- 
tures of tree-sized lepidodendrids. 

Stigmarian lycopsids represented in our study 
were the principal genera of the Westphalian (mid- 
Late Carboniferous) coal swamps. The five genera 
coexisted in Westphalian (mid-Pennsylvanian) coal 
swamps for about nine million years (Hess & Lip- 
polt, 1986) in an area from Europe to mid-con- 
tinental United States. Additional occurrences in 


